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ON MR KING-HELE’S THEORY 
OF THE EFFECT OF THE EARTH’S OBLATENESS 
ON THE ORBIT OF A CLOSE SATELLITE 


P. F. Message 
(Received 1960 February 8) 


Summary 


The relations are obtained between a and 2, the parameters defining 
King-Hele’s plane II, and the osculating inclination and longitude of node. 
The secular motion of the node is thereby obtained in terms of mean osculating 
elements in a form comparable to that arising in Brouwer’s treatment of the 
canonical equations. A modification is made to King-Hele’s theory which 
removes the singularities in the Je terms of dQ/dys, and it is also found that 
there are no terms of order Je? in « or Q. 





In view of the intrinsic interest of Mr King-Hele’s theory (1) of the motion 
of a close satellite, and since it has played an important part in the study of the 
Earth’s external potential by means of artificial satellite observations, it might 
perhaps be of interest to consider the relations between some of the novel 
parameters that he uses and the more usual osculating orbital elements. 

He defines a plane II, of constant inclination « to the equator, rotating about 
the Earth’s axis in such a way that the satellite moves in it always in the same 
sense about the Earth’s centre. (We shall see that while this can be done exactly 
to the order J, where the second harmonic term in the Earth’s potential is 
(2J )/3r-8P,, a modification is required to order Je.) In the plane II he uses & 
for lougitude measured from the line of greatest slope which lies in the northern 
hemisphere, and r for the radius vector. If Oxyz is a moving frame of reference 
such that O is the centre of the Earth, Oxy is always the plane II, Ox its intersection 
with the plane of the equator, the positive direction being that of the satellite 
when crossing the equator from south to north, and Oy the line =o, then we 
have for the position vector of the satellite, referred to this plane, 


r=(—rsiny, rcosy, 0). (1) 
Then, OX YZ being an inertial frame, such that OXY is the plane of the equator, 


King-Hele’s quantity Q is the longitude of Ox measured backwards (i.e. in 


the opposite sense to %) along the equator. The formulae relating coordinates 
in the two frames are 


X= xcosQ+sin Q(ycosa—zsin«), 
Y= —xsinQ+cos Q(y cosa—zsin«), 
Z= ysina+zcosa. 
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Therefore the position vector is, referred to OX YZ, 
r=(—rsinycosQ+rcosy%cosasin Q, 
rsinyg sin Q+rcos%cosacos Q, 
rcos#sin «). ‘2) 
The velocity of the satellite relative to the frame Oxyz is 
(v)=(—rsing—rbcosy, rcosp—rbsiny, 0), (3) 
and relative to a fixed frame it is 
v=(v)+Oar+O(Ve) (4) 
where @ = (0,0, — (2) referred to OXYZ. (Note that « is constant to order J.) 
Now the osculating orbit plane is the plane defined by the vectors r and v, 
so its unit normal is given by n=N/|N|, where 
N=rav=ra(v)+rO0-r(r. 0) + O(Ve, 0?) (5) 
so that 
N?= {ra (v)}?+27°{r a (v)}. 0+ O(Je, O°). 
Now referred to Oxyz, 
ra (v)=(0,0, 77), 
and referred to OX YZ, 
ra (v)=(—rdsinasinQ, —rssinzcosQ, rbcos«), 
and thus 
N? = rif? — arb Q cos x + Of. Je, 2 Q 


IN| =r pm ee cee +OJe, do (6) 
dis dys 


so that 


n= < —sinasinQ (: + = cosa} —cos sina S (sin ys cos 2 — cosy cos « sin Q), 
dis dys 


IQ 5 i ‘ 
—sinacos 2 (: - = cos 2) + cos sin « — (sin ys sin 2 + cos cos « cos 2), 
au v7) 


d¢ ; ¢ 
+ cos 2( + = cos x) — a + costysintaT} + +O {ye (FG y } ‘ (7) 


Thus in terms of 7, the inclination of the osculating plane, wal KX, the longitude 
of the ascending node, measured from OX, we readily see that 

n=(sinisin {{, —sinzcos {, cos?). (8) 
Comparing these two expressions for n we see that 


q) q 2) 
t=a+ 4D inasin®$ +O Je, & : 
dys dys 


( 
N= —-Q- Gz cossing +O Aue (=) \, 
or, using King-Hele’s result dQ/dj =JL*R* cosz+ O(Je), 
i=a+ }JL?R*sin 2asin? y+ O(Je), 
R= —-Q— LR? cos asin 2ys+ O(Je). (10) 
Now « and Q have no periodic parts of order J, and so these expressions give us 
at once the periodic perturbations of this order in 7 and §Q, and also show that 
the mean value of 7 is given by 
ip =a+}JL*R* sin 2a + O(Je). (11) 
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When the satellite crosses the equator from south to north, the two angles 
and — Q are each equal to the longitude of the satellite measured from OX, and 
therefore they are equal to each other. Thus the mean rates of change of these 
two angles are the same. Now from King-Hele’s equations (88) and (71) I derive 


dQ 2m [/f/R\? . wf/R\ (10... 5 3 iy 
, ir F cos 4 (>) +J (>) (=2sin —s +3o(F (4—7sin?«)}, 
(12) 


T being the mean period between successive returns to the same node. Using 
this, together with the expression for ig above, we obtain for the secular motion 
of the node in terms of the mean inclination i, and the mean radius vector 7, 


f 2 4 4 
ote =— ou COS ly i" (>) +J? (>) (2 sin” ty — 5) + 7 D(5) (4—7sintia, 
(13) 
the J and J? terms of which are in agreement with the result deduced from 
Brouwer’s treatment of the canonical equations for the perturbations (2) (which 
does not include the terms in D). 

The terms of order Je in dQ/d/ as given by King-Hele’s equation (48), 
contain singularities at Y=o and y=. Either one of these, but for general 
values of 8 not both, may be removed by a suitable choice of the disposable 
constant arising in the solution of the equation (47). Now a change dx in the 
value of « must be accompanied, if the value of the colatitude 6 is to be unaltered, 
by a change cosecacoty%d« in Q, which on differentiation with respect to % 
gives an expression with the factor cosec*y, which is therefore associated with 
a change in «. The irreducible singularity in dQ/dy, arising through a factor 
cosec” ys, is a consequence of the fact that we cannot define a rotating plane of 
constant inclination « to the equator such that the satellite always lies in it and 
always, moves in the same sense on it round O unless its maximum northerly and 
southerly latitudes are the same, which they are not if both the eccentricity of 
the orbit and the perturbations are taken into account. Consequently « must 
be allowed to vary to order Je. If we make this modification to King-Hele’s 
theory, then, and at the same time allow for possible terms of order Je? in « and Q, 
we must put 


dQ 
dis 
dx 


ds 


If we follow King-Hele’s procedure, and also put «) for the mean value of «, 
we find that the equations of motion require, to order Je, 


=e(1+eA+2eH+JA) 


=e(ey+eK+JT). (14) 


dy m dy _ Se 
tanya +2A cose (JF 2tany7) = 20s (y+), (15) 


a particular integral of which, having no singularities, is 


A=$.cos (+B), 


y =% sina, sin (f+ 8). 
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The equation for the terms of order Je? is 


rn 


which can be solved by # = K =0, showing that there are no terms of this order 
in « or 2. To the order of J* and D the equation is 


o, (17) 


dH : 
2H + tan ys “ir + 2 tanycosec a). K — cosec a» 
5 


dA dv . 
tan ar) + 2A —cosec %, (5 —2tany.I ) 


+ 


== 2R*L* (: — £ sin? a — — sin? x) sin? ys ) — 3R?L7(2 sin? x sin? ys + cos? x) 
3 3 


| yee — 

+ 4 © RLA3 — 7 sin? % cos? ys), (18) 

7 
which has the solution [' =o, and A as given by King-Hele in his equation (64). 

Thus we have 

% = HX — FE SIN %q COs (ys + B), (19) 

and the only modification to King-Hele’s formulae for Q is in the periodic terms 
of order Je. The mean value of dQ/dt is unchanged, and so this modification 
leads to no change in the estimates of the values of J and D that have been 
derived from observations of this quantity for artificial satellites. 


Gonville and Caius College, 
Cambridge : 
1960 February 5. 
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A SOUTHERN HEMISPHERE SURVEY OF THE RADIANTS OF 
SPORADIC METEORS 


A. A. Weiss and 7. W. Smith 


(Communicated by E. G. Bowen) 


(Received 1960 January 4) 


Summary 


A further survey of sporadic meteor activity, from December 1956 to 
August 1958, has been made with 67 mc/s narrow-beam equipment operated 
at a limiting line density close to 10" electrons/em. Analysis of the echo 
rates suggests that the sporadic radiant distribution may be sharply divided 
into two components. One component is concentrated towards the apex and 
corresponds to a distribution of orbits uniform in heliocentric co-ordinates 
except for a moderate preponderance of direct over retrograde orbits. 
The other component lies close to the plane of the ecliptic, with radiants 
strongly concentrated in two regions centred on longitudes +68° from the 
apex. No change is apparent in the radiant distribution as the limiting 
sensitivity is lowered from the 3rd to the 7th visual magnitude. 





1. Introduction.—Radio echo surveys of the rate of occurrence of sporadic 
meteors over extended periods have been made at Adelaide (Weiss 1957) and at 
Jodrell Bank (Hawkins 1956a). In the former survey the resolving power of the 
aerial was inadequate for the determination of the extent of the concentration of 
sporadic radiants towards the plane of the ecliptic, but groups of radiants near the 
Sun, the anti-Sun and the apex of the Earth’s way were inferred from supplemen- 
tary measurements of the directions of reflection points. The equipment used in 
the Jodrell Bank Survey had much greater resolving power, but the large statistical 
fluctuations resulting from the low echo rate (rarely exceeding 15/hr) tended to 
obscure the form of the diurnal and seasonal variations, with the result that full 
use could not be made of the resolving power of the equipment. Hawkins, by 
treating the whole of the data on a statistical basis, deduced a strong concentration 
of radiant points in the plane of the ecliptic except in the vicinity of the antapex, 
and a large preponderance of direct over retrograde orbits. 

A second survey has now been completed at Adelaide, using a high-resolution 
equipment of the same type as that used by Hawkins, but with a much higher echo 
rate (up to 10o/hr). Although this work is to some extent a repetition, in the 
southern hemisphere, of Hawkins’s survey in the northern hemisphere, an 
important new feature is that the statistical fluctuations in the echo rate are now 
sufficiently small for the diurnal and seasonal variations to be easily recognized. 
The detailed study of these temporal variations leads to a sporadic radiant distri- 
bution which differs in some respects from that obtained by Hawkins. 

2. Collecting area and sensitivity of the equipment.—The equipment, which also 
serves in the detection and measurement of meteor shower radiants, has already been 
described (Weiss 1955, 1960). 

By virtue of the property of specular reflection, whereby an echo can be 
obtained from the ionized meteor trail only if the normal to the trail at the reflection 
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point passes through the observing station, the collecting areas for the two narrow- 
beam, low-elevation aerials are narrow strips on the celestial sphere, lying along 
great circles which intersect near the zenith. Contours of constant sensitivity for 
each aerial are drawn in Fig. 1. ‘These contours have been prepared from the 


South 











North 


Fic. 1.—Contours of constant sensitivity for the Adelaide survey equipment, drawn on the 
celestial sphere. Contours for the N and S channels are identical; only the zero contour is drawn 
for the N channel. Point A represents an arbitrary position of the apex of the Earth’s way, and S 
an arbitrary point source, to illustrate the ys, @, € system of co-ordinates. 

echo rates received from the concentrated Sun and anti-Sun sources, as a function 
of source zenith angle, in conjunction with the variation of the echo rates from 
several showers as their radiants move across the collecting areas. ‘They apply 
to collection in the major lobe only; minor lobes in elevation are excluded by 
range selection of echoes and the minor lobes in azimuth have power gains less than 
4 per cent of the maximum. 

It will be seen from Fig. 1 that the times of appearance of echoes from a discrete 
source in the two channels are a unique function of its declination and right 
ascension ; this fact forms the basis of measurement of radiant coordinates with this 
type of equipment. With the exception of radiants close to the south celestial 
pole, a radiant crosses the S collecting area in a direction roughly normal to its 
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axis, a situation never encountered with the N collecting area. The S collecting 
area alone contains the south celestial pole, consequently radiants with declinations 
less than — 65° do not appear at all in the N channel. 

The line density of the faintest meteor segment detected is approximately1o™ 
electrons/cm, corresponding to a visual magnitude M,~+7. This limit has 
been determined by two methods which give concordant results (Weiss 1960): (a) 
the dependence of the sporadic echo rate on equipment sensitivity; (b) a compari- 
son of the echo rates of the Geminid stream, as detected by the present equipment 
and by the 27 mc/s Adelaide wind equipment, which has been calibrated absolutely. 
The sensitivity was not held constant over the period of the survey, both because of 
severe electrical interference, and of equipment modifications. The corrections 
to the echo rate to take account of these changes in operating conditions did not 
exceed 50 per cent and were usually much less. 

Hawkins (1956a) has given a similar limit of 10" electrons/cm for the Jodrell 
Bank survey, which seems inconsistent with the relative echo rates. However, he 
has subsequently (Hawkins 1956b) revised the Jodrell Bank limiting sensitivity to 
M,= + 42. 

3. The observations.—With the exception of a few short breaks, the observa- 
tions were continuous from December 1956 to December 1957, and again from 
April 1958 to August 1958. No records were obtained from January to March 
1958. Some 3 per cent of the operating time included echoes from accepted 
showers; records for these periods were rejected from the analysis. 





P 1985T January 1957 June 


a 








5 








80} 

















te 20 
Local Time (Hr) 


Fic. 2.—Typical diurnal variations in the echo rate. Each curve is averaged over a calendar 
month, ~——-—-—WN aerial channel; S aerial channel. 
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Mean hourly rates, estimated by averaging the echo rate per }-hr over a period 
of a month, are the basic material used in the following analysis. A first impression 
of the diurnal and seasonal variations in the echo rate may be obtained from the 
diurnal variations shown in Fig. 2 for the four seasons of the year; curves for two 


of the corresponding months in 1958 illustrate the repetition of the main features 
from year to year. 


oe § Channel 


xem—x N Channel 


D J Fr { J A 0 
1957 Month 








Fic. 3.—Seasonal variation in the monthly average of the total daily echo rate. 


S Channel 
----N Channel 


Fcho Rate 


Hourly 


Mean 








Local Time 


Fic. 4.—Diurnal variation in the echo rate, averaged from December 1956 to November 1957. 


Fig. 3 indicates the seasonal variation in the total daily echo rate, corrected for 
equipment sensitivity. The average diurnal variation throughout the year, as 
found for the period December 1956 to November 1957, is sketched in Fig. 4. 








No. 1, 1960 hemisphere survey of the radiants of sporadic meteors 


4. Analysis of the echo rates 

(a) S aerial channel.—In this channel the diurnal echo rate exhibits three 
peaks, each with different seasonal behaviour. ‘The peaks near midnight and 
midday reach maxima in June and December respectively; the peak near 06 hr 
L.T. is a maximum in March and a minimum in September. The echo rate is a 
minimum near 19 hr in all months. On the supposition that this seasonal 
behaviour may be attributed to a sharp division of thesporadic radiant distribution 
into two components, contours of constant hourly echo rate have been prepared 
from the monthly diurnal variations. These are mapped in Fig. 5. Because it 
was known that the equipment sensitivity did not remain constant, the observed 
echo rates for some months have been corrected by small factors (from 0°76 to 
1°50) needed for a reasonable compromise between regularity in the diurnal and 
in the seasonal variations. 





{a) E-Component 











Month. 




















(b) D-Com 


Month 


10 14 22 02 


Locat Time 


Fic. 5.—Contours of constant echo rate from December 1956 to November 1957. Numbers 
against the contours are hourly echo rates. Contours have been prepared from diurnal variations 
similar to those in Fig. 2, on the assumption that the radiant distribution may be sharply differentiated 
into two components. 
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The symmetry of these contours, and also of the average echo rates of Fig. 4, 
about 07 hr (or 19 hr) L.'T. is what would be expected if the perihelia of the orbits 
are uniformly distributed in solar longitude, and if at the same time the ionosphere 
exerts no influence on the detection rate. This confirms the similar conclusion 
reached by Hawkins in his northern hemisphere survey. 

The peak in the echo rate near o7 hr L.T. is due to antapex—apex streaming 
imposed by the motion of ¢he Earth on a component of the sporadic radiant 
distribution which covers the whole of the celestial sphere. We will call this the 
D-component, and for convenience will refer to the apparent concentration round 
the apex as the apex source. The midday and midnight peaks originate in a second 
component, whose radiants are located in two regions with centres on the ecliptic 
at longitudes + 68° from the apex (Section 5). This is the E-component, with 
Sun and anti-Sun sources. The seasonal behaviour of the maximum echo rate 
from these two components is illustrated in Fig. 6. The inference from this 
diagram is that the E-component (Sun and anti-Sun) sources are more highly 
concentrated than the D-component (apex) source. 


Rate 


 E-Component 
e@ Sun Source 
® Anhi-Sun Seurce 


e 
A 


n 





Hourly 


Maximum 


D-Component 
@ Apex Source 





4. i. L i 
20 0 40 50 
Zenith Angle of Source 





Fic. 6.—The dependence of the echo rate from the Sun, anti-Sun and apex sources upon the 
zenith angle of the source. 


The echo rates for the S channel contain no information regarding the width in 
latitude of these sources. From Fig. 1 it is apparent that the axis of the S collecting 
area almost coincides with a line of constant hour angle (approximately + 1 hr), 
and that the time of passage of a point radiant across the collecting area increases 
slowly from 1 hr for = +25° to2 hrford’=—55°. ‘The time of appearance of 
echoes in the 5 channel (relative to transit time) and their duration is thus insensi- 
tive to the declination of any point source which is not too close to the south 








No. 1, 1960 hemisphere survey of the radiants of sporadic meteors 11 


celestial pole. Consequently the time of passage of a diffuse source centred on the 
ecliptic is determined almost completely by the width of the source in longitude. 
Hence the only influence of the gyrations of the ecliptic on the S channel rate is to 
alter the apparent strengths of the sources from month to month; the general 
form of the diurnal variation is unchanged and the three sources are clearly 
resolved at all times of the year. 

(b) N aerial channel.—On the other hand, the peculiar form of the diurnal 
variation on the N channel, in which the anti-Sun source is always resolved but the 
apex and Sun sources are blended over most of the year, can only be understood if 
the D-component has considerable width in latitude whilst the diameters of the 
E-component sources are comparatively small. Lines of constant hour angle 
intersect the N collecting area at angles up to go”, and the length of time over which 
echoes are detected from a diffuse source is now determined by the width in latitude 
as well as in longitude. 

This is most clearly illustrated by the situation at the March equinox. If we 
postulate, for example, a high density of radiants in the D-component at elonga- 
tions from the apex as small as 30°, this component in March will continue to give 
a high echo rate for some 6 hr after the apex itself has passed out of the N collecting 
area. Consequently echoes will be received from the D-component for some time 
after the centre of the Sun source (at =o) has passed through the collecting area, 
and the apex and the Sun sources cannot be resolved. The abrupt fall in the echo 
rate after midday, and the complete resolution of the apex and anti-Sun sources in 
all seasons, point to a close concentration of the E-component to its two centres. 

At the September equinox the zenith angle of the apex source at transit is a 
maximum. The time of passage of the apex source is not so protracted and the 
echo rate is low, and one would expect the N channel rates to be dominated by 
the E-component sources. This also agrees with observation. 

The diffuseness of the D-component also affords a natural explanation of the 
observation by Hawkins in his Jodrell Bank Survey that the echo rate in his sector 1 
increases from 20 to 08 hr L.T., whilst the echo rate in his sector 2 decreases. 
Hawkins’s sector 1 is equivalent to our N collecting area but is not so greatly 
inclined to the N—S meridian (25° compared with our 35°-5); his sector 2 is also 
inclined to this meridian at an angle of 25°, which is larger than the 14° for our S 
collecting area. 

5. The E-component.—From the contours of Fig. 5 the width in R.A. (and 
hence in longitude) of each of the two sources of this component is approximately 
shr. This agrees with the width of 75° found by Hawkins and Prentice (1957) 
for meteors brighter than M,~ +3. It is not possible to determine the width 
in latitude from the echo rates, but evidence from both channels suggests that the 
axes of the sources are comparatively short. Hawkins and Prentice give approxi- 
mately equal widths in latitude and longitude. 

The radiants of the two sources have been determined month by month from 
the times of appearance of echoes in the two channels. Radiant positions in 
ecliptical co-ordinates are illustrated in Fig. 7. The accuracy of the radiants does 
not warrant corrections for zenith attraction and diurnal aberration, even if the 
velocity were known, and they have not been applied. Although an analysis 
based on probability ellipses shows that the distribution of the radiant positions 
about the mean is approximately Gaussian, there is some suggestion of an annual 
variation in all four co-ordinates. It is not known whether this is a real effect or 
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Fic. 7.—Latitude and longitude of the Sun and anti-Sun sources from December 1956 
to August 1958. 


due to errors introduced by the blending of the apex and Sun sources on the N 
aerial channel. The mean radiant positions are: 


Sun source A-A, =-—208+61 B= —2°7+6°7° 
Anti-Sun source A—A,.= +2471 + 6°4° B=+5°5+7°6° 


The centres of the two sources most probably lie on the ecliptic at elongations 
+68° from the apex. This is in excellent agreement with the elongations of 
+ 66° found by Hawkins (1956a). 

The orbits of the meteors constituting the E-component have been described 
by Hawkins. They are probably direct, elliptical with high eccentricity, and with 
semi-major axes comparable with those of the short period meteor streams. 

6. The D-component.—Because the apex source is not resolved on the N 
channel, it is not possible to determine the centre of the radiant area at the time of 
maximum activity in the same way as for the E-component; the identification of 
this point with the apex rests on the results of previous surveys. In confirmation 
of Hawkins’s result, from 15-21 hrs L.T. there is no time displacement between N 
and S channels, which suggests an avoidance of the ecliptic over this period. This 
is consistent with antapex—apex streaming in the D-component, but the following 
qualification is necessary. The seasonal variation in the echo rate near the time of 
passage of the antapex is indicated in Fig. 8, in which is plotted the hourly echo 
rate averaged over the period 18-20hr L.T. Since the celestial pole lies inside the 
S collecting area, this suggests that high declination radiants exist in unusual 
numbers in the early part of the year. This has already been noted (Weiss 1957) 
as a time of relatively low density of meteors in the Earth’s orbit. 
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There are in Section 4 two lines of evidence that the D-component is only 
weakly concentrated tothe apex. Although the present survey affords no evidence, 
it is reasonable to suppose that for this component the radiant density is a function 
of the apparent elongation ¢ from the apex; the visual results of Hawkins and 
Prentice (1957) lend weight to this supposition. This radiant density law 
includes the special cases in which the heliocentric velocities of all meteors are the 
same and the meteor orbits are distributed uniformly in heliocentric co-ordinates. 
These are the models U already introduced in interpreting the results of the previous 
Adelaide survey. Although neither of these assumptions is justified, these models 
U forma convenient starting point for the interpretation of the contours of Fig. 5 (0). 


L 


oe 5 Channel 
xe ex N Channel 








- = oe 
Month 


Fic. 8.—The seasonal variation in the echo rate on the two aerial channels from 18-20 hr L.T., 

around the time of passage of the antapex. 

Computation of complete theoretical contours for comparison with Fig. 5 
is not justified. The extent to which models U fail to account for the observations 
can be judged adequately by confining attention to the echo rates at those times 
when the apex or antapex lies on the axis of the S collecting area (approximately 
o7andighrsL.T.). Echo rates at these times are plotted in Fig. g as a function of 
angular distance 4, where go —y¢ is the elevation of the apex along the great circle 
representing the axis of the S collecting area (see Fig. 1). Additional information 
is available at 13 hr L.'T. in June and December, when the apex—antapex axis is 
approximately perpendicular to the S collecting area; the echo rate is then 16/hr. 

The S collecting area (Fig. 1) is a strip of almost constant angular width centred 
on the great circle perpendicular to the direction of maximum echo sensitivity, 
which is at elevation 10°, azimuth 14°N of E. The behaviour of the echo rate 
across a traverse of this strip is (from shower data) independent of the elevation 
of the traverse. It is then reasonable to approximate this collecting area by a 
strip of constant angular width, with uniform sensitivity along any traverse, and 
sensitivity as a function of elevation go— @ as deduced from Fig. 6. Denote this 
sensitivity by W(@). Let us now consider the times when the apex falls on the 
axis of this collecting area, and denote by D(e) the apparent radiant density along 
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any great circle joining the apex and antapex. D(e) so defined is independent of 
season. Provided the collecting strip is not too wide, the solid angle subtended 
at the centre of the celestial sphere by a small element « to « + de or the collecting 
strip is independent of € except in the vicinity of the apex or antapex where, for 
given de, the solid angle first increases and then decreases to zero as «0 or 180°. 


Modei 
sU tt. 








Fic. 9.—Comparison of observed and predicted echo rates as a function of elevation, go-%, of 
the apex along the axis of the S collecting area. Predicted rates are normalized to the echo rate of 
16/hr observed when the apex—antapex axis is perpendicular to the S collecting area. 

To a sufficient approximation, the echo rate when the apex lies on the axis of the 
S collecting area is 
W(0) D(0—1) a0 


~T, 


N()= | 


where ¢ is a function of season. The apparent radiant density D(e) is given by 


V’ siney deg 

V sine de* 

€y is the true elongation of the radiant whose apparent elongation is e, V is the helio- 
centric velocity and V’ the geocentric velocity of the meteor. We have 


D(e)= 


r 


09 e 
€)=€+ arcsin oo 


dé, V ° Ws Ve 
— =1+[(1-— —=sin’e — cose 
de ( y2 V 

V2? = V2+V,?+2V,V cose, 


V, is the velocity of the Earth. 
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Computed echo rates for V = 42, 35 km/sec are shown in Fig. 9. These have 
been normalized to an echo rate of 16/hr when the apex—antapex axis is perpen- 
dicular to the axis of the collecting strip; the computed echo rate is then 


a] 2 
N,~90 = D(e=90°) : W(6) dé. 
—n/2 
It is evident that the observed echo rates for the D-component may be satisfied 

formally by a uniform heliocentric radiant distribution with hyperbolic velocities ; 
further calculation indicates V~ 50km/sec. This is reminiscent of the hyperbolic 
orbits inferred from visual counts under an identical set of assumptions. Since 
the hyperbolic velocities are no longer acceptable, the observations imply an excess 
of direct over retrograde orbits. The actual ratios of direct to retrograde orbits 
are found to be 2-5 for V = 42, 4:0 for V=40, 15 for V=35 km/sec. These ratios 
are much smaller than that found by Hawkins (1956 a), who gives a ratio of 30 for 
V =40km/sec; Hawkins and Prentice (1957) find a still larger ratio (~ 50) from 
their visual observations. 


— ' 








Local Time 


Fic. 10.—Diurnal variation of echo rate as a function of equipment sensitivity, 1957. Each 
curve ts a ratio of the echo rates measured with a common receiver and two video outputs, one held 
at constant high sensitivity whilst the other varies between medium and low sensitivity A—Oct. 26 
to Nov. 12; B—Nov. 21 to Nov. 27; C—Dec. 2 to Dec. 25; D—Nov. 27 to Dec. 2; E—Novw. 12 
to Now. 16. 


That this discrepancy in the ratio of direct to retrograde orbits is not an effect 
of different limiting sensitivities between the three surveys is established by 
Fig. 10. In this diagram are plotted the ratios (smoothed in a running group of 3) 
of echo counts made on the N channel with a common receiver and two video 
output stages, one held at constant high sensitivity whilst the other was varied 
between medium and low sensitivity. At lowest sensitivity the maximum echo 
rate was 8/hr, and at highest sensitivity, 75/hr. There is a tendency for the echo 
rate at low sensitivity to be comparatively lower near 18 hrs L.T. than near 06 hr, 
which would imply an increase in the proportion of direct orbits amongst the 
fainter meteors. 

7. Comparison with the 27 mc/s survey.—The 27 mc/s survey, extending from 
1952 to 1956, was made with a broad-beam aerial system directed towards the 
zenith. The resolving power was so low that, on the basis: of the echo rates 
alone, it was not possible to distinguish between a distribution of radiants confined 
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to the plane of the ecliptic with an apparent apex—antapex ratio of 4:1 (model £), 
and a distribution uniform in heliocentric co-ordinates, except for a moderate 
preponderance of direct over retrograde orbits. Supplementary measurements of 
the directions of reflection points of trails indicated a concentration of sporadic 
radiants to the apex, with additional strong concentrations near the helion and 
antihelion positions. 

The model E, with radiants closely confined to the plane of the ecliptic, is 
ruled out by the present survey. We then see that the two surveys agree as to the 
main features of the sporadic radiant distribution. In particular, neither survey 
supports a very great preponderance of direct over retrograde orbits in the 
D-component. 

In the present survey the equipment sensitivity was not sufficiently well con- 
trolled for any firm conclusions to be possible regarding the variation of density of 
meteors round the Earth’s orbit. It is however apparent that there is no evidence 
for annual variations much larger than the +20 per cent found in the earlier 
Adelaide survey. 

We have already had occasion to refer more than once to the excellent agreement 
between the present results and the northern hemisphere radio and visual surveys. 
Although some discrepancy remains as to the extent, but not the existence, of the 
excess of direct over retrograde orbits, there is little doubt that the discrepancy is 
not an astronomical effect. ‘There is no major change in the sporadic radiant 
distribution as the threshold brightness is lowered from the 3rd to the 7th 
magnitude. 


Division of Radiophysics, C.S.ILR.O., 
Department of Physics, 
University of Adelaide: 
1959 December 31. 
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Summary 


Observations at the Cambridge Observatories showed that the broad-band 
i” and B photometric systems may be reproduced by observations in narrow 
bands about 200 A broad centred at 5410 A and at 4390A respectively. At the 
Royal Observatory, Cape of Good Hope, measurements were made of the 
radiation flux outside the Earth’s atmosphere in four narrow (~ 200 A) wave- 
length bands from stars of known magnitude and colour index, by comparing 
bright stars with an artificial star, composed of a calibrated lamp and blue 
filter. A 4-inch refractor was used, with interference filters and a d.c. photo- 
meter. The results are given in Table I. They are in good agreement with 
some other recent measurements. Systematic errors are discussed in detail. 

Narrow-band visual magnitudes were measured for 221 stars, and com- 
pared with broad-band visual magnitudes determined by Johnson, by Hogg 
and by Eggen, and with early and recent Cape results. When linear colour 
equations are removed the standard deviations of the residuals are between 
+o™-o914 and +0™-036. The internal random errors of the best individual 
series are believed to be +0™-o10 (narrow-band) +0™-+007 (Johnson) and 
+o™-907 (Cape). It follows that the large differences usually found between 
similar broad-band visual magnitudes cannot be ¢xplained by differences of 
band-width and equivalent wave-length combined with differences of stellar 
spectra. 

There is good agreement between the present narrow-band colour indices 
and gradients on the Greenwich system measured at Mt Stromlo. The 
relation between absolute and Greenwich gradients was found to be 
d=G,,+ 1°11, in excellent agreement with the best previous results. 





1. Introduction.—The total radiation of a star cannot be measured, for the 
Earth’s atmosphere absorbs large parts of the ultra-violet and infra-red radiation, 
and no receiver is perfectly absorbing at all wave-lengths. ‘*‘ Monochromatic ’”’ 
measures in regions only 1 A wide, say, are also impracticable, for the sensitivity 
would be inconveniently low, and the measures would not represent the brightness 
of the star as this is normally understood because stellar spectra do not have a 
smooth black-body energy distribution but are affected by absorption or emission 
lines or bands. Narrow-band measures are an acceptable compromise, for, as 
Woolley (1) has shown, if measures of spectra are restricted to a band-width of 
less than 200A in the visual region the average intensity in the band will be 
proportional to the intensity at its centre, within o-1 per cent. Only the most 
intense stellar spectral features would seriously affect measures made in a region 
200 A wide. 

‘There are now accurate broad-band photometric measures of several hundred 
bright stars, and less accurate measures of many thousands more. It would, 
therefore, be of considerable value to theoretical astrophysicists to determine 
the flux of radiation of certain wave-lengths received, outside the Earth’s 
atmosphere, from stars of apparent magnitude V =o-00 and of various colour 


2 
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indices. The present problem is to select suitable narrow-band regions of the 
spectrum, to measure the intensity of stellar radiation in them, and to find with 
what accuracy narrow-band measures may imitate broad-band measures. The 
observations required were made at the Cambridge Observatories and at the 
Royal Observatory, Cape of Good Hope, and are described below. 

In all absolute determinations of stellar radiation it is necessary to measure 
and to allow for the atmospheric absorption. An important practical advantage 
of narrow-band measurements is that the effective wave-length is defined almost 
entirely by the equipment; it is scarcely affected by the spectral distribution of 
the stars observed. ‘Therefore the atmospheric absorption is not sensibly 
dependent on the colours of these stars. ‘The spectral sensitivity of the detector 
used cannot always be measured with the greatest precision, but if, as here, 
narrow-band measures are made, it can have but little effect on the equivalent 
wave-length of the equipment. Sources with widely different colour temperatures 
may then be compared more reliably. 

2. The selection of suitable narrow-band regions.—\f narrow-band measures 
are made at the equivalent wave-length of the broad-band V system, say, the 
relative brightnesses of stars of various colour indices will be changed very little. 
There will be a quadratic term in the relative colour equation of the two systems, 
due to the band-width of the V system, but this may be smaller than observational 
errors. On the other hand, if the narrow-band measures are made at too short 
a wave-length, red stars will appear relatively faint, blue stars relatively brighter, 
than in the V system. ‘This then provides an observational method of selecting 
suitable narrow-band regions. Ideally, many stars covering a wide range of 
colour index should be observed. 

Observations were made with the 36-inch reflector of the Cambridge 
Observatories to find the required regions of the spectrum. ‘I'wo monochromators 
were used, the first using quartz lenses and prisms, the second a grating and 
mirror in the Ebert—Fastie arrangement (2). Slits were placed in the spectrum 
to transmit the desired wave-length regions. ‘The ‘‘blurring’’ of the edges 
of the regions transmitted caused by the seeing tremor disk or by maladjustment 
of the star image in the focal plane diaphragm was small compared with the width 
of the regions—about 200 A. 

Using the first photometer, narrow-band observations were made of stars 
in certain of the Kapteyn Areas at +15° which had been measured on the 
International System (3), and in the galactic clusters IC 4665 and 4139 measured 
on the U,B,V system (4,5). ‘The second photometer was used to make 
narrow-band observations of stars in the two clusters NGC752 and Praesepe, 
also measured on the U, B, V system (6, 7). 

The centres of regions 200 A wide which gave the best fit to the International 
System were at 4365 A and at 5435A. The first of these does not agree well 
with the value of the effective wave-length of the I, region, 4270 A, determined 
by Seares and Joyner (8). Their value was the weighted mean wave-length to 
which Seed 27 plates were sensitive, the plates being exposed in a small 
spectrograph to ‘‘northern skylight’’. This method gives too much weight 
to the shorter wave-lengths, for even in a milky sky far more blue light will be 
scattered from the atmosphere into the spectrograph than yellow. ‘The better 
agreement for the I,, region between the present result, 5435 A, and the result 
of Seares and Joyner, 5430 A, is readily explained by the much smaller band-width 
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of that region, making their result less sensitive to the spectral distribution of 
the skylight used. The regions giving the best agreement with observations on 
the U, B, V system were centred at 4390 A and 5410A. The ultra-violet region 
cannot be matched by a narrow-band region, except very roughly, because of the 
presence of more intense and more numerous spectral lines and other features 
there. 

The atmosphere at Cambridge is rarely suitable for the most precise 
photometry, particularly if it is desired to compare a calibrated lamp with the 
stars. ‘The observations were therefore continued at the Royal Observatory 
at the Cape of Good Hope. 

3. The measurement of stellar radiation in narrow-band regions.—At the Royal 
Cape Observatory an artificial star, composed of a calibrated lamp and a colour 
temperature conversion filter, was compared with bright stars. A 4-inch refractor 
was used, and a d.c. photometer with an E.M.I. photomultiplier and four 
interference filters, y, 5,, 6,, and 6,, which were chosen using the Cambridge 
results as a guide. y had its peak transmission at 5410A. 6, had its maximum 
transmission at 4390 A, but was relatively transparent at 4340 A and it was feared 
therefore that Hy would affect measures of A-type stars. The two other filters, 











Fic. 1.—The response of the photometer used at the Cape. Secondary maxima in the transmis- 
sions of the ‘b” filters were removed with a dichroic filter, which is responsible for the cut-off at 
41504 and the dip and rise at 4900A to 5000A common to all three. Secondary maxima in the‘ y’ 
filter transmission were avoided by using a photomultiplier insensitive at 7000A, and incorporating 
yellow glass in the filter. 

b,, with maximum transmission at 4600 A, and 5, at 4220 A. were chosen to avoid 
Hf, Hy and Hé. The transmissions of these filters were measured at the Cape 
with a double monochromator and d.c. photometer, the light falling on the 
filters just as when they were in use in the telescope. The relative sensitivity 
of the photomultiplier was obtained from an E.M.I. publication; it is sufficiently 
precise to use typical manufacturers’ data, for the filters alone define the equivalent 
wave-lengths of the photometer rather accurately. Although the sensitivity 
of the photomultiplier falls rapidly near 5400A the equivalent wave-length of 
the y filter is changed only 25 A if the change of sensitivity is neglected. The 
response of the photometer is shown in Fig. 1. Each of the 6 filters was used 
with a dichroic filter to remove secondary maxima. 


2* 
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Observations were made to determine the magnitudes and colours of 16 
‘‘standard stars’’. Each of these was compared with the next at nearly equal 
altitudes on at least 20, and in some cases on 25 nights, in the order a—b-a-b, 
The average standard error for a comparison based on 20 nights was + 0"-003; 
the closing error in the magnitudes was 0™-014+0™-o12. This closing error 
was removed by arbitrarily applying corrections of o™-oo1 to all but two of the 
comparisons. Magnitudes were then adopted for the standard stars, in agreement 
with these comparisons, and adjusting the origin of the scale to agree with that 
of other visual magnitude systems. The colour indices were adjusted to be equal 
to each other and to (B— V) at (B—V)= +0°'5. 

The standard stars were also compared at unequal altitudes to determine the 
extinction. When each standard star was at about 30° altitude in the East it 
was compared with the standard star approaching the meridian, in the order 
a—x—a—x-a. ‘The values of secz for the times of the individual observations 
were obtained from graphs drawn for this purpose, and the zenith extinction 
was computed. Average values for the zenith extinction were: 

y: o™197, b,—y: 0-093, b,—y: o™-122, b,—y: O™-139 
but these must be corrected (Section 4.4) for systematic errors, less than 
o™-o1 in each case. On the total of 55 nights, fewer than 15 per cent of the means 
for the night differed by more than 0™-025 from these values. ‘This includes 
intrinsic scatter, and errors of observation, which were mainly caused by 
scintillation. 

The extra-atmospheric deflection for each standard star when near the 
meridian was then calculated, and compared with the deflection obtained on 
pointing the telescope to the artificial star, to determine the apparent magnitude 
and colour indices of the artificial star. The mean values for these, before 


applying corrections for systematic errors (all less than o™-o1), were 


y: —1™-563, 5,—y: +0™-618, b,—y: +0™599, b3;—y: +0™-691. 
The standard error of each is +0™-003. The observations of the artificial star 
were usually but not always made immediately before or after a determination 
of the extinction, for they had to be fitted into the programme without interfering 
with the equal altitude comparisons of standard stars. In the intervals left by 
these observations 1005 measurements were made of 205 other stars to try to 
discover how accurately narrow-band measures may reproduce broad-band 
measures, and to find the relations between the narrow-band y magnitudes and 
b—y colour indices and broad-band V, B—V measures. 

The flux of radiation of certain wave-lengths received outside the Earth’s 
atmosphere from a star of apparent magnitude V =0-00 was then determined 
from the following measurements : 

(i) the colour temperature and luminous intensity of the lamp, 

(ii) the transmission of the colour temperature conversion filter (see 
Section 4.2, below), 

(iii) the distance of the objective of the telescope from the lamp, 

(iv) the apparent magnitude and colour indices of the artificial star on the 
natural system of the photometer, after corrections for systematic errors 
had been applied, and 
the relations between this natural system and the standard U,B,V 
system of Johnson (9g), defined by the narrow- and broad-band 
magnitudes and colour indices of the stars observed. 
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The results are shown in Table I. The relative values in most of this table 
are believed to have a standard error of +1 per cent, the absolute values 
+2 per cent, but the errors may be as large as +3 per cent for stars of colour 
index B—V greater than +1-2. It will be seen in Fig. 5 below that only four 
very red stars were measured both on the narrow-band system and by Johnson 
on the U,B,V system. These are not sufficient to define the relations between 
these two systems more precisely. 


TaBie I 


The flux of radiation received outside the Earth’s atmosphere from a star of magnitude 
V=o-00 


(Unit = 10~’ ergs cm~? sec (100A)7!) 
(B-V) 53904 4571A 4401A 
4°01 7°06 8-02 
3°99 6°73 7°45 
3°98 6-41 6-91 
3°96 6°12 6°40 
3°94 5°82 
3°92 5°48 
3°90 5°06 
3°89 4°68 
3°87 4°37 
3°85 4°12 
3°33 3°39 
3°82 3°67 
3°80 3°45 
3°78 3°23 
3°76 3°02 
3°75 281 
3°73 2'60 
3°71 2°41 
3°79 2°23 
3°68 2°06 
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4. Systematic errors 


4.1. Calibration of the lamp.—The lamp compared with the stars, known as 
lamp .\, was calibrated in 1957 August:at the National Physical Laboratory, 
Teddington (N.P.L.), in 1958 May at the Council for Scientific and Industrial 
Research Laboratories, Pretoria, South Africa (C.S.1.R.), in 1959 January by 
comparison with other lamps calibrated at the N.P.L., and in 1959 March by 
comparison with a ribbon filament lamp calibrated at Professor Kienle’s 
laboratory at Heidelberg. At the C.S.I.R. Laboratories the luminous intensity 
alone was measured; in every other case both the colour temperature and the 
luminous intensity were measured. ‘The current in the lamp was adjusted by 
measuring the potential drop in a well-aged series resistor of about 0-12 ohms. 
The same series resistance, potentiometer, and calibrating standard cell were 
used at all calibrations and when the lamp was compared with the stars. 

At the N.P.L., lamp A was aged and the current in it was then adjusted so 
that its colour matched a tertiary sub-standard lamp of colour temperature 
2360 °K. The potential drop in the series resistance, 0°0937 volts, was then 
measured with the potentiometer and standard cell. The luminous intensity 
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was also measured, and was found to be 1°88 candelas. 1 candela is defined as 
the luminous intensity of 1/60 cm? projected area of a black-body at the 
temperature of melting platinum, 2042°K. At the C.S.I.R. Laboratories the 
lamp was run as before, and the luminous intensity was found to be 1-go candelas. 
By comparison with other lamps from the N.P.L. in 1959 January, when the 
stellar observations had been completed, the luminous intensity was found to 
be 1-917 candelas. ‘The stellar comparisons do not confirm the small brightening 
of the lamp indicated by these laboratory calibrations. ‘This does not necessarily 
prove an inconsistency, however, for the small difference in question could just 
about be accounted for by possible small secular changes in the behaviour of 
this type of lamp, going beyond the normal tolerances on the laboratory 
calibrations. Comparisons with the ribbon filament lamp yielded an intensity 
3°2 per cent greater than the final N.P.L. result, illustrating the greater difficulty 
of absolute calibrations. It may be noted that the N.P.L. lamps were referred 
to a black-body at the platinum point, and the ribbon filament lamp to a 
black-body at the gold point. 

As noted above, three determinations of the colour temperature of the lamp 
were made, at different times, in different regions of the spectrum, and by 
comparison with lamps calibrated with different black-bodies. ‘The two N.P.L. 
calibrations, made at the beginning and at the end of the stellar observations, gave 
colour temperatures 2360 °K and 2380°K respectively. ‘This apparent rise in 
colour temperature may be due to ageing of the lamp, but is not confirmed by 
all the stellar comparisons. ‘The emissivity of tungsten is not constant, but in 
the visible spectrum rises from long wave-lengths to short, more rapidly from 
7000 A to 5500A than from 5500A to 4000A. The colour temperature of a 
tungsten filament therefore depends on the wave-length at which it is measured. 
The N.P.L. results given above refer to the visible spectrum as a whole, but 
here we are interested only in the region from 5500 A to 4ooo A. Measures of 
the emissivity of tungsten may be used to show that the colour temperature in 
the blue and green region of the spectrum is approximately 15 “K lower than 
in the whole of the visible spectrum (10, 11). The N.P.L. results given above 
are therefore equivalent to results of 2345°K and 2365 °K in the blue and 
green region. ‘The comparison with the ribbon filament lamp, after the stellar 
observations, gave the result 2375 “Kinthe blueand green region. The difference 
of 10°K from the second N.P.L. result is a measure of the difficulty and 
uncertainty of absolute calibrations based on comparisons with different 
black-bodies. The mean of the N.P.L. results for the blue and green region 
was 2355 °K; the adopted colour temperature including the ribbon filament 
comparison was 5 °K higher, viz. 2360 °K in the blue and green region of the 
spectrum. These discordances in the calibration of the lamp, though small, 
are larger than the uncertainty in the comparison of the lamp with the stars. 
In any repetition of this work it would be advisable to make frequent photo- 
electric monochromatic comparisons of the lamp with a black-body source 
at a known temperature, and so avoid many of the present uncertainties. 

The weighted mean luminous intensity of lamp A, at the mean epoch of the 
stellar comparisons, was, for the N.P.L. and C.S.I.R. calibrations alone, 


1/4(1°88) + 1/4(1-g0) + 1/2(1-917) = 1-903 candelas. 


Because a tungsten filament is brighter in the green part of the spectrum than a 
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black-body matched for brightness, as at the N.P.L., over the whole of the visible 
spectrum, this is to be increased by 0-3 per cent when determining the intensity 
at 5390. Comparisons with the ribbon filament lamp indicated a result 
3°2 per cent higher. Equal weight was given to the two measures by increasing 
the N.P.L. result by a further 1-6 per cent. The final value adopted for the 
luminous intensity of lamp A was 1°939 candelas. 

The emission of a black-body at 2360°K was obtained from tables (12) 
based on the following values of the radiation constants : 


¢, = 2hc? = 1-1909 x 10~° erg cm* sec radian ~, 
Cy =hc/k = 1-4380 cm °K, 
where 
cA 
(e°2/4T _ 7) 


B(A)= ergcm~* sec! radian. 

These relative values were then normalized to make the luminous intensity of 
lamp A equal to 1-939 candelas, using the Relative Visibility Factor, or the 
sensitivity of the standard human eye, which is tabulated in various handbooks, 
and the relation between candelas and watts: 680-2 candelas of light of maximum 
visibility (§550.A) are equivalent to 1 watt radian-*. These relations are con- 
sistent with each other and with the definition of the candela as the luminous 
intensity of 1/60 cm? of a black-body at the temperature of freezing platinum, 
2042 “K. 


4.2. Calibration of the filter.—The transmission of the blue colour temperature 
conversion filter was measured at the N.P.L., and at the Royal Cape Observatory, 
by the writer. An ideal colour temperature conversion filter should do nothing 
but change the spectrophotometric gradient of the lamp, a quantity which is 
almost independent of the wave-length at practicable temperatures. The 
logarithm of the transmission of the filter should therefore be a linear function 
of reciprocal wave-length. ‘The best relation of this form which fitted the measures 
was 


logo 7 = — 1°54(2"490— 1/A). 
In Fig. 2 the residuals 


logy 7 + 1°54(2°490 — 1/A) 


are plotted against 1/A to show the deviations of the filter from the ideal, either 
more transparent (+ve) or less (—ve). It will be noted that although the 
agreement between the N.P.L. and Cape measures is excellent from 5300A 
to 6000 A there is an apparently systematic difference between the two sets of 
measures from 4200A to 5200A. However, the greatest difference, at 4800 A, 
is only 0-3 per cent of the light incident on the filter. 

The measurements of the transmission of this filter made at the Cape were 
obtained with a double monochromator and the photoelectric photometer from 
the 4-inch telescope. A Hilger 4-prism spectrograph formed one half of the 
monochromator, with the plateholder replaced by an electric lamp with a small 
filament. Its dispersion ranged from 4 A/mm at 4000A to 20 A/mm at 6000 A. 
The band pass varied between 4 and 20 A. The position of the lamp was measured 
carefully, and the wave-length transmitted was found within +0-4 to +2°0A. 
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A Hilger Constant Deviation Spectrometer was the other half of the double 
monochromator. ‘he photometer was removed bodily from the 4-inch telescope 
and its focal plane diaphragm served as the final slit. 





LOGig T, — LOGig 5 
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Fic. 2.—The transmission of the colour temperature conversion filter relative to the “‘ ideal’”’ 
transmission, logy, T 1°54 (2°490—1/A). 


The difference between the N.P.L. and Cape results might be explained by 
a systematic error of not more than 5 A in the measurement of the wave-lengths. 
The Cape measures were made in two runs, overlapping between 4335 A and 
5142A, the ‘‘camera”’ of the 4-prism spectrograph being readjusted between 
the two runs. An error in the Cape measures of the wave-lengths would show 
as a scatter between 4335 A and 5142A. ‘There is a slight scatter, but this is not 
as great as the difference between the Cape and N.P.L. results. A scale error 
in the photometer used would be expected to be greatest when measuring low 
transmissions, at long wave-lengths, but here the agreement is excellent. ‘Tests 
showed that scattered light was negligible; the apparent transmission of one 
of the interference filters was found at one wave-length to be less than 2 parts 
in 10°, while the maximum transmission was about 20 per cent. ‘The colour 
temperature conversion filter was uniform, and the measurement of different 
areas, both central on the filter, could make no difference. 

In the determination of a smoothed transmission equal weight was therefore 
given toeach measurement. As the Cape measures were not made at wave-lengths 
for which the Planck Function is tabulated, the N.P.L. values were used at first 
in the calculation of the energy from the lamp incident on the telescope objective 
as a function of wave-length. ‘These results were plotted against reciprocal 
wave-length, and a smooth curve was then drawn which bore the same relation 
to these points as did the smoothed transmission function to the N.P.L. values 
of the transmission. ‘The intensity of the energy reaching the telescope was read 
from this graph at intervals of 25 A, and these values were weighted according 
to the relative sensitivity of the photometer with each filter in use. 

4.3. The distance of the lamp from the objective of the telescope-—This was 
determined from a survey of the Royal Observatory grounds (13), and checked 
by direct measurement. ‘The mean of the two results was 111-20m and the 
difference between them 6 cm. 
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4.4. The apparent stellar magnitude and colour indices of the artificial star.—The 
artificial star was compared with a standard star near the zenith once every 
two or three hours. The atmospheric absorption was determined by comparing 
two stars, one near the zenith and the other at an altitude of about 30° in the 
East. Lower observations would give greater ‘‘leverage’’, but to be reliable 
would require the atmosphere to be uniform over a much larger distance, and the 
traces obtained on the Brown Recorder would have shown greater scintillation. 
The extinction was determined from the observations assuming that it was 
dependent linearly on secz, and independent of the azimuth and the colour 
indices of the stars observed. These assumptions are justified respectively if 
the air mass is equal to secz, as it very nearly is up to sec z = 2:0, if the atmosphere 
is well mixed, as it usually is by the south-easterly wind at the Cape, and if, as 
here, narrow-band measurements are made, making the effective wave-length 
of the star-atmosphere-telescope-filter-receiver combination very nearly 
independent of the colours of the stars. 

Four very small corrections have to be made to the apparent magnitude and 
colour indices of the artificial star found in this way, because: 

(i) the zenith distance at which the extinction star was observed was changing 
rapidly ; 

(ii) flexure and magnetic effects may alter the sensitivity of the photometer 
in different directions ; 

(iii) the air in the path of the light from the artificial star absorbed and 

scattered it; and 

(iv) the artificial star faded slowly when unattended. 


(i) The sidereal time of the middle of each observation was recorded to 
o-1 min, and the value of secz for the star was found from graphs drawn for 
this purpose. ‘The zenith distance of one standard star of each pair compared 
was small and changed slowly, but the other, the ‘‘ extinction star’’, was rising 
rapidly; on the average secz decreased by 0-012 per minute. Each colour, 
y, 5, by, bs, was observed for about 50 seconds in that order. The y observations 
were accordingly made 75 seconds before the time of mid-observation, 6, 
25 seconds before, 6, 25 seconds after, and b, 75 seconds after. 

The correction to the y magnitude of the artificial star is small, so it is 
sufficiently accurate to assume that the extinction star was at an altitude of 30°, 
sec = 2°0, at the time of the y observation. The value of secz for the standard 
star was, on the average, 1:08. The true difference of secs would then be 
2:00—1'08=0'92, but the observed difference, that for the time of mid- 
observation, 1°985—1-08=0-go5. The y zenith extinction was therefore not, 
on the average, o™-197 but 


o"*197 Xx aa =0""-1938. 


The difference between the corrected and uncorrected values of the y zenith 
extinction is —0™-0032, and the correction to the y magnitude of the artificial 
star, which was compared with standard stars observed through an average air 
mass of 1°08, is therefore 


1:08 x —0™:0032 = — 0™-0035. 
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Similarly it may be shown that the corrections to the blue magnitudes of the 
artificial star are 

b 


and the corrections to the colour indices are 


1: —0™-0017, b,: +0™-o018, b,: +0™-0059, 


b,—y: +0™-0018, b,—y: +0™-0053, bs—yv: +0"™-0094. 

(ii) Flexure following movement of the telescope and photometer may 
change the optical path or the illuminated area of the photomultiplier cathode, 
and will also change the orientation of the photomultiplier in the Earth’s magnetic 
field. ‘These effects may alter the sensitivity of the photometer and were therefore 
investigated. 

The photometer was attached to the telescope by a welded angle iron frame. 
Flexure was certainly slight, and movement of the Fabry image small compared 
with its own diameter of 3 mm. The average sensitivity of the cathode over an 
area of this size may be assumed to vary almost linearly if the movement is small. 
The photomultiplier and its housing were rotated through 120° after each night 
on which 50 or more observations were made, or after a similar number of 
observations on two nights or more. In this way a comparable number of 
observations was made in each orientation, and the effect of flexure is virtually 
eliminated. 

If the E.M.I. photomultiplier used were perfectly symmetrical about its 
longitudinal axis, rotation about that axis would not alter the effect of the Earth’s 
magnetic field on its sensitivity, but rotation about any other axis might. ‘The 
cell was tested by mounting a radium spot firmly a short distance from the 
photocathode. It was found that the sensitivity varied as the telescope was 
moved to point in turn to the artificial star, the average standard star position 
and the average extinction star position, but the variations were different for the 
three different orientations of the photomultiplier and its housing relative to the 
telescope. The peak-to-peak variation in any orientation of the photomultiplier 
was of the order of o™-o1, but when the average in all three orientations was 
calculated it was found that the sensitivity when pointing to the standard star 
position was only about o™-0007 greater than the average of the sensitivities when 
pointing to the artificial star and to the average extinction star position. Because 
the artificial, standard and extinction stars were observed through air masses of 
approxjmately o, 1, and 2 units, the correction to the magnitude of the artificial 
star is —o™-o014. No correction is required for the colour indices. 

This effect may be explained if the stream of electrons liberated from the 
cathode is not symmetrically placed in the dynode assembly; more electrons 
may be “‘ spilled ’’ off the edges of the dynodes in some orientations than in others. 
The cathode of this photomultiplier is displaced about 2 mm from the mechanical 
centre of the glass envelope. 

(iii) The artificial star was about 111 metres from the telescope. It is 
reasonable to suppose that the absorption and scattering of light by the air in 
this path is the same as that by a similar thickness elsewhere, for the observations 
carrying greatest weight were made in windy, turbulent conditions. At 15 ~C 
the scale height of the atmosphere is about 8-5 kilometres and the extinction of 
the artificial star was approximately 111/8500 of the zenith extinction. The 
corrections are 

y b,-y b,—y b,-y 
— 0™-9026 —o™-o012 —o™-9016 —o-o018. 





No. 1, 1960 Absolute measures of stellar radiation 27 


Ozone contributes about o™-03 to the zenith extinction at 5400 A and is absent 
at ground level, but the corrections are not appreciably affected by this 
refinement. 

‘The observations of Stebbins and Kron at the Lick Observatory (14) over 
an air path five times as long as that used here showed that this simple approach 
is justified by its accuracy. 

(iv) More than an hour after switching on the lamp, when the accumulator 
had reached an almost constant potential and the controlling resistors had reached 
a steady temperature, the current in the lamp fell very slowly. as the accumulator 
discharged. ‘The lamp then faded at a rate between o™-o2 and o”-03 per hour 
if the current was not readjusted. ‘The artificial star was therefore observed as 
soon as possible, about four minutes, after readjusting the current. In this interval 
it would fade by about o™-0017. 

‘The systematic corrections to the magnitude and colour indices of the artificial 
star are summarized in Table II. 


TasLe II 
Summary of corrections to the magnitude and colour indices of the artificial star 
y b-y b.—-y bs—-y 
Observed, uncorrected: — 1°563 +0618 +0°599 +o°691 
Changing zenith distance: —0°0035 +0°0017 +0°0053 +0°0094 
Magnetic effect: —O'0014 ws fe ae 
Absorption of light from the artificial star: —0°0026 —o-0012 —o0'0016 —o-0018 
Fading of the artificial star when unattended: —0-0017_ —0°0003 —0°0003 —0'0003 
Observed, corrected : —1°572 +0'618 +0°602 + 0°698 
5. Comparison with other workers’ results 

5.1. Absolute measures.—The result given in Table I for the radiation 
flux received from a star of apparent magnitude V =o-oo and colour index 
(B—V)=0'60 is 

(3°85 +.0°08) x 10-7 ergscm~* sec! (100 A)~! at 5390 A. 
This may be compared with the result of Woolley and Stibbs (15§) for a star 
of magnitude I, =o-oo and of solar type, 

(3°37 £0°13) x 10~* ergscm~*sec~! (100 A)~! at 5430A. 
Woolley and Stibbs used Kuiper’s (16) value of the apparent magnitude of the 
Sun, I, = —26°84+0°04. If the more recent value of Stebbins and Kron (14) 
is used, I, = — 26°73 + 0°03, and if allowance is made for the difference of zero 
point between the I, and V scales (9, 17), o™-017, the result of Woolley and 
Stibbs for a star of V = 0-00 becomes 

(3°67 + 0°11) x 1077 ergscm~*sec~! (100 A)~! at 5430A. 
The remaining difference of 4:8 per cent cannot fully be explained as the result 
of differences of line absorption in the solar spectrum between 5380A and 
5480 A and in stellar spectra of solar type near 5390 A; such differences are not 
likely to be more than 1 percent. Ifthe colour temperature of the Sun is 6500 °K 
at 5400 A the difference to be expected from the change of intensity of its continuous 
spectrum alone is only 0-7 per cent. 

Dunkelman and Skolnik (18) have measured the solar radiation flux, and 


found 


19°7 microwatts cm~* A-! at 5390A. 
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If this is combined with the writer’s value, the apparent magnitude of the Sun 
is found to be V = — 26-772 or I, = — 26°755, in good agreement with Stebbins 
and Kron’s value, I, = — 26°73. 

5.2. Magnitudes and colour indices.—Narrow-band magnitudes and colour 
indices of 221 stars are given in Table III. The internal agreement of measures 
of stars other than standards shows that the standard error of one observation 
is +o™-o020 in the magnitudes, and about +0™-025 in the colour indices. 


Taare III 


Magnitudes and colour indices of 221 stars 


The stars are identified by their numbers in the Revised Harvard Photometry (H.A.50), 
and by their names. The narrow-band magnitude, y, and colour indices, b—y, are given, 
and the number, N, of observations. Standard stars are indicated by S._ In the last column, 
Notes, stars observed recently at the Cape (20), by Johnson (9), Hogg (21), or by Eggen (22) 
for magnitudes and colour indices, or by S. C. B. Gascoigne at Canberra (24) for relative 
gradients, are indicated by C, J, H, E and G respectively. 
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TasLe [1]—continued 
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Tas.e I11—continued 
HR Name y b,—-y bo—y b,-y 
8911 «x Psc 4°95 +0°175 +0°08 +002 
8969 t Psc 4°18 +0°53 +0°51 +0°52 
8984 A Psc 4°515 +0°295 +0°22 +0°165 
8988 w*® Agr 4°49 +0°145 +0°025 —0°045 
9098 2 Cet 4°55 +0°145 +o°o1 —0'065 


9103 3 Cet 5°045 +1285 + 1°63 + 1°93 5 


The standard errors of the tabulated magnitudes and colours are believed to be 
+o™-o10 and +0™-o13 respectively. ‘They are rounded off to the nearest 
o™-005. Colour indices may usually be measured more accurately than 
magnitudes, but in this case because of the small aperture of the telescope the 
principal source of random errors was scintillation, and the colours, depending 
on comparisons of ratios of deflections, cannot be measured as accurately as the 
magnitudes, which depend on comparisons of single deflections. 

These magnitudes and colour indices have been compared with broad-band 
measures made at the Cape (19, 20), and by Johnson (9), Hogg (21) and Eggen 
(22). The relations between broad- and narrow-band visual magnitudes are 
sensibly linear, but those between colour indices or blue magnitudes are 
non-linear and in some cases have inflections, because of the greater band-width 
of the region measured for B magnitudes and because spectral lines are stronger 
there and more numerous than in the visual region. 

The comparisons of the narrow-band magnitudes and colours with two series 
of measurements made at the Cape, and with Johnson’s measures, are shown in 
Figs. 3 to 7. The V and B—V measures in the first Cape series (19) (Fig. 3) 
were determined from the B,, and C,, measures of Cape Mimeogram No. 1, 
using the relations (20) 

V =B,,—1:31C,,, —0°035 + (0012 —0°0437(C,, —0°52)*) +E 
and 
B—V=1-037C,, + 0°195 

where E is an empirical term for stars with B—V less than —o-1. These 
magnitudes were determined photographically by the Fabry method ; the colours, 
measured photoelectrically in 1949, 1950 and 1951, required corrections for 
non-linearity in the galvanometer, lamp, and scale combination. These measures 
are not regarded at the Cape as a first-class series, and have been superseded 
by magnitudes and colour indices measured photoelectrically with the 18-inch 
reflector at the Royal Cape Observatory (20). This series, known in this paper 
as the second Cape series, is now almost completed. Clearly the agreement with 
the second Cape series (Fig. 4) is much better than with the first. In Fig. 4 
the circle represents y Vel, a star with a Wolf—Rayet spectrum. The only other 
star showing a residual greater than o™-04 from the mean colour equation, 7 Sgr, 
has not yet been fully observed at the Cape. There appears to be no scale error 
(Fig. 6) and the residuals do not appear to increase significantly for the faint 
stars. There are some variations of the residuals with right ascension (Fig. 7) 
but a comparison made eight months ago, when the y magnitudes were completed 
and the Cape magnitudes more than half complete, showed no variation with 
right ascension except that which is now shown between 2" and 8". 
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Fic. 3.—Comparison of narrow-band measures with V and B—V measures derived from Cape 
Mimeogram No. 1. 
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Fic. 4.—Comparison of narrow-band measures with the latest Cape V and B—V measures. 
y Vel, a star with a Wolf—Ravet spectrum, is shown as a circle. 
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Fic. 5.—Comparison of narrow-band measures with V and B—V measures of Johnson. The 
stars shown as circles are «CMa, B Ori, and eleven stars fainter than V=4'8. 
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A preliminary comparison with magnitudes measured by Johnson (9) (Fig. 5) 
was less satisfactory, but a large part of this disagreement was traced to eleven 
stars fainter than V = 4:8, which presumably had not been measured accurately 
by the writer with the 4-inch telescope and narrow-band filters, and two very 
bright stars, « CMa and § Ori, shown as circles in Figs. 5, 6 and 7. The 
agreement between the narrow-band and Johnson’s magnitudes for the 
remaining 65 stars is as good as between the narrow-band and second Cape series 
of magnitudes; the differences between the results appear to be independent 
of the magnitudes of the stars and very nearly independent of their right ascensions. 
The scatter of the residuals appears to be a little larger between 2" and 8"; as 
this is also shown in the comparison with the second Cape series the writer’s 
magnitudes are almost certainly at fault. The errors in the first Cape series 
obscure these small variations. 








JOHNSON 











3 
MAGNITUDE (Vv) 


Fic. 6.—Comparison of narrow-band y magnitudes with V magnitudes measured at the Cape 
and by Fohnson versus magnitude after corrections for colours are applied. The stars shown as circles 
are described in the legends to Figs. 4 and 5. 


Comparisons of the narrow-band measures with the V, B— V results of Hogg 
(21) and the V, P—V results of Eggen (22) show differences similar to those 
found in the comparison with the first Cape series. 

The best linear colour equations between the broad- and narrow-band 
visual magnitudes are given in Table IV. The standard deviations are those 
of the residuals between the narrow-band and other magnitudes after the colour 
equation is removed. It includes all errors of measurement, and the small random 
differences to be expected as a result of the different band-widths and equivalent 
wave-lengths measured and the peculiarities of stellar spectra. The internal 
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Fic. 7.—Comparison of narrow-band y magnitudes with V magnitudes measured at the Cape 
and by Johnson versus right ascension after corrections for colours are applied. The stars shown as 
circles are described in the legends to Figs. 4 and 5. 

standard errors of the latest Cape magnitudes and Johnson’s are believed to be 
+0™-007. ‘The residuals between these and the narrow-band magnitudes are 
very little larger than is to be expected from their internal errors; a further 
external error of +0™-007 is sufficient fully to account for the residuals. If 
such good agreement may be obtained between broad- and narrow-band visual 
magnitudes, slightly different broad-band systems should agree with each other 
at least equally well. In general they do not, but this cannot now be explained 
by the ‘‘slight differences of band-width and equivalent wave-length, and the 
properties of stellar spectra”’ 
Taste IV 


Best linear colour equations relating broad- and narrow-band visual magnitudes 


Source of V Ref. Colour equation Standard deviation 


Cape Series I (19) V=y—o-015 —0'035 (B—V ) +o0™-032 
Cape Series II (20) V =y—0°'024-—0'050(B-—V) +o™-o14 
Johnson (9) V=y—o-015—0'055 (B—V) +0M-923 
Johnson (rejecting 13 stars) +o™-o14 
Hogg (21) V=y—o-015 —o-o15 (B—V) +oM-928 
Eggen (22) V =y—0°'020—0°'070(P—V) +o0™-036 

A comparison of the 5, magnitudes, at equivalent wave-length 4571 A, and 
the narrow-band magnitudes m,,;. of Hogg and Hall (23) shows a colour equation 
of the expected sign and magnitude: 


M4559 — 6; = — 0°262 + 0:023(B—V). 


The zero difference is large because Hogg and Hall adjusted their magnitudes 
to be equal to B or V magnitudes at (B— V)=0-0, while the b, — y colour indices 
were made equal to (B—V) at +0°5. The residuals, + o™-023, are rather larger 
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than would be expected from the stated internal errors of the two series, 
+0™-008 in the m5; measures, +0™-o10 in the y magnitudes and +0™-o13 
in the b,—¥y colour indices from which the b, magnitudes were found. 





-90.! T T T T T T T 
Masso 7b = 0.262 +0.023(8-V) 


T 7 7, 
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Fic. 8.—Comparison of narrow-band b, and m,;,, magnitudes of Hogg and Hall. 


5.3. Gradients.—The b—y colour indices are compared with gradients on 
the Greenwich system determined at Mt Stromlo (24) in Fig. 9. The straight 
lines have the theoretical slope 

AC/AG = > - x; [92 
and are drawn through the mean positions of the groups of points representing 
the hotter stars. Cooler stars deviate systematically from these relations, for the 
effects of spectral lines are greater at short wave-lengths and lower temperatures. 
It will be noticed, for instance, that their effect is much greater in the colour 
index b,—y than in 5, —y. 
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Fic. 9.—Comparison of narrow-band colour indices b—y and relative gradients G of Gascoigne. 


The relation between Greenwich gradients and spectrophotometric gradients 
was determined from the colour temperature of the lamp and the transmission 
of the colour temperature conversion filter, the colour indices of the artificial 
star formed by the lamp and filter, and the Greenwich gradients corresponding 
to these colour indices, obtained from the straight line relations in Fig. 9. 
Greenwich gradient zero corresponds to the following absolute gradients: 


b-—y: +110, b,—-y: +1°16, bs—y: +1°12. 
The second of these results must be rejected, for the b, filter transmits very 
readily at 4340 A and the b,—y colour indices are affected by Hy. The mean 
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of the other two results is in excellent agreement with the best previous results, 
@=G,,+1'11 (25), but the agreement is partly fortuitous, for there is no 
guarantee that the effects of line absorption should be identical at 5390 A, 4571 A 
and 4311 A, in the spectra of hot stars, though these do have but few lines and the 
differences are likely to be small. 
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THE PHOTO-DETACHMENT OF H- 
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Summary 


The photo-detachment cross section of H~ is calculated numerically 
using exchange free waves. The improvement got by using exchange wave 
functions is shown by making a comparison with experimental results. 





The process 
H~-+hv>H+e- 
is of interest because of the réle it plays in determining absorption in stellar 
atmospheres. 

The physical situation is that, initially, the ion is in its ground state, the 
electrons having spins in opposite directions, and, after the transition has taken 
place, the hydrogen atom is in its ground state, the electron spins remaining the 
same. 

This is a problem in quantum theory which can be treated using first order 
time-dependent perturbation theory. This leads to the usual dipole formulae* 
for the absorption coefficient, 


[¥en F2)(2 +22) Py (Fy, F2)d7,d7, 
dipole length 


2 


87 e7y 
c 


eht |p. 9.2 2 
ed | ¥*,(r,, n(x +22 )¥altoraddndr 


2mm*cv 


dipole velocity 
. _ «>... 
[Peston +38 )¥a(ra raddrid 
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dipole acceleration 


x 
2nm?cv® 


where hv = E,, — Ex, E,, Ey being the energies of the system. 
The functions ’ are solutions of the Schrédinger equation 
a: Bice > 
{-8Vit-1Ve- 2-242 —EhW Ge, r,)=0 (2) 
m 1%, "1: 
in atomic units (m=e=h=1). Here the coordinates r, and r, refer to the 
positions of the electrons and the proton is treated as having an infinite mass. 


* These formulae differ from those quoted by Chandrasekhar due to differences in the normaliza- 
tion of ¥’g (see equation (3b) below). 
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Since the electron spins are anti-parallel, the functions ‘’ must be symmetrical, 
ie. Y(r,,r.)=+'¥(r.,r,). Also the orthonormality relations, 


[¥eaern re) P(t, P2)d71,d7,=S, a (3a) 


for bound states and 


[veut r,)' T(r), F2)d7,d7,=5(E, — Ey.) (35) 
for continuous states must hold. 


The function ‘’, corresponds to the initial bound state of the system. The 
final states of the system ‘’, are positive energy solutions of (2) and correspond 
to the system of a hydrogen atom in its ground state plusafree electron. ‘The same 
functions have to be considered in calculating the cross section for elastic scattering 
of electrons by H atoms. 

As yet, it has not been possible to solve the Schrédinger equation (2) exactly ; 
in order to get any solutions certain approximations are necessary. The dipole 
formulae (1) give identical results only for exact solutions. ‘The work done so 
far in this problem shows that the absorption coefficient is sensitive both to the 
functions ‘’ and to the particular dipole formula used. 

Solutions for the H~ ion are got by the usual bound-state variational method. 
The accuracy of the wave function is dependent on the number of parameters used 
and it is most reliable for values of r, and r, giving the largest contributions to the 
energy integral 


I I 


| We(r,, ra) | ae +a hs (ry, ty)dryd ry 


% % Ti 


The usual collision theory treatment of scattering of electrons by hydrogen 
atoms consists in making a partial wave analysis of the wave function as follows: 


(ry Fe) =2 Fi(ha, F2) 
=U 

where 

. I . . . 

Fifa, f2)= {4 (ri)Fi(r2) +4 (radF(es)} 
and ‘(r) corresponds to the normalized ground state wave function of the 
hydrogen atom, 
Y(") = (2). 


The functions F;(r) may be written 


Fi(r)= ¥,,o(F) = filr) 


where Y, ,(#) is the normalized spherical harmonic, 


2 12 
Y;,o(?)= ce *) P,(cos @). 
It can be shown that in tlfe exchange approximation, a full account of which is 
to be found in an earlier paper (1), the functions f,(r) satisfy the equation 
d? [(l+1 , ‘ , 
(ja P+ Ven + Bh fl) = KOON, (4) 


r2 
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where the potential V(r) is given by 


V(r)= a( 1+ set 
and the exchange operator K(r) by 


K(yfie) = aren] — 0+) f _xe-*fi(s)ds+ 2 (3 { “emsul*tf(s)ds 


2l+1 


+ re oi, “Hi(s)as) | ‘ 


The functions f,(r) are normalized to have the asymptotic form 


fir) ~ (3) sin (hr — die -+%). 


Using atomic units this gives positive energy functions satisfying the norma- 
lization condition (36); k? is numerically equal to the kinetic energy of the ejected 
electron in Rydbergs (13-60eV). 

If K(r)= 0 and V(r)=o then the solutions f(r) would be spherical Bessel 
functions and in this case »,=o. This approximation corresponds to using 
plane waves for the final state. The phases obtained when K(r) and V(r) 
are non-zero may be used to determine the elastic scattering cross section. 

The central field equation is got by putting K(r)= 0. Numerical solutions 
of this equation for /=o0, 1, for the energy range 0 < k? < 1°75, have been obtained 
by Chandrasekhar and Breen (2). Solutions of the exchange equations for /=o0, 
1 and 2 have been calculated numerically for the range o< k? <1 on the English 
Electric Computor DEUCE by the present author. A full account of this work 
will be given elsewhere. The equations solved include those for the antisymmetric 
function, obtained from (4) by replacing K(r) by — K(r). These antisymmetric 
solutions have the property 

V("1, 2) = —‘Fi(F2, 61). 

The solutions of the exchange equations will be used in the investigation of 
the free-free transitions of H~ which will be undertaken shortly. The photo- 
detachment and free—-free processes are both of importance for absorption in stellar 
atmospheres. 

The exchange approximation, whilst it allows for the possibility of exchange 
of electrons, otherwise neglects the distortion of the atom by the scattered electron. 
This distortion can be shown to introduce an extra potential, the polarization 
potential, which has the asymptotic form 


2r4 


where « is the polarizability of the atom. At small energies and large distances 
this effect could be important. 

Because of selection rules the only solutions of (4) which contribute to (1) 
are those for /=1. Chandrasekhar (3) has shown that the integrals involved in 
(1) can be simplified to the form 


{ : Wr)filr)dr (s) 
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where W(r) is a simple function and in the dipole length and velocity formulae 
consists of polynomials and exponentials. 

Chandrasekhar (4) in his most recent work used a 20-parameter H~ wave 
function due to Hart and Herzberg together with plane waves. He found that, 
in the dipole velocity formula, there was no appreciable change from his earlier 
work with an 11-parameter function. 

Chandrasekhar and Elbert (5) used this 20-parameter function and solutions 
in the central field approximation for the dipole velocity formula and came to 
the same conclusion. In this work they gave an indication that exchange might 
be important. 

The main part of the present work concerns itself with the use of exchange 
functions in (5) and their effect on the absorption coefficient. The integrations 
were done numerically on DEUCE and the particular cases covered were : 

(1) The dipole velocity formula with 11- and 20-parameter H~ functions. 

(2) The dipole length formula with the 20-parameter function and both 

exchange and central field free functions. 

In Table I the absorption coefficient is given for the particular cases cited 
above including the results of Chandrasekhar and Elbert. The wave length is 
given by the formula 


_ 911-2671 


. k?+] 


A 


I being the electron affinity of the ion in Rydbergs; J takes on the values 0:05512 
and 0-055289 for the 11 and 20-parameter functions respectively. 

Recently Smith and Burch (6) have been able to measure the frequency 
dependence of the photo detachment reaction. Their results were given in 
arbitrary units relative to the value at 5,280 A. They estimate that their measure- 
ments are correct to about 2 per cent. Branscomb and Smith (7) in an earlier 
experiment made an absolute measurement which they estimated to be correct 
to within 10 percent. This absolute measurement was used in putting the relative 
cross-sections on an absolute scale. The best theoretical calculations all lie 
within 10 per cent of this cross-section, and the dipole velocity result with the 
20-parameter H~ wave function and exchange free waves bears the closest 
resemblance in shape to the experimental curve. 

It seems likely that the error in the best theoretical results might be less than 
the 10 per cent error possible in the absolute experimental results. Figs. 1 (a) 
and 1(c) compare the theoretical results (20-parameter H~ function, exchange 
free waves, dipole length and dipole velocity formulae) with the experimental 
results, the absolute values being determined by the theoretical results and the 
curves being fitted at 5,280 A. In Fig. 1 (b) a similar comparison is made with 
Chandrasekhar and Elbert’s work. 

The closest agreement is seen to be in the dipole velocity results of Fig. 1 (a). 
At low energies agreement is not so good; this is most probably due to the neglect 
of polarization in the exchange approximation. The dipole length formula 
gives the poorest agreement, this is probably due to the fact that the largest con- 
tribution to the integrals comes from functions at large distances, i.e. where the 
H~ function is less reliable and where the effect of polarization may be important 
in the free waves. 
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Fic. 1.—Photo detachment cross sections. A comparison of theoretical and experimental results. 
All theoretical calculations use the 20-parameter H~ wave function and the experimental values are 
normalized to agree with those calculated for 5280 A. 

(a) Calculations using the dipole velocity formula and exchange free wave functions. 

(5) Calculations using the dipole velocity formula and central field free wave functions. 

(c) Calculations using the dipole length formula and exchange free wave functions. 
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The favourable agreement with experiment gives an indication of the reliability 
of the free wave functions for determining the elastic scattering cross-section for 
electron collisions with hydrogen atoms. 

The best values for the absorption cross-section are most probably the experi- 
mental values given in Fig. 1 (a), normalized using theoretical results. 
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Summary 


The collision broadening and shift of the line Ca 6573A due to an 
external pressure of helium have been measured. The half-intensity damping 
width of the line was found to be 1:42 + 0°07 x 10~*° cm™! per atom per cc of 
helium, and the shift o-19+0°05 x 10°" cm™'! per atom per cc of helium 
towards the violet. The ratio of broadening to shift is 7-5. This result is 
shown to be consistent with the hypothesis that short-range repulsive forces 
between calcium and helium atoms are predominantly responsible for the 
broadening and shift. No deduction about collision broadening and shift in 
stellar atmospheres can be made from these laboratory results, but it is 
suggested that a study of more highly excited lines may yield astrophysically 
relevant results. 





Introduction.—The experimental study of collision broadening and shift in 
spectral lines is of interest for three main reasons. First, the problem is of some 
importance in physics, and the recent flow of theoretical papers on the subject 
(e.g. 1, 2) indicates that there are still features of the theory which are obscure 
and uncertain. Secondly, with the aid of a tested and reliable theory, information 
concerning interatomic forces may be deduced from measured values of the 
broadening and shift. While it is possible to estimate theoretically the order of 
magnitude of the long-range attractive forces of the van der Waals type, almost 
nothing is known, except for pairs of very simple atoms, of the short-range 
repulsive forces between atoms. As was shown in a previous paper (3), and is 
further illustrated in this paper, these short-range forces may in certain circum- 
stances be an important cause of collision broadening and shift. Thirdly, the 
damping constant of Fraunhofer lines is an important parameter in the theory of 
line formation in stellar atmospheres. It has so far proved impossible to measure 
this quantity directly in the laboratory, because of the difficulty of creating in 
the laboratory an atmosphere of atomic hydrogen at low pressures in imitation 
of the conditions in stellar atmospheres. Collison shifts also may be important 
in some stellar spectra, particularly those of white dwarf stars. To make 
reliable deductions of the stellar damping constant and shift from laboratory 
data an experimentally verified theory of collision broadening and shift is required, 
together with some detailed information concerning interatomic forces. The 
present paper represents a further contribution to this end. 
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Experimental method.—The intercombination line of calcium, 6573, 
45" 1S, — 454p *P, was formed in absorption in an atmosphere of helium, and the 
collision broadening and shift were determined. 

The experimental technique and the apparatus were identical with those used 
previously in the investigation of Ca A4227, 4574S, —4s4p'P, (3). The absorption 
tube had to be heated to about 650 °C in order to give a concentration of calcium 
vapour sufficient to form an absorption line of suitable intensity. Helium 
pressures of 5, 40 and 70cm of mercury were used. The continuous spectrum 
was provided by a tungsten strip-filament lamp, and the necessary high resolving 
power by a Babcock grating mounted in the new Oxford solar spectrograph 
and used in the fourth order. ‘This gave the dispersion of 4:17 mmperA. The 
reference line employed in the measurement of the shift was Hg 4358 A in the 
sixth order of the grating spectrum, and was excited in a mercury 198 isotope lamp. 
This line fell about 10cm distant from the line Ca 6573 on the plate. Ilford R40 
plates were used. These required an exposure time of 15 minutes when the 
spectrograph slit was 0-°077mm wide. ‘The plates were calibrated for intensity 
measurement in the usual way. 4 plates for the measurement of the shift between 
5 and 40 cm mercury and 4 for the shift between 5 and 70cm mercury were 
measured. For the determination of the broadening at 5, 40 and 70cm mercury 
7, 4 and 3 plates respectively were used. 

Measurement of the shift-—The measuring procedure was the same as that 
used in (3). ‘Table I shows the results obtained. ‘The quoted error is in each 
case the standard deviation of the mean. The pressure has been converted into 
particle density with the use of the measured temperature, which was 930 °K 


for the 5 and 40cm mercury plates and go03 “K for the 5 and 70cm mercury plates. 


TABLE | 


Shift between 0°52 x 10'* Shift between 0-53 x 10! 
and 4°16 x 10'* atoms per and 7°49 x 10'* atoms per 
ce of helium 1o~* cm™! ce of helium 1o-* cm“! 
— O'5 +, 39 
+12°0 + 10°3 
+12°0 + 37 
+ &2 +25°9 

Mean +7°9+2°!1 Mean +11°3+4°4 
(—3°4 +0°9 mA) (— 4:9+1°9 mA) 


From these results we find that the shift per helium atom per cc is 
B/N = +0°19+0°05 x 10-*°cm~!cm® 
where the positive sign indicates that the shift is towards the violet. 


Measurement of the damping coefficient.—The line profiles were found from 
microphotometer tracings of the plates, and from these profiles the atomic 
absorption coefficient (apart from a constant factor) was obtained as a function 
of wavelength. The absorption coefficient contours were fitted to Voigt profiles 
at the o-5 and or intensity widths, and the tables of van de Hulst and Reesinck 
(4) used to deduce the damping and Gaussian components of the profiles. The 
individual plate measures are shown in ‘Table II, where the errors are again 
standard deviations of the mean. 


4 
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TasLe II 


5 cm helium pressure 40 cm helium pressure 70 cm helium pressure 
0°53 x 10'8 atoms helium per 4°16 10!%atomshelium per 7°55 x 10'* atoms helium per 
cc cc cc 
Damping Gaussian Damping Gaussian Damping Gaussian 
half-intensity half-intensity half-intensity half-intensity half-intensity half-intensity 
width cm~!_ width cm~? width cm-! width cm— width cm~! width cm~ 
0°047 0076 o°127 0046 0°163 0°072 
o'o71 0°046 0°099 0°060 o'r71 0°049 
0°048 0°065 0098 0°067 0°139 07108 

0060 0°057 0085 0089 
0071 0°044 
0062 0048 
0058 0062 
Means 
0059 0°057 o"102 0066 o'158 0'076 
+0°004 +0°004 +0008 +0°008 +0°008 +0012 


Fig. 1 shows the plot of the damping half-intensity width against number 
density of helium atoms. ‘The equation of the least squares straight line through 
the points is 

Av= 1-42 x 10-°N + 0-048. 











00) 
0.0 2.0 40 6.0 80 10.0 





0 


Fic. 1.—Plot of damping half-intensity width against number density of helium atoms. 
Ordinates are damping half-widths in cm and abscissae are in units of 10%* atoms of helium per cc. 


The intercept 0-048cm~! at N=o is the apparatus function contribution 
to the damping width, so that the collision broadening gives a half intensity width 
of 

2y =1°42+0°07 x 10°-*®Ncm— 
where y is the damping constant of the line in cm~!. The mean of all the values 
of the Gaussian half-intensity width is 0-067 +0-005 cm~!. ‘This is due to the 
Doppler width of the line and the apparatus function of the spectrograph. 
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Discussion of the results —A valuable parameter for the comparison of theory 
and experiment is the ratio 2y/|8|. ‘The above measurements yield the result 


10°6 > 2y/|B| > 5-6 

with the most probable value 2y/|8]=7-5. The value predicted by the 
Lindholm theory for long-range attractive forces whose potential varies as the 
inverse sixth power of the distance between the interacting atoms is 2y/|8|= 2-76. 
Furthermore, it is very probable that the attractive force would give a red shift. 
On the other hand, short-range repulsive forces would most probably give a 
violet shift, and, assuming that the potential of this force varies as the inverse 
twelfth power of distance, the Lindholm theory predicts 2y/|8]=6-70. The 
results lead to the conclusion that, if we accept the Lindholm theory as reliable, 
then in Ca6573A the entire collision broadening and shift are due to the short- 
range repulsive forces. This contrasts with the case of Ca4227A formed in a 
helium atmosphere, where short-range and long-range forces were found to be 
of comparable importance, and produced a very small resultant shift. 

No information can be deduced from these measurements about collision 
broadening and shift in the Fraunhofer line A6573. It is certain that long-range 
forces are more important when the perturbing atoms are hydrogen than when 
they are helium, but how important they are relative to the short-range forces 
it is impossible to say. In order to gain useful information for astrophysical 
purposes it is necessary to find spectral lines for which long-range forces are the 
predominant cause of broadening in a helium atmosphere, for then the use of the 
laboratory results for a helium atmosphere to predict the damping constant and 
collision shift for an atomic hydrogen atmosphere can be put on a sound theoretical 
basis. It is probable that such lines will be the more highly excited lines of, 
e.g., the calcium spectrum, and it is proposed to investigate these in the future. 


University Observatory, 
Oxford : 
1960 April 1. 
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THE MAGNESIUM b LINES IN LATE-TYPE STARS 


T. 7. Deeming 
(Communicated by the Director of the Cambmndge Observatories) 


(Received 1960 February 1) 


Summary 


The absorption lines of Mg t at A 5167~-84A have been measured photo- 
electrically in 539 Go-Ks stars. ‘The results show an increasing absorption 
with advancing spectral type and a negative absolute luminosity effect. No 
significant differences are found between the strong and weak line groups of 
Miss Roman, or between weak-CN and strong-CN stars. High velocity 
class III stars have slightly stronger Mg than normal. It is suggested that 
the Mg b lines can be used as a good criterion of the NIK luminosity class 
for stars with B—I|° between o-7o and 1°30, and as a fairly good criterion 


of MK spectral tvpe between GS and K5 





1. Much less work has been done on stellar spectra at wave-lengths above 
5000 A than in the range 3500-5000 A and the only previous work of any extent 
on the magnesium b triplet at AA 5167-3, 5172°7, 5183-6 is that of Thackeray 
(1,2), who measured by photographic photometry equivalent widths and 
central intensities of a number of lines, including A 5183-6, in 8o late-type stars. 
He seems to have been the first to give a substantial quantitative account of the 
negative absolute magnitude effect shown by this line, although this had been 
reported briefly by Ohman (3) a few years before Thackeray's first paper appeared. 

The three-channel spectrophotometer at the coudé focus of the 36-inch 
reflector of the Cambridge Observatories which has been described by Griffin 
and Redman (4), offers an admirable method for the observation of almost any 
strong feature in a stellar spectrum for wave-lengths from about 4000 to 6000 A. 
Here it has been used for the more extensive measurement of the Mg b lines 
with a view to their possible use as criteria of luminosity, spectral type or 
population type in late-type stars. 

The Mg b lines can be detected in main sequence stars as early as Ao but this 


programme is restricted to stars in the range Go-Ks5, where the lines are rather 


strong and easily measurable. 

2. The spectrometer and the technique of observation and reduction are 
exactly as described by Grithn and Redman (4). Integrated intensities have 
been measured in each of three adjacent regions of the spectrum with wave-length 
limits 

A 5110°4-5139°8 A, 

B  5160-2—518g°0 A, 

C §214°4-5243'0A. 
These wave-lengths were checked during each night’s observation and kept 
constant to better than +1 A, but it was found experimentally that errors in 
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wave-length adjustment of up to 5A did not materially affect the results. The 
effects of stellar radial velocity on the results were therefore assumed to be 
negligible. 

The magnesium b triplet lies in the central channel, and if A, B and C 
respectively are the integrated intensities in each of the channels then the ratio 
(A+ (C)/B forms a measure of the absorption in the central channel referred to 
the mean of the other two. The excess of absorption in the central channel is 
mostly due to the magnesium lines. Following Griffin and Redman we use 
not A, B and C themselves but A/A,, B/B,, C/C,, where A,, B,, C, are 
measures made with a lamp, interpolated to the time of the star observation. 
The *‘ magnesium ratio”’ r is (A/A, +C/C,)+B/B,. The advantages of using 
this particular combination of measurements have been fully discussed by 
Griffin and Redman and all their arguments apply equally to this programme. 

539 stars with spectra between Go and Ks have been measured and the ratio r 
formed for each star. ‘These results are quoted in Table I which shows 
(1) Henry Draper number, (2) Mg b ratio r. The value quoted is generally the 
mean of six separate observations taken in pairs on three different nights. The 
r.m.s. error of one of these observations is + 0-05 and of the mean about + 0-02, 
i.e. one per cent. (3) Adopted spectral type, mostly by Miss Roman (5, 6) on 
the MK system; otherwise from the General Catalogue of Radial Velocities (13). 
(4) Colour index B—V by Johnson et al. (7, 8, 9, 25). (5) Total space velocity 
in km/sec relative to the local standard of rest, from Miss Roman’s papers. 
(6) ‘The notes indicate members of the strong lines (s) and weak line (w) groups 
of stars of Miss Roman. Giants for which Griffin and Redman find a CN ratio 


more than 0:08 above normal for the appropriate spectral type and luminosity 
class are called strong CN stars and are indicated by CN+. Giants with a 
CN ratio more than 0-05 below normal are called weak CN stars and are indicated 
by CN-. 


TABLE I 


HD 
166 
417 
443 
613 
1400 
2774 K2 Ill 
2901 K2 Ill 
2925 Ko II] 
3346 Ks Ill 
3457 K4 Ill 
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N WN 


MN NNN ON 


to 


to 


4656 





T. #. Deeming Vol. 121 


TABLE I—continued 


HD Spectrum (B-V) 
4817 : cKs5 

5234 : K2 Ill 
5286 . sgK1 

5395 P G8 III-IV 
5516 G8 III-IV 
5916 G8 III-IV 
6186 Ko III 
6497 K2 III 
6582 Gs VI 
6833 Kir Ill 
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12533 
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cKs 
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TaBLeE I—continued 


HD Spectrum (B-V) 
18449 K2 III 
18970 2 Ko II-III 
19476 Ko III 
19656 Kr III 
19787 Ka2 Ill 
19735 Ks III 
20123 Gs Il 
20468 K2 II 
20277 sgG8 
20618 sgG5 


20644 K3 II-III 
20893 K3 Ill 
21120 G8 III 
21552 K3 III 
21754 Ko II-III 
22072 dG7 
23183 . Ko III 
23841 K2 III 
25604 , Ko III 
25893 ° dK2 


25975 , K1 Ill 
26311 ° cKs 
27022 , Gs Ill 
27371 ° Ko III 
42° 939 ; dKs 
27348 . G8 III 
27382 : Kr Ill 
27697 . Ko III 
27971 “4 K1 III 
28100 . G8 III 


28292 . K2 Ill 
28305 . Ko III 
28307 
29038 
29139 
29317 
30504 
30834 
31421 
31767 


31782 
32923 
33554 
33556 
34334 
34255 
34559 
34575 
35186 
35620 
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TABLE I—continued 


HD Spectrum (B-V) 
37160 ° G8 IIIp 0°95 
37171 e Ks III 1°55 
37394 ‘ dK 0°84 
37601 : sgGo9 

37984 ‘2 Ki Ill 1°16 
37981 ? Kr IV 

38230 ‘68 K1 V 0°83 
38656 i G8 III 0°94 
38751 "4 G8 III 
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TABLE 1—continued 


HD Spectrum (B-V) 
58207 Ko III 

58367 . G8 III I°OI 
58972 . K3 III 

59148 K2 Ill 

59294 “4 K2 Ill 

60318 Ko Ill 

60522 Mo III 

62044 Ki Ill 

62285 : Ks Ill 
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G6 III 
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TABLE Il—continued 


HD 
82741 
82885 
83240 
53425 
83805 
34453 
85503 
87837 
$9269 
89484 


= 


Spectrum (B-V) Notes 
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G8 III w 

Ko IV 

K2 III CN+ 
K4 III 

dGs 

Ko III 


> 


wWNUMMwWBN SD Ww 
-=~IwWww Duis ss] Ww 


NNN NNN NW DN WN 


Ki III 
sgKo 

G8 III-IV 
sgKo 

G8 II-III 
K3 Ill 

Ka III 
K3 III 

K2 III 
dGo 


90250 
90572 
99537 
QgIOoIr 
g1612 
92095 
92424 
92523 
94054 
94132 


NY NNWR 
wr Ww 
ouu 


~ . 
7333 
w ° 


N N 
> ww > 


“NM WN 


© 


Ki III 
K3 III 
Kr III 
dG8& 
K2 III 
K1 III 
Ko III 
sgG7 
Ki III 
G7 IV 


94264 
94247 
94600 
94549 
94669 
95345 
95659 
96436 
96833 


97561 


N NNN NN KN DN WN 


97907 
98262 
98524 
99196 
99491 
99648 
100470 
100696 
101454 
IOI501 


K3 III 

K3 III 

Kir III-IV 
K4 III 
dKo 

G8 II-III 
Ko III 
Ko III 

Ko III 
G8 V 


- Dw 
oS oOo 


onriu N 


NNN NN WN N ND N 


w NNW N 


ba | 


101673 
102224 
102328 
103095 
104556 
104979 
105043 K2 Ill 
105475 dGog 

105631 , dK 

105963 ° dK2 


K3 III 
Ko III 
K3 Ill 
G8 VI 
G8 V 

G8 III 


w 


Oxy 


nw 


wu Ap = Nm 
au 


NNN NN N NON 





No. 1, 1960 The magnesium b lines in late-type stars 


TaBLE 1—continued 


HD Spectrum (B-V) 
106714 Ko III 
106760 K1 III 
107328 , Ko III 
107383 G8 III 
107469 dKo 
107950 : G7 Ill 
10822 G8 III-IV 
108381 Ki III-IV 
109317 : Ko III 
110833 dKo 


111028 , Ki IV 
111067 47 K3 III 
112033 "22 Gs III 
112989 , Kip 

113092 : G8 III 
113226 . Gog III 
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HD 
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139341 
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TABLE I—continued 


Spectrum (B-V) 
G8 III 0°96 
K4 III 

Ko III 1-02 
K3 III 

K2 Ill ey, 
K4 III “40 
Ks III 

Kop 

dK4 

G8 IV 
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Ks III 
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TABLE I—continued 


HD Spectrum (B-V) 
155410 K3 III 
156283 K3 Il 
156681 . K4 II-III 
157999 : K3 Il 
158633 dK1 
158899 K4 Ill 
159966 Ko III 
160346 dK3 
161074 K4 III 
161096 K2 III 
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162211 
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TABLE I1—continued 


HD Spectrum (B-V) Notes 
182488 ; Ko V 0°82 
182572 . G8 IV 0-78 
182762 , Ko III 
184010 sgG8 
184398 K2 II-III 
184467 Ki V 
184406 K3 Il 
185194 G8 III 
185351 Ko III 
185622 cK5 
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TABLE I—continued 


HD Spectrum (B-V) 
199191 . G8 III 
199169 . K4 III 
199253 . Ko III 
199580 . K1 IV 
200905 . Ks Ib 
201051 , sgK1 
201091 . Ks5 V 
201092 2: K7 V 
201251 . K4 II 
202109 . G3 II 


203344 Ko III-IV 
203504 Kir Ill 
204771 Ko III 
205435 G8 III 
205512 Kir III 

dG8 

Ko III 
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Gs Ib 


Ko Ib 
207134 , K3 III 
209747 K4 III 
209960 K4 III 
209945 “ Ks III 
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TaBLe [—continued 
ND 
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Spectrum (B-V) Notes 
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3. Results 


3.1. Variation of r with spectral type and colour index.—Fig. 1 gives the mean 
run of the Mg ratio 7 with spectral type for each of the four principal luminosity 
classes, dwarfs, subgiants, giants and supergiants. The MK luminosity classes 
II-III, and III-IV have been included in the giants and the class IV-V has 
been included in the dwarfs, so that only fairly well attested subgiants have 
been treated separately. All luminosity classes show a marked increase of r 
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Fic. 1.—The mean run of r with spectral type for each of the principal luminosity classes 
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with advancing spectral type. The r.m.s. scatter from these mean curves is 
about 0-06 per star, except for dwarfs, where it is about o-12. The chief cause 
of the greater spread for dwarfs appears to be less accurate spectral classifications. 
For dwarfs also there are signs of a systematic difference between MK types 
and others, in the sense that the MK type tends to be a little earlier. It should be 
noted that if the class V curve in Fig. 1 uses only stars with MK types a smoother 
run is obtained, although only 23 stars are then available from G8 to Ks. Further, 
in that case the curve lies higher, the difference in r being about o-o2 at G8, 
0:05 from Ko to K2, and perhaps o-15 at K4. The evidence therefore is that the 
use of MK types gives a better separation of class V from class III stars. 

The colour index (B—V) is in some respects a more satisfactory parameter 
than spectral type and it is important to find the behaviour of r with (B—V) 
for those stars for which the colour index is known. ‘There are 158 such stars 
on the programme. Fig. 2(a) shows r plotted against (B—V) for stars of all 
luminosity classes, and in order that a fair comparison should be made, Fig. 2 (b) 
contains a diagram with the same stars as appear in 2 (a), but with spectral type 
as abscissa. It should be emphasized that this diagram contains less than 
one-third of the total number of stars used to obtain the curves of Fig. 1. 

In the case of luminosity class III stars there is a linear relation between r 
and (B—V), within the errors imposed by the scatter of individual star values: 
For Class III, r=0-622(B—V)+1°677, 0°85<(B—V)<1°55. 

The r.m.s. scatter in r per star about this line is 0:08, which is slightly greater 
than the scatter for the same luminosity class in the (r, spectral type)-relation. 

There are also approximately linear relations between r and (B—V) for 
class IV and class V stars, although there are not enough data to determine them 
well. In the case of class V stars there is a suggestion of a maximum for r in 
the vicinity of (B— V’) = 1-0, but there are not enough points to give any certainty 
of this. We may note that there is possibly a maximum at K4 or Ks in the 
(r, spectral type)-relation, although again this needs much more data for 
confirmation. 

3.2. Variation of r with absolute magnitude.—The discovery by Thackeray 
of the negative absolute magnitude effect shown by the Mg b lines is confirmed 
by the present measurements. ‘The increased absorption in dwarfs can be seen 
in Figs. 1 and 2. If we anticipate the later calibration of our ratio scale in terms 
of Thackeray’s equivalent widths, of A 5183-6, the only line of the triplet which 
he measured, we find that the increase in equivalent width passing from a typical 
Kol II star to, say, 54 Psc, a Ko V star, is 0-3 A, or about 25 per cent. 

In the range of surface temperatures covered by our measures, and within 
the range of absolute magnitude covered by classes III to V, the b lines increase 
in strength both with decreasing temperature and with decreasing intrinsic 
luminosity of the star. At constant spectral type class V stars are hotter than 
those of class III, but at constant (B— V) the two are of approximately equal 
temperature, so that we might expect a plot of r against (B— V) to show a greater 
separation of dwarfs from giants than a plot of r against spectral type. This is 
indeed what is observed, ‘The change in equivalent width from a class III star 
of colour (B— V)= 1-00 to a class V star of the same colour is 1-8 A, an increase 
of over 100 per cent. This assumes that we are allowed to extrapolate the linear 
relationship between r and equivalent width beyond the domain of overlap 
of Thackeray’s measurements with ours. 
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Fic. 2.—(a) Mg ratio r as a function of (B—V). 


(b) The same stars as appear in (a) but with spectral type as abscissa. Crosses indicate 
supergiants; dots, giants; filled circles, subgiants; and open circles, dwarfs. 
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Fig. 1 shows that subgiants have absorption intermediate between that of 
giants and dwarfs and that supergiants have the same or somewhat weaker 
absorption than giants. We note that when the abscissa is colour (Fig. 2 (a)) 
there seems to be a better separation of giants and supergiants, but there are then 
only four supergiants available. Until we pass K3 the giant and dwarf curves 
in Fig. 1 (spectrum as abscissa) run more or less parallel, but the subgiants 
appear to be nearer to dwarfs in the late G types and to giants in the early K’s. 
And in Fig. 2 (a) (colour as abscissa) a somewhat similar trend is shown, with 
subgiants merging with dwarfs near (B—V)=0-°6. It should be pointed out 
that we have Mg ratios for only 39 subgiants, that MK types are available for 
about half of these only, and (B— V) values for even fewer. Lack of data is clearly 
the chief difficulty here. 

3.3. The Mg ratio r as a criterion of spectral type and of absolute magnitude.— 
We have seen that r is a function of spectral type and of absolute magnitude so 
that we can use it either to determine spectral type if the luminosity class is known, 
or to determine luminosity if we know the spectral type or colour. ‘The latter 
problem is generally the more important. 

The purely observational r.m.s. error of an Mg ratio from six measurements 
is +002. For luminosity class III this would produce an error of about 0-04 
in the deduced spectral class or o™-03 in (B—V). But the scatter arising from 
measurement errors is small relative to that arising from the intrinsic scatter of 
individual stars, which corresponds to about 0-12 of a spectral class (perhaps 
rather more for dwarfs), or o™-12 in (B—V). 

The error of determination of a spectral type, provided the luminosity 
class is known, compares favourably with the errors of the MK system itself or 
with the accuracy obtainable from analogous measures of the G band (4) and 
further, except possibly in luminosity class V, the b lines avoid the ambiguity 
in the G band type arising from the maximum near K2, so that they could with 
profit be incorporated into a scheme of spectral classification, at least as far as 
Ks. In general there does not seem much useful purpose in using the b line 
measures to deduce a (B—V) colour. One needs prior information of the 
luminosity class and the resulting colour determination is not particularly precise 
relative to the errors obtained in straightforward good photoelectric photometry. 
For instance Morgan and Johnson (7) claim a probable error of +0™-oog in 
(B—V). Given the climate it would appear preferable to measure (B—V) 
directly. 

We do not have enough MK spectral types for the dwarfs, nor do we 
have enough measures of (B—V) to obtain comparable results for dwarfs, but 
there seems to be no a priori reason why we should not expect similar accuracy 
for the determination of spectral types of class V stars also, at least to K2 or K3. 

The possibility of using r as a luminosity criterion is more attractive, especially 
if we know the colour of the star concerned. It is evident from Fig. 2 (a) that 
if we know the colour of a star in the range (B— V)=0-70 to 1-30 then the Mg 
ratio will provide an almost unequivocal determination of luminosity class. 
Almost, because we are not yet certain of the behaviour of subgiants in the diagram, 
and because the separation of points on the supergiant side of the class III stars 
is not good. (In part this may merely be an indication that the supergiants do 
not form a well separated sequence of their own and that there are many 
intermediate stars which have generally been given luminosity class II on the 
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MAK system.) ‘lhe very large scatter in the later-type dwarfs make prediction of 
luminosity rather uncertain if only spectral type is known. However the 
separation when colour is used instead is so good that it would seem better to 
estimate spectral type from a measured (B—V) and to use the colour with a 
measured ratio r to give a luminosity class, rather than to estimate spectral type 
directly and try to obtain a luminosity classification on that basis alone with no 
direct colour measurement. 

It is unfortunate that not enough good absolute magnitudes exist for stars 
in general for us to determine directly whether the scatter in our measurements 
can be ascribed to a real spread in the absolute magnitudes of the class II] stars 
in the solar neighbourhood, a ‘‘ cosmic scatter’’. ‘The use of spectroscopically 
determined absolute magnitudes to resolve this question is not permissible as 
there may be some effects which cause changes in line intensities without there 
being any real changes in absolute magnitude, and such effects may influence 
both the Mg b lines and the lines customarily used for absolute magnitude 
determination, so that a spurious correlation might be obtained. It is found 
that there is no correlation between the Mg ratio (with the effects of spectral type 
and luminosity class removed) and the spectroscopic absolute magnitudes of 
Young and Harper (10), although this cannot be taken as conclusive evidence 
that the scatter in x is not due to scatter in absolute magnitude. Further measures 
of ry in a number of galactic clusters may help to elucidate the *‘ cosmic scatter ’”’ 
problem. 

It seems then that if we have a measurement of (B— V) in the range 0-7 to 1°3, 
the Mg b ratio provides a good criterion of luminosity class, but that for the 
moment we are not justified in supposing that it will give accurate determination 
of absolute magnitude. If we know spectral type but not colour index we can 
still make a separation of luminosity classes but the distinction is in this case 
not quite so clear. 

3.4. Population type differences, high velocity effects—N. G. Roman in her 
paper on the spectra of F5—Ks5 stars (§) discusses certain groups of stars which 
she says are characterized respectively by weak or strong Fraunhofer lines, the 
two groups tending to have high and low velocities respectively. Griffin and 
Kedman have pointed out that for stars near Ko the designation ‘* weak-line’”’ 
in fact means strong A 4227 Cai and G band, ‘“‘strong line’’ means weak A 4227 
and G band. Miss Roman also discusses stars having unusually weak or strong 
CN, giving them respectively the designations ‘‘weak CN ”’ or “‘4150’’. Since 
Griffin and Redman’s measures appear to give more complete information on 
this latter point, for present purposes we shall take ‘‘weak CN”’ to refer to 
giant stars with CN anomalies less than —o-o5 and ‘‘strong CN ”’ to refer to 
stars with CN anomalies greater than +0-08. In ‘lable I ‘‘ weak line’’ stars 
have been marked w, ‘‘strong line’’ s, ‘‘weak CN’’ CN-, ‘“‘strong CN” 
stars CN+. ‘The differences from the mean for the appropriate spectral class 
have been taken for the class Il] stars which belong to one or other of these 
groups. ‘The mean value of the ‘‘ Mg anomaly’’ for each group is given in 
Table Il. ‘The r.m.s. scatter in each case is such that we cannot regard any 
of these groups as having an ‘‘ Mg anomaly”’ very significantly ditterent from 
zero. ‘There is a hint of very slightly weak Mg among ‘“‘strong line’”’ stars, 
to parallel the weak A 4227 of Cai, and a suggestion of very slightly stronger Mg 
for stars having weak CN. 





No. 1, 1960 The magnesium b lines in late-type stars 


Tasle II 
Mean anomalies for weak line, strong line, weak CN and strong CN stars 


Stars Mean anomaly r.m.s. scatter per star 
weak line -0°003 + 0°006 0°049 
strong line —0°016+0°009 0066 
weak CN +0°014 +0°009 0064 
strong CN +0:005 +0°010 0'059 


lig. 3(a) shows a plot of the Mg anomaly, evaluated using spectral type, 
against space velocity (relative to the local standard of rest) according to 
Miss Roman. ‘There is no sign that the Mg ratio is systematically less for high 
velocity stars; on the contrary there is a slight tendency for high velocity stars 
to show greater r (67 high velocity stars of luminosity class III give an ‘‘ Mg 
anomaly’? +0°034+0-009). Such a trend would be consistent with the idea 
that high velocity stars show more of the characteristics of under-luminous stars. 
In Fig. 3(6) the Mg anomalies have been taken with respect to (B—V). 
Unfortunately this restricts the number of stars considerably, but again there is 
a slight tendency for a greater r at higher velocity. 

To summarize, there are no certain variations of the Mg ratio with other 
spectral characteristics usually denoted by ‘‘weak line’’, ‘‘strong line’’, 
‘*weak CN’’, ‘‘ 4150’; there is possibly a slight increase with space velocity. 

4. Correlation of Mg ratio with Thackeray’s 1939 measurements.—Our 
measurements cover a much larger wave-length range than did Thackeray’s, 
and inevitably we take in absorption lines of other elements, chiefly Fe, and in 
the later types we include a certain amount of molecular absorption by TiO, 
MgH, C,. The two sets of measurements are therefore comparable only to a 
limited extent, but since the line 5183-6 which Thackeray measured is 
attributable almost solely to Mg1, we may use his measurements to check whether 
our measures are affected to any considerable degree by these other lines or bands. 

We have 55 stars in common with Thackeray and Fig. 4 gives a comparison 
of his measures with ours. The correlation appears to be satisfactory and supports 
the view that what we have been measuring is in fact primarily Mg1 absorption. 
The regression line can be represented by 


W = (3°26 + 0°03)r — (5°88 + 0°07) 


where r is the Mg ratio and IV is the equivalent width of 4 5183-6 in Angstrom 
units. This diagram omits two values (for « Tau and 5 Lac) for which Thackeray 
introduced corrections for ‘TiO absorption. This relation may be found useful 
if our Mg ratios are to be compared with the predictions of a model stellar 
atmosphere. 


5. Theoretical interpretations of Mg intensities. Until now we have not been 
concerned with any reasons for the effects discussed. The analysis has been 
on an empirical basis and the results derived for the variations of r with 
luminosity or with spectral peculiarities of various kinds will remain valid 
whatever the physics is behind the ‘‘ Mg ratio’’. In what follows we shall 
be concerned with the interpretation of r and the physical basis of the variation 
of r in stars of different types and luminosities. 
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Fic. 3.—(a) Mg anomalies, evaluated using spectral type, plotted against space velocity as 
computed by Miss Roman. 


(6) The same, evaluated using B—V colours. 
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5.1. First we wish to know whether the absorption that has been measured 
is primarily due to magnesium or whether there is an appreciable contribution 
due to molecular absorption by TiO which is known to occur strongly from K5 
onwards, or from MgH which has been reported in K stars by Ohman (3) and 
by Davis (12). Code (11) reports that in population II K giants there exists 
a dip in the continuum near the region we have measured which is not present 
in late-type dwarfs. ‘The origin of this dip is not assigned and Code does not 
make it clear whether this is a feature peculiar to population II spectra. Further- 
more the diagrams of continua which he publishes are for two stars of strongly 
differing spectral types, HD 103095, G8 VI and HD 128902, K4III, so that 
one cannot be sure what the result of comparing a giant and a dwarf of the same 
spectral type would be. 

It happens that we have measured both these stars, finding an Mg ratio 2°43 
for HD 103095 and 2-49 for HD 128902. We would expect this intensity dip 
to tend to increase the Mg ratio: in fact HD 103095 has a somewhat lower Mg 
ratio than we should expect of its spectral type and HD 128902 whose Mg ratio 
should perhaps be somewhat increased, shows an unusually low ratio for a 
K4III star. Presumably this dip is due to some molecular band absorption, 
probably TiO; in any case TiO bands should show in the tracing of HD 128902, 
if they are going to affect the Mg ratio appreciably. Since we find little if any 
effect of this dip, we conclude that it, and molecular absorption, have not affected 
the Mg measurements to any appreciable extent. This is confirmed by the 
correlation we find between r and Thackeray’s equivalent width. (The point 
at W=2-7A, r=2:58 in Fig. 4 is 6And and although it seems to suggest that 
the Mg ratio is too high for the equivalent width, this particular value for W is 
rather uncertain according to Thackeray.) 
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5.2. Assuming that r represents a measurement of the behaviour of the Mg 
b lines, we now wish to interpret this behaviour in terms of the physical properties 
of the atmosphere of the various types and magnitudes of stars which we have 
measured. The behaviour of Mg 5183-6 with spectral type has been predicted 
by Aller (13) and by Russell (14) on the basis of the Boltzmann and Saha formulae 
together with some allowance for variation of continuous absorption. Russell 
predicts a maximum at GgV and Aller one at 7=4250°K. If we interpret 
this temperature as the reversing layer temperature and follow Allen ((15), p. 189), 
the spectral type would be K2 V-K3V. The only suggestion of a maximum 
that the present measurements show is at (B—V)=1-10 or at K4V but this 
is based on only two stars, 61 Cygni A and Bat (B— lV’) = 1-19 and 1-38 respectively. 
61 Cygni B is the only dwarf in ‘l'able I with an MK type later than K5s—-Morgan 
classifies it as K7V. Ohman (3) gives a reproduction of a Mount Wilson 
spectrogram of this star, together with a microphotometer tracing showing the 
(0,0) band of MgH. ‘The spectrum is considerably depressed to the short 
wave-length side of A 5200 and this would perhaps have the effect of reducing the 
ratio r by reducing the light in channel A. ‘Thus we should not put too much 
weight on to 61 Cyg B in deciding whether or not there is a maximum. Certainly 
we can say that there is no sign of a maximum at Gg as Russell predicts and any 
maximum is unlikely to be as early as K3 as Aller predicts, although the star 
with the highest measured 7, HR 753, has an MK type of K3 V. ‘The theoretical 
difficulty is to know how the continuous absorption coefficient varies along the 
spectral sequence and what part, if any, molecular absorption plays in 
determining it. 


c.2 


5.3. Asimilar difficulty is encountered in trying to explain why the Mg b lines 


are stronger in dwarfs than in giants. The reduction in electron pressure in the 
giants and the consequent greater ionization produces a change in the right 
direction but appears to be several orders of magnitude too small. From a detailed 
comparison of the spectra of « Boo and 7oOphA, Miss 5S. A. van Dijke (16) 
deduces that lines of excitation potentials less than about 2-7 volts are stronger 
in giants while lines of E.P. greater than 2-7 volts are stronger in dwarfs. ‘The 
actual dividing potential depends on the element and its state of ionization, but 
is in general in the range 2°5-3:0 volts. The excitation potential of the lower 
states of the Mg b triplet is 2-7 volts which is just on the dividing line, so that 
we cannot easily say whether the Mg b lines are in agreement with the tendency 
indicated by Miss van Dijke’s measurements. The lines of Mgi which she 
measured, with E.P. 4:3 volts were found to be stronger in the dwarf. Aller and 
Stoddard (17) find for Mg 3838 a similar variation with spectral type and absolute 
magnitude to that which we find for the Mg b lines, and point out that this 
disagrees with the theory of Pannekoek (18). The line 3838 belongs to the 
triplet AA 3829-4, 3832-3, 3838-3 which arises from the same lower level as do 
the b lines (19), but the other triplets from this level are too far in the ultra-violet 
to be observed successfully. Nevertheless the similar behaviour of A 3838 and 
A 5183 suggests that the negative absolute magnitude effect arises principally 
from differences in the populations of the 3P® level, or from differences in. the 
continuous absorption, rather than from some anomaly in the region of the b lines. 
possibly due to TiO or MgH. The Cat triplet analogous to the b lines is at 
ANG162°2, 6122°2, 6102-2. ‘These lines are reported by Miss Burwell (20) to 
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be stronger in dwarfs than in giants, the differences increasing with advancing 
spectral type. No data appear to exist on the triplets Sri AA7070°1, 6878-4, 
67g1°1, and Bat AA7g05°8, 7392°4, 7195°2. 

As was pointed out by Ohman (21), it is known from studies at eclipses 
that Ca and Mg appear very high up in the solar chromosphere (e.g. S. A. Mitchell 
(22) gives 6000 km for Mg 3838, 3500 km for Mg 5184, 1000 km for Ca 6162). 
In addition Christie and Wilson (23) in a study of the 1934 eclipse of ¢ Aurigae 
studied the behaviour of the Mg triplet AA 3829, 3832, 3838, and found that the 
behaviour with height in the atmosphere of the K star is markedly different 
from that of other lines, showing a minimum at 0-05, and a maximum about 
o-08R,. (pp. 448 to 453). Miss van Dijke (16) considers that there is evidence 
for a stratification of elements in the atmosphere of « Boo, similar to that 
observed in the Sun, although the data on this point do not include magnesium. 
\ll this suggests that the difference of Mg b intensities between giants and dwarfs, 
and the similar behaviour of the 3830 A Mg1 triplet and the 6130 A Cat triplet, 
may possibly be related to peculiarities in the distribution with respect to height 
in the stellar atmosphere of the atoms concerned. A similar idea was put forward 
by Lindblad (24) many years ago. 

5.4. The foregoing discussion has been restricted to giz unts and dwarfs 
because for subgiants and supergiants there are insufficient data to draw any 
certain conclusions. The peculiar position of « Peg at (B— V) = 1°58, r= 2-32 in 
Fig. 2(a) suggests that perhaps interstellar reddening could be an fates 
factor for supergiants. However according to Morgan, Harris and Johnson (8) 
e Peg is not noticeably reddened with respect to the other supergiants g Peg and 
é Cyg. 

The behaviour of subgiants is more interesting. Although there are 
relatively few stars available for discussion, there is a tendency for them in earlier 
types to resemble dwarfs, as far as Mg absorption is concerned, and in the later 
types to resemble the giants. This is borne out by Fig. 2 (a), where although 
there are only 11 subgiants, they form a well defined sequence. ‘This behaviour 
is perhaps what would be expected on the basis of current ideas on evolution, 
as the later-type subgiants are supposed to be among the oldest stars in the solar 
neighbourhood and are presumably structurally similar to giants (A well-known 
example is 6 Eri which unfortunately was not measured because it was too far 
south to be accessible with the 36-inch telescope at Cambridge.) ‘lhe early-type 
subgiants cannot long have left the main sequence, so that they should still bear 
a similarity to dwarfs. 

‘The increase of Mg absorption in high velocity stars is not well enough 
established to merit extensive comment here. We would merely point out that 
if lines in the blue—violet part of a spectrum are weak then there is a tendency 
for the assigned spectral type to be too early and the measured colour will be 
too blue for the true effective temperature of the star. In this case the Mg 
absorption would appear to be too high for the measured (B — V) or spectral type. 

6. Stars with peculiar Mg b lines.—There are a few stars which stand out 
from the general run of r with spectral type and (B— V’) in Figs. 1 and 2. We give 
notes on some of these, but we would emphasize that this list is not exhaustive. 
Note that space velocities are based on assumed absolute magnitudes which may 
need correction. 
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(i) HD 6833 

Kr III, (B— V) =1-18, CN ratio 2-03, Mg ratio 2-20, space velocity 235 km/sec. 
The weak Mg at first suggests a supergiant but this is ruled out by the weak CN. 
It may be that we have the beginnings of an indication of metal deficiency in an 
extremely high velocity star. 

(11) HD 29038 

K3 III, (B— V)=1-18, CN ratio 2-03, Mg ratio 2:65, space velocity 101 km/sec. 
A class III star of this spectral type and colour would normally have a much 
stronger CN and much weaker Mg. A class V star would have this CN but 
stronger Mg. Class IV stars are rare in spectral types as late as K3. 

(iii) HD 31782 

Ko lll, (B—V)=0-80, CN 1-93:, Mg 2:40. The colour, CN, and Mg all 
suggest that this star is a dwarf and not a giant. 

(iv) 54 Her (HD 152879) 

K4IIl, (B—V)=1-41, CN 2:23, G band 2°86, Mg 2-25, space velocity 
74 km/sec. The CN and G band are normal. The Mg is very low both for 
this spectral type and for this colour. 

(v) BHer (HD 148856) 

G8III, (B—V)=o0-91, CN 2-23, G band 2-71, Mg 2:06. The colour and 
G band are normal, the CN is slightly strong for the spectral type, the Mg ratio 
has the lowest value measured for any star. 

7. Suggestions for further work.—I\t would be helpful to measure the Mg ratio 
in a greater number of stars, particularly in subgiants, which could easily be the 
subject of a large investigation, designed to clarify their position in relation to 
class III and V stars. It would also be of interest to measure the analogous red 
calcium and strontium triplets to see whether they behave in the same way as 
the Mg b lines. 

The results show that as a luminosity criterion the Mg ratio is much more 
effective when used in conjunction with colours than with spectra. There are, 
unfortunately, only a comparatively small number of stars in the sky for which 
(B—V) colours are available. Accurate colour work cannot profitably be 
attempted in this country but in a good climate the data could be obtained fairly 
rapidly for comparatively bright stars such as are considered here, without 
elaborate apparatus. 
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IN ‘THE PLANETARY NEBULA NGC 7027 
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Summary 


The abundances of O°, O+ and O*? are determined from [O1], [O11] 
and [O 111] forbidden lines and the abundances of O**,,O** and O*° from 
Our, Otv and Ov recombination lines. In the recombination calculations 
cascade from all upper levels is taken into account. It is found that O** 
is the most abundant oxygen ion and that the abundance of O*° is small. 


The chemical abundance ratio 1s found to be N(O))N(H)= 1-2 x 10 


1. Introduction. Since gaseous nebulae are optically thin in the observed 
emission lines the numbers of emitting atoms or ions are proportional to the 
observed line intensities. In this respect the problem of determining the 
chemical composition of a nebula is simpler than that of determining the 
composition of a stellar atmosphere. ‘The main difficulties in the nebular 
problem are, first, that there are large departures from thermodynamic 
equilibrium and a great deal of atomic data is therefore required in order to 
interpret observed line intensities and, secondly, that for many chemical 
elements the ionic abundances can be deduced from observations for only one 
or two stages of ionization and approximate estimates must be made for all 
other stages (1, 2, 3, 4). 

Aller, Bowen and Minkowski (5) have made a particularly detailed study 
of the line intensities in the spectrum of the planetary nebula NGC 7027. 
Atomic oxygen is represented by lines of [O01], [O11], On, [O1nrj, Orn, Ov, 
and Ov. ‘Vhe [O1], [O11] and [O11] forbidden lines are excited by electron 
impact. Using calculated collision cross-sections and transition probabilities 
the abundances of O°, O* and O*? may be deduced from the forbidden line 
intensities. Most of the permitted lines are excited by radiative recombination 
but the Bowen fluorescent mechanism (6) (resonant absorption of Het Lyz 
in the transition O** 2p?°P,--2p3d°P,, possibly due to a chance coincidence in 
wave-length) excites the strong Ou lines. The main work of the present paper 
is the calculation of recombination spectra for Out, Orv and Ov and the 
deduction of abundances for O+?, O+t4- and O+°%. Our results determine the 
total oxygen abundance without any significant correction for unobserved stages 
of ionization and provide information about the ionization equilibrium which 
is useful for the determination of abundances of other elements. 

For the calculation of recombination spectra one must consider not only 
direct captures on a given level but also captures on all higher levels followed 
by cascade. Calculations of reasonable accuracy have been made previously 
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only for Hi, Het and Hert (see (7)). Im order to be able to make such 
calculations for complex ions we have developed an approximate general method 
for the calculation of recombination rates (8) and have obtained improved methods 
of performing the cascade calculations (9). 

2. General remarks on the calculation of recombination spectra.—Let the energy 
levels be denoted by subscripts i=1, 2, 3,... in order of increasing excitation 
energy. We consider the usual cases, Case A for an optically thin nebula and 
Case B for a nebula optically thick for lines ending on the i=1 ground state. 
Let the radiative transition probabilities be denoted by A;;. For Case A the 
probability that population of 7 is followed by a direct radiative transition to j is 
P,, = Aj;/A; (2.1) 


where 


(Case A). (2.2) 
k=1 
For Case B transitions to 7=1 are counterbalanced by absorptions from i=1 
and the effective value of P,; is therefore P;;=A,;/A; where now 
i-1 
A;= > Ax (Case B). (2.3) 
k=2 
Let C;; be the probability that population of 7 is followed by a transition to j 
when all cascade routes are taken into account and let C;,=1. ‘Then (g) 


i-t 
C= > PiCai- (2.4) 
k=j 
Let N,N .a, be the number of captures on k per unit volume per unit time, 
N, being the electron density and N, the density of recombining ions. ‘The 
number arriving on 7 due to capture and cascade is 
co 


NN Cui (2.5) 


k=i 
and of these a fraction P,; give the 7-»+j line. The number of quanta per unit 
volume per unit time emitted in the line is 


9= NN Hi; (2.6) 
where the effective recombination coefficient is 
@ 
Xi; a Py 2 Cri (2.7) 
=4 


For Cases A and B we use the notation x, x;. It sometimes happens 
that some of the z,% are much greater than the x; this is the case if 


¢—1 
Ay > > Ax (2.8) 
k=2 


or if a similar inequality holds for upper states which are important for populating 
level i by cascade. If «,{") is only a little larger than «;(” then the resonant line 
quanta, k->1, need be scattered only a small number of times to effect the 
conversion from Case A to Case B but when «,;“ is much greater than «,@) many 
more scatterings are required for the Case B limit to be approached. 

At first sight it would seem that the question of whether Case A or Case B 
applies could be settled by estimating the optical depths for resonant line scattering, 
but that the problem is more complicated may be demonstrated by considering 
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the efficiencies of fluorescent mechanisms. ‘The essential ideas concerning the 
Out fluorescent mechanism may be understood by considering O*? to be an 
atom with three levels, i=1, 2 and 3. Heit Lyx happens to have a wave-length 
practically identical with that of the 3-1 line. Absorption of a single He 11 Lyx 
quantum may be followed by re-emission of 3->1 radiation or by emission of 
observable 3-2 radiation. It will be shown (Section 7.4) that these two processes 
have probabilities P,, = 0-983 and P;,=0-017. Let W (Het, Lyx) be the number 
of Heit Lyx quanta produced in the nebula and let (3,2) be the number of 
3>2 quanta escaping from the nebula. Define the overall efficiency of the 
fluorescent mechanism as 
= A (3; 2) (2 ) 
~ AN (Hen Lyx)’ ~ 
If each Hei! Lyx quantum can escape from the nebula after being scattered n 
times by O*? it is readily shown that 
& =(1— P*,,) (2.10) 
which approaches unity for m large. ‘The condition that # should be close to 
unity is essentially the same as the condition that Case B should apply for the 
recombination spectrum. 
For NGC 7027 the O11! resonant line optical depth is 7,, ~~ 10' (Section 7.1). 
If scattering by O** were alone considered we would have n~7,,7~ 10° and 
hence A=1. Allowance for scattering by He* would cause n to be larger still. 
But from the observed spectrum .# may be deduced and we obtain #7 =0-43 
(Section 7.4). Using a slightly different method, it is shown elsewhere (10) 
that A~o-3 for both the Out and Nu fluorescent efficiencies in a number of 
nebulae. Putting # = 0-43 in (2.10) gives n= 33. ‘The large discrepancy between 
n~ 10° and n= 33 must be due to some physical processes having been overlooked. 
The following might be considered: (i) absorption of Ol resonance quanta 
in ionizing H atoms; (ii) the Zanstra effect (11) of scattering with redistribution 
of frequency; (iii) the wave-length difference of o-o1g A between the centres 
of the He11 Ly« line and the O11! absorption line. 
3. Energy levels, transition probabilities and recombination coefficients 
3.1. Energy levels.—Let E; be the energy of level ¢ measured in Rydbergs 
(13°60eV or 109737 cm~') and let EF; be negative for bound states and zero 
at the series limit (foo). ‘lhe energy parameter e; is defined by 
E; =27¢, (3.1) 
where z is the residual charge; s=m+1 fortheion X*”. ‘The effective quantum 
number »; is defined by 


I 
E=-s (3.2) 
V; 


and the quantum defect yu, by 
y= 1 — ; 
where n; is the principal quantum number. Fora given spectral series the quantum 
defects may be expressed as a series in €;: 
bh; = a+ be; + ce? 4 eves (3-3) 
Assuming LS coupling, the recombination coefficients «(ySL/) are pro- 
portional to (2/+1) and the level populations V(ySLJ) determined by 
recombination and cascade will also be proportional to (2/ +1). We may 
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therefore neglect J quantum numbers in the recombination calculations and 
consider only centres of gravity of spectral terms: 
d(2J + 1)E(ySL/) 


MySL)= 2 (3- 
at nia (28+1)(2L+1) +) 


3-2. Transition probabilities—TVhe quantal formula for A;; is 
4, 
are 
A, = {| —) 2'a~,; 5 
ij 6a, ij (3 5) 


4n 
(¢ ‘) = 2°6774 x 10° sec! (3.6) 





where 


6a, 


(x being the fine structure constant, c the speed of light and a, the Bohr radius) 
and where 


iy = (€;— 62? — (3-7) 
S;; being the line strength (12) in atomic units and «, the statistical weight of 
level 7. For a single jumping electron 
Ai, =(€,—€, PC(t-7)|2(nl\|r|njl;)P. (3.8) 
The coefficients C(i-j), which have been discussed in a previous paper (8), 
have the properties: 
(i) wjC(t->j) = wjC(j>1) (3-9) 
(ii) For the transitions O* * 1s?2snl->1s*2sn'l’ and O** 1s?2s"nl->18728?n'l', 
C(L+l’)=1, /(214+1) (3.10) 
where /., is the larger of /,/’. 
(ii) For the transitions O+? 2pnlSL->2pn'l'SL’ the coefficients C(plL—pl'L’) 
are given in ‘lable I. ‘hey satisfy the sum rule 
> C(pILpl'L')=C(Lel’). (3.11) 
LV’ 
(iv) For 2pnlSL—2p*SL’, 


C(plL—p*L') =2C(plL-—ppL’). (3.12) 


TaBLe I 


C(pILp—l'L’) 

1/9 
1/3 
5/9 
2/3 
1/6 
1/2 
1/150 
1/10 

14/25 


Ne Nee Ne OD 


w 


Most of the integrals (n,/,|r|n;/;) were obtained from the Bates and Damgaard 
tables (13); in this approximation 2(m,/,|r|n,/;) depends on the effective quantum 
numbers v;, v; but does not depend explicitly on z. For a few transitions lying 
outside the range of the Bates and Damgaard tables the tables of Burgess and 
Seaton (8) were used; for [v;—v,| not small these two sets of tables give results 
which are practically identical. 
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For transitions to 2p* it is important to include the term involving ¢= 1 + du/év 
in the normalization of the 2p radial function (14). When this factor is included 
it has been shown (8) that good results are obtained for photo-ionization from 
2p? and equally good results may be expected for bound-bound transitions 
2p?—2pnl. In practice this means that the values of ./(2p*--2pn/) obtained 
from the Bates and Damgaard tables must be divided by (14) 

C=(v—-1)(v+2) v(v+T) 
where vis the effective quantum number for the 2p electron in the 2p? configuration. 
3.3. Recombination coefficients. —For recombination to level /, 


. 
> Ww; 


x . (=) : (mkT’) 1S e fuk? | | (/i 7a,(i je hukT d( hp) (3.14) 
c" \7 I 


WwW, 
ba i 


where w. is the statistical weight of the recombining ion, J,the threshold ionization 
energy and a,(@) the cross-section for photo-ionization by radiation of frequency 7. 
The photo-ionization cross-section is given by (8) 


praes’ (21+ a ae ( . ‘] A 2 
yp AZ Vib, € bb pe I 
3 \ / 2 Gul \ Prauii 


where /; is in Rvdbergs, z*e’ is the ejected electron energy in Rydbergs, 

Cry =C(tet tL’) (3.16) 
where 7, specities the state of the ion, 7’ the angular momentum of the ejected 
electron and L’ the total angular momentum of ion plus ejected electron. ‘The 
quantity ,’,,, is the extrapolated quantum defect for the 7,n'l’L’ series. ‘The 
g(vl, el’, ) are obtained from the tables of Burgess and Seaton (8). 


From (3.14) and (3.15), 


slitnta 2r\ gy, 2 
, (=) “i = \12)(A) (3.17) 


> 


4 L 


where (27! 2a‘a,7c/3) = 2: < 10° cm sec}, A= (ARes*/ kT) = 157 890 37/ T and 
DA=VS Cop | (1+ v2 PleQilve ere POO de. 
] l’ ~ OV 
For photo-ionization from 2pn/L the final state is 2pe'l/’L’ with l'=14 
I / 
and L’=/, l’+1. It was found that the amount of calculation could be 
substantially reduced without significant loss of accuracy on replacing p’,,, by 
>¢ VVEVL 


es (3.19) 


Using the relations (3.9), (3.11) we then have 
s ( ‘ : lea Eel )? Cil -] eg (ule Ue’, )|*. (3.20) 


L 
Using the Burgess and Seaton expression for g, one obtains the expansion 


> Ca 


where 
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4. The O11 recombination spectrum.—The O*? ion has ground configuration 
1s*2s"2p? with spectral terms *P, 'D and !S. Assuming all O*% ions to be in 
the ground state, 1s?2s*2p, recombination will populate O+? states 2pnl4Z and 
2pnl*L. We require the effective recombination coefficients for the O11 lines 
belonging to the transitions 2p3p->2p3s. 

Table II gives the O*? quantum defects «(2zpn/SL) obtained from observed 
energy levels (15). ‘This includes the constants a and 6 in the linear extrapolation 


formulae »=a+be. Extrapolated values of u, given in brackets, are used when 


the energy level has not been observed. States n/ with />3 may be taken to be 


hydrogenic. 
Tase I] 
O*? quantum defects p(2pnl.SL) 


Singlets 


p 


A 


og 








0°45! 
0°425 
(0°412) 
(0°407) 


0°395 
—0°362 


Triplets 


s p 


A. 
‘ 2 


3p 3S 3p 8p 3p 





0°506 


0°627 0°393 0°368 0°424 0°049 


0°371 
(0°362) 
(0°357) 


o°610 
0604 
(0-600) 


0°347 
—0°'316 


0°594 
—o'186 


0°308 
(0-285) 
(0°273) 
0°247 
—o'840 


o°401 
(0-392) 
(0-387) 
0°377 
—0°313 


0°049 
(0°049) 
(0-049) 


0°049 
0°000 


Table III gives the elements C;, 5,5, of the cascade matrices. For Case B 
triplets we exclude transitions to 2p**P and for Case B singlets we exclude 
transitions to 2p?'D but include transitions to 2p?'5S. 


Table IV gives values of D,(A). For the triplets the integrals in (3.18) were 
evaluated numerically for temperatures of 1 x 104 and 2x10'°K. Only one 
rather weak singlet line is observed and the D,(A) for the singlets were therefore 


calculated a little less accurately, for 1 x 104°K, using the expansion (3.22). 
6 
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Taste III 
The cascade matrix Cj; 
Singlets 
Case A Case B 


A = * 
a, 


3p 'P 3p 'D 3p 'S 3p 'P 3p 'D 
‘O704 ‘ISI -0766 "1814 3804 
"0594 *1202 "0516 *1729 *3757 
‘0537 “1112 "0432 *1635 *3511 








"0000 ars a ‘0000 ies a 
‘OlIg : ‘ ‘0286 "0742 "1384 
‘O150 ‘ ' “O190 ‘oO: "0530 
‘O1d4 ° . ‘0218 


*1137 
"0553 
"0305 "6456 


‘0156 "0003 
"0040 “OO12 
"0076 


‘oo1dS 
‘0000 ‘000 "0257 ‘0024 
‘OOO! ‘0002 ‘O282 "0034 
"0004 "0094 ‘ ‘0004 
0020 ‘o106 ' "0024 
"0026 ‘O1I9 : *0033 


"0004 "0004 . "0004 "349 *46072 
‘OOI!I ‘C099 : ‘OO13 4740 
‘OO15 ‘O104 ‘ ‘OO1S *3325 °4774 


6h "0004 ‘0094 "0026 "0004 "a "4072 


From the 3p states the only possible transition is to 3s and therefore P,, ,,=1 
The expression for the effective recombination coefficient becomes 





re o ce OD AAC. gysz’ 
a(3pSL’, 3s) = 1-677 x 10-Mt-12 es x( os mSeSt (4.1) 
k VE 
with t=10-'7. Putting k=2pn/SL we have to evaluate the sums 


¥ PulA)Cn, spsu: (4-2) 
a 


r= n 
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TaBLe I1I—(continued) 
Triplets 


Case A Case B 


A A 


3p *P 3P°S = 3p*P— 3p *D 
1437 ; ‘1136 "3512 5353 
"1093 ° "1037 "2904 +4802 
"0939 . “0907 2856 *4484 





— : “0000 ie oa 
‘O1S4 ’ "3542 *2356 "0792 
‘0210 "0253 *3407 ‘2829 ‘0846 
0234 0206 3345 *2970 "0950 


“0000 ae a “0000 — 

‘0399 "0577 "1300 "4418 "3212 
"0225 ‘O3I4 "1255 5352 *2408 
‘0157 ‘0200 °1277 *5597 ‘2211 


_ "0000 ioe — “0000 
"0092 *4177 ‘O113 ‘0692 "4045 
‘0161 "4451 ‘or61 ‘IOI "5041 
‘O19I "4505 “O1S85 "1126 "5231 


"0059 "0007 6065 *3524 "0407 
‘0136 ‘0008 *4663 4381 "0445 
‘o156 ‘0021 "4454 *4326 ‘0769 
‘O221 ‘0031 *4250 "4209 ‘1010 


"0039 ‘0023 =P 6241 *3759 
‘O152 0038 "0067 7672 ‘1978 
‘0210 ‘0060 "0268 °7074 ‘2340 
"0243 "0085 0374 6646 "2562 


=e .0000 oa oe “0000 
"0013 "9194 ‘0016 ‘0096 ‘9260 
"0020 ‘8801 "0074 "0258 “8946 
‘0027 8364 "O105 "0374 8562 


"0024 4823 °1147 ‘2721 ‘6129 
"0041 "4799 "1043 "2882 "5924 
"OO51 +4776 ‘1008 ‘2920 *5867 


“002 *4823 *1147 ‘2721 6129 
“4812 "1103 ‘2793 6038 
*4799 ‘1078 2822 "5989 


6h "0024 *4823 "1147 ‘2721 “6129 


for each spectral series of given /SL. The contributions from n <6 are obtained 
from Tables III and IV. 


The elements of the cascade matrix tend to finite limits as n->0o and from 
Table III it is seen that the elements vary slowly with n for all series in which 


6* 
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TABLE IV 


Values of DA) for Or? 


Singlets, 7 





Triplets, T= 





Triplets, 7 





P 


0°929 t 4°3 ~ : °28 
0°628 "48 2-0: 31 se 
0°54 ' 2°05 I 2°6: 76 0°49 
07435 ] “4 , 2°46 2°55 o-81 


the elements are fairly large. Also ‘Table IV shows that the D,(A) do not vary 
rapidly with m and the variation of the summand in (4.2) is largely determined 
by the factor 1/v,*. ‘The sum for n>7 may therefore be estimated on fitting 
D,C,, ; to some slowly varying function, 


D.C, :=f(vn) (4-3) 


and replacing the sum by an integral, 


‘ 
2 


- y.* py." 


F PaCnat ~ gf), (PLE) ay, (4-4) 


_ 
n= U 


2 


We consider the functional form to be chosen for f(v,). Using the results of 
Burgess and Seaton (8) and asymptotic expansions obtained by Burgess (16) 
it may be shown that, for large n, 


f(v,,)~ const. x v, 7? (4.5) 
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where p lies between 7/3 and about 0-7. The procedure employed in practice 
was to use either 


fy, = Se + : (4.6) 


D 


n 


f(v,)= 


Cc 
L 
n 


according to which gave the best fit for n= 3, 4, 5 and 6 and to determine the 
constants by fitting at m=5 and 6. It was found that the final results were not 
sensitive to the form chosen for f(v,). 

The above methods were used for all states n/ with / < 3 and for 5g, 6g and 6h. 
The contributions from n> 6, /> 3 were made using hydrogenic recombination 
data and making use of the fact that cascades from these levels to 3p must pass 
through one of the ng levels. From ‘Table II] we see that Chg g,.,, varies 
slowly with n and may be replaced by some mean value Cy spss Lhe 
contribution to the sum in (4.1) is therefore 


Se l 


N\ 7 wsrDuszl aAlS L,3pSL’ 





v=Tl=4 L=hlil n 
2 2 2! (21+1)D,, 
- 3(28 7 1)C 5 apau z > 2 = (4-8) 
n=7 l=4 nm 
since, with we, = (2S+1)(2L +1), the sum over L gives a factor 3(2S + 1)(2/+ 1). 
For hydrogenic ions 


a. 
2/+1)D,,(A)= — — S,,(A) (4.9) 


(16, 17) and therefore 
~ ~! (2/ + 1)D,)(A ) ie f : = (A) - . \ (21+ 1)D, AA) 


n=7 [=4 n° 3 ne? n” n=7l1 


: (4.10) 
‘ n> 4 
Values of }S,(A)/n® were obtained using the tables of Seaton (17) and values 
of D,,(A) for /=o0 to 3 were interpolated using equation (3.22) together with 
results previously obtained by Burgess (15). Values of 


> (2/+1)D,, (4.11) 
Cc 


' ‘ : A 
calculated for n=7 to 12 were fitted to expressions of the type —+ 5 +5 
: nn 


which were used to complete the summation from m= 13 to oo. Results of these 
calculations, given in ‘Table V, should be correct to within 5 per cent. 


TABLE V 


A= 142'10 
(T=1 x 10‘, z=}3) 
n—1 (o]/+ \ 
<i (al+ DD,A 


2 
he 
n=7 l=4 nm 
x S (2l+1)D, A 
a 1+1)D,,(A) 
n=7 l=( n 
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Table VI gives the effective recombination coefficients for the 3p—>3s multiplets. 
For 4 of the 6 multiplets the intensities for Case B exceed those for Case A by 
an order of magnitude, the reason being that in Case B a large contribution to 
the 3p populations comes from cascade through the nd states but in Case A 
the nd states have large probabilities of decay to 2p”. Of the two exceptional 
multiplets, 3p 'S->3s!P is weak for both Cases and 3p *?D->3s °P is strong for both 
Cases. ‘The reason for the weakness of 3p 'S->3s'P is that, of the 2pnd terms, 
only 2pnd'P gives direct transitions to 2p3p'S but in both Cases we assume the 
nebula to be thin for 2pnd'!P-+>z2p?'S, which effectively depopulates 2pnd'P. 
The reason for the great strength of 3p *D->3s°P is that 2pnd3F cannot undergo 
direct transitions to 2p?°P and is therefore very effective, even in Case A, for 
populating 2p3p °D. 

TaBLe VI 
O 111 effective recombination coefficients 
Singlets 
Values of 10'*agp 11, 3¢ sp (O**) 
Case A Case B 
T=10'°K 
1°07 
3 85 
4°98 
Triplets 
Values of 1o'a3, a2, 3 
Case A Case B 
L T= T=1 x 108 T=2x10'°K 
S 2°91 27°3 16°4 


r g*I0 79°3 48-1 
D 124°3 153°9 86°3 





In Section 7.5 we shall estimate the recombination contribution to the 
[O 111] forbidden line intensities. Case B gives an upper limit to «(2p? !S->2p?!D) 
since it is assumed that each recombination on an excited singlet configuration 
eventually leads to population of 2p?!S. We obtain 


A(!S+'D) 
A('S+'D)+ A(S>'P) 
x fa(2p?4S)+ ¥ SS a(2pal Ly}. 

3 UL 


n 


a(2p?4S->2p?!D) = 





For T=10' °K this gives «(2p?!S—2p?!D) = 87-0 x 10-4 cm3sec—1. 

5. Otv and Ov recombination lines. —Calculations of effective recombination 
coefficients have been made for the lines O Iv 1s?2s?3d?D,..>15725"3p *P 3. and 
Ov 187253,p 1 P—>1572535 1S. 


For Olv we have 
ie % ee 
3d, 3p I 3d,3p > o,f nl3d 
nl 
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where P37 3,=1 for Case B and P3, 3, =0-002 for Case A. The sum over nl 
in (5.1) will be much less sensitive to whether Case A or Case B is assumed. 
This sum was evaluated for Case B using hydrogenic data, an approximation 
which will introduce little error since only the ns and np series are appreciably 
non-hydrogenic and these contribute only 1-4 per cent to the total summation. 
The calculations were made for 7'= (16/9) x 10K, the reason for this choice 
being that the hydrogenic “,,; used in the O11 calculations could again be used 
(7 =10' and z= 3 giving the same value of A as 7T'= (16/9) x 10 and z=4). Also 
T= (16/9) x 104 happens to be close to the correct electron temperature for 
NGC 7027. From these calculations we obtained a3, 3,8) = 2-1 x 10712 cm® sec} 
and a 34.3,“~4x 107 cm'sec™!. These figures must be multiplied by 0-6 to 
obtain the coefficients for the line 7D;,.>?P3,. 

Since the Ov 3!P->3!S line is very weak the measured intensity will be 
uncertain and an elaborate calculation would not be justified. For small (7/z*) 
a large contribution to the effective recombination rates comes from highly 
excited states, which are closely hydrogenic. Let us define 


‘n, 3p — 2 4, jC nl, ap Ants (5.2) 
Then for Case B, 
%3 p,34s> } D3 = a (5-3) 


n=%3 
where «,, is calculated for a hydrogenic ion. Using values of C,,; 5, calculated 


for Ov singlets and values of «,,, calculated from hydrogenic data with (T/z")=0 
(16) we obtain: 


n 3 4 5 6 
a 0*40 0°13 018 O19 


i?) 
We took C,, 3, =0°20, independent of n, and evaluated Sx, for T=1-7 x 108 °K 


3 
and z=5 (17) to obtain a; p.3 .=4x107%cm'sec~, a figure which should be 
correct to within a factor of two. 


6. Data for NGC 7027 


».1. Relative intensities.—Intensities are denoted by J(A) relative to the H8 
intensity /(4861)=100. ‘The relative quantum emission rates are defined by 
g(A) = (A/4861)1(A). + Weuse J, for observed intensities, /. for intensities corrected 
for reddening and J, for calculated recombination intensities. 

The intensities of Aller, Bowen and Minkowski (§) have been discussed in 
a previous paper (7). It was found that a reddening constant, c= 1-22, could 
be chosen so as to obtain agreement between J, and J, for the stronger lines 
of H1, Het and Hert but that there was then poor agreement for the weak lines 
and in the infra-red. Suggested explanations in terms of physical processes in 
the nebula appear to be unsatisfactory and the possibility has been considered 
that the intensity anomalies are due to systematic errors in the measurements. 
In the present paper we shall assume that this is the case. We take the true 
relative emission rates to be given by q(A)=&q,(A) where €=€,(/,)&(A) and 
where the correction factors £, and &, are given in (7). For all the lines 
considered in the present paper, values of J,, J, g,, € and q are given in Table VII. 
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6.2. The electron temperature and electron density —The following ratios 
depend on the electron temperature 7’, and electron density V, (18, 19): 


g(6300) + q(6364) 





p’[O 1 | 





p’ [Or] 
p[Ou]= 


p’(O1t] = 





4(4363) 
p’(Nu]= q(6548) + q(6584) 


q(5755) 








No. 1, 1960 oxygen in the planetary nebula NGC 7027 89 


Fig. 1 gives the 7T,, N, relations obtained using the observed ratios and atomic 
data from (19), (20) and (21). The ratio p’[O1] may be uncertain since /(§577) 
is weak and could include a contribution from airglow radiation but it is 
satisfactory that p”"[O1] in fact gives good agreement with the other curves. 














Fic. 1.—7,, N, relations for NGC 7027 from the forbidden line ratios (equations (6.1)) 
and from the surface brightness, SB (equation (6.2)). We adopt Te=1'7xX10'°K and 
Ng=2:°2 X 10° cm-3 


In Fig. 1 we include the curve calculated from the relation 
log N,= 5-97 + } log S,— } log (Ar) + } log (t®?/bye7,) + 3c (6.2) 
where S, is the observed H§ surface brightness in ergscm~*sec~!, 4 the nebular 
radius in seconds of arc, r the distance in parsecs and t=10~47',. We use values 
of S,, A and r from (22), 6,e’, from (16) and ¢ from (7). 

The different 7,, N, relations are seen to give reasonably concordant results 
but p[O11] gives a density lower than that obtained by the other methods; this 
difference may be due to local condensations (19). From the results of Fig. 1 
we adopt 7,= 1-7 x 10! K and N,=2:2 x 1ofcm~3. 

6.3. Abundances from forbidden lines.—The effective recombination coefficient 
for a line nn’ in the spectrum of an ion X*” we denote by «,,,,(X*™"). The 
absolute emission rate for the H§ line is 

9(4861) = N.N(H*)aq, .(H°) (6.3) 
where, using Case B for H1, 2, .(H®) in cm sec~! is 2-99 x 10°" for T,=1 x 104 °K 
and 1°595 x 10-'4 for 2x 10 K (16). For the t+; forbidden line of some ion B 
we have 

9,,=N(B)Aj- (6.4) 
The total ion density \(B) is expressed in terms of the excited state density V,(B) 
using relations of the type (19) 


N(B)P(i) 
S 


N( B)= (6.5) 


where 


(6.6) 


the sum being over all ground configuration states. 
The quantum emission rate relative to q(4861)= 100 is 
qij _ NAB) Ais 


100 N,N(H*)ay 9(H") ’ 
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We are interested in the total numbers of ions in the ionized region and may 
therefore equate N(H~) to the total hydrogen density N(H). From (6.5) and 
(6.7) we obtain 
N(B)_ qj SN, 
N(H) 100 Aj;P(i) 

Abundance ratios O°/H, O*/H and O*? H are given in Table X. For [O1] 
and [Oi], P(c)/S is approximately proportional to N, and the abundance 
ratios insensitive to .V, but, because the [O11] lines are more highly forbidden, 
O*/H is more sensitive to N.. 

7. Interpretation of the O11 line intensities 

The O11 lines are excited both by fluorescence and by recombination. We 
shall show how the contributions from these two processes may be separated. 
From the recombination contribution we obtain the O*? abundance and from the 
fluorescent contribution we obtain the overall efficiency of the fluorescent 
mechanism. We also consider the contribution to forbidden line intensities 
due to recombination. But first of all we calculate the optical depths for Out 
resonance lines. 

7.1. Optical depths.—The optical depths for absorptions from the 2p?!D 
and 'S metastable states are deduced from observed absolute intensities. ‘The 
optical depth for 7k absorption is 


a4, o(H). (6.8) 


*R, 
Tik= N,«;;,dR (7.1) 

J Ry 
where A, is the inner radius of the nebula, R, the outer radius and «,;, the 
absorption coefficient. Suppose that 7 emission is observed and that this is 
forbidden-line radiation for which the nebula is optically thin. Let Fj; be the 
flux of ¢-»j quanta at the surface of the nebula (F;;= S;;/hy; 


; where S,; is the 
surface brightness). “Then 


- fN,4,;dV 
47R,? 


F 


ij 


where the integral is over the nebular volume. From this we obtain 


g 
N,dR=GF,,/A (7-3) 


ars 
7 Ry 


where @Y is a geometrical factor. For N; uniform in R,<R<R, we have 
G=3R2/(R2+ R,Ryo+ R,?). From (7.1) and (7.3), 


/ 
Tin = Gei_Fiy| Ay. 7.4) 


For thermal Doppler broadening the absorption coefficient at the line centre 


7e M a 
=——| —; f(R,t) (7-5) 
- 7yN 4 ’ 4° 
MV; p. ( 2k ) 
where 7p, is the frequency and ./ the ion mass and where, in the notation of 
Section 3.2, 


f (Rk, i) = he, —€ CR) z(t) P. (7.6) 


Table VIII gives optical depth for lines ending on 2p?!D calculated using the 
measured HH surface brightness (22), multiplied by 10° to correct for extinction, 
and relative emission rates of Table VII. It is seen that the nebula is optically 
thick for some, but not all, of these lines. 
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Optical depths 7(2p* ‘D-+-1pnl'L) for O 111 lines in NGC 7027 
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10°4 . 20° 70°6 
3°93 0°24 oR 21°4 


Using relative emission rates of Table VII and calculated transition 
probabilities (21) we obtain 


N(2p2!D) _ 
Nap? ts) 4600 (7-7) 


and using equations of the type (6.5) we obtain 


Nia *P) Bs 
N(@p™'D) = 713: (7.8) 


The nebula is therefore optically thin for transitions ending on 2p715 and 
optically thick for those ending on 2p?3P. 

The optical depth for the absorption line 2p?°P,->2p3d*P, is of importance 
for the fluorescent mechanism. Expressing the population in the form (6.5) 
we obtain, using collision strengths from (23) and transition probabilities 
calculated from (24), 


P@P,)= x? +0°46x + 0°08 } 
PO@P,)=3x2+ 15x (7.9), 
P(2P,) = 5x2 +0-18x | 
with x=10-?N,/7T,1?2. This gives 1-0: 3°0:3°9 for the ratios of the 2p?3Py , 2 
populations. For the optical depth we obtain finally 
7(2p?3P,, 2p3d 3P,) = 2:2 x 104. (7.10) 
7.2. Separation of fluorescent and recombination emission rates.—We denote 


fluorescent emission rates by g;, recombination rates by g, and corrected observed 
rates by g=q,;+ 4, From Table VII we have 


q(3p 3S>353P) Lo (7 11) 
qap*D>35°P) > 


The corresponding ratio for recombination is 0-18 (Case B) or 0°023 (Case A) 
and the ratio for fluorescence is 





9(3p?S->35 *P) - C3q°P) 3p98 
g(3p °D 38 3P) C3a°P, 30°D 





=14°9 (7°12) 


it being noted that }A(SLJ--S’L’J’) does not depend on J. The fluorescent 
> 


ratio (7.12) is equal to the ratio of the 3d>3p transition probabilities and should 
be determined accurately. Since the observed ratio is intermediate between 
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the recombination and fluorescent ratios the relative contributions of the two 
processes may be determined. For Case B we obtain: 
q(3p°S-+35°P) (3p *D->35*P) 
4°4 
3°4 


Fluorescent 65°3 
Recombination 06 


Observed 5°9 78 


Since q,(3p*5—>3s%P) is always small, the deduced value of ¢,(3p*D—>3s*P) is 
practically independent of whether Case A or Case B is assumed. 

Having separated the fluorescent and recombination contributions by 
considering the relative intensities of the two multiplets we may check that theory 
agrees with observation for the fine structure intensities within the multiplets. 
For recombination we use the fact that N(SLJ/) is proportional to (2/ +1) 
(Section 3.1). It is here assumed that the states O*%2p?P,.,?P5. are populated 
in the 1:2 statistical weight ratio. ‘lhe density of NGC 7027 should be high 
enough to ensure that this is a good approximation. ‘The position when this 
condition is not satistied has been discussed in (25). For the fluorescent 
mechanism we solve the problem of cascade from 3d°P, taking J quantum 
numbers into account. 

‘The comparison of fine structure emission rates is given in Table IX. ‘The 
agreement between theory and observation is seen to be very good. For 
3p °S>3s°P the upper state has no fine structure and the agreement therefore 
provides a check only on the intensity measurements and on the basic assumption 
of LS coupling. ‘The 3p*D-—-3s°P multiplet is more interesting. Here the 
fluorescent and recombination mechanisms give quite different fine structure 
intensities. It is seen that the combination of g, and g, which gives the observed 
intensity ratio for the two multiplets also gives the observed relative intensities 
for the components of the ®>D--*P multiplet. 


TABLE IX 


Quantum emission rates for O11 2p3p—-2p3s in NGC 7027 


J 


0%3 
o'8 
ome) 
03 
o"4 
3°4 
o'6 (B) 
0°06 (A) 
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In ‘lable 1X we include the emission rates for the singlet line, 3p !P—+>3s'P, 
to which the fluorescent mechanism does not contribute. It is satisfactory that 
the observed emission rate is intermediate between those calculated for Cases A 
and B. 

7.3. The O** abundance.—The abundance ratio is determined from 

N(O**) _ 9(3p*D->3s9P) a4, 2(H®) (7.13) 

N(H) 100 %3y »p,39¢p (OF?) * _— 
Using ¢,(3p *D, 3s°P) = 3-4, values of «, .(H°) quoted in Section 6.3 and values 
of %3),sp, 3sep for Case B from Table VI we obtain N(O**)/N(H)=0-64 x 107%. 
This result is practically independent of the assumed value of 7, (we obtain 
0-65 x 10-8 for T,=1 x 10'°K and 0-63 x 10-* for T,=2 x 104°K) and is not 
sensitive to whether Case A or Case B is assumed for Ou triplets (the O** 
abundance would be increased by a factor of 1-24 if Case A were assumed). 

7.4. The efficiency of the fluorescent mechanism.—TVhe probability that absorp- 
tion of a single Hert Lyx photon by O*? will give fluorescent lines is 





CW) = oo167 
2 Csdy,sy. = 0167. 
This is the quantity P,, used in Section 2. 
The total emission rate for fluorescent quanta is 
IO I11)= = 9 (3d 3P, >3p i 
L 
The overall efficiency of the fluorescent mechanism is 
—  Q¢Omt) 
~ 9(Hermt Lyx) ° 
From the cascade theory of Section 4 we obtain 
) — thyeb 3m -3p 
9 (O UL) = 1-659 (3p FS--35°P) 
and for the He11 Lyx emission rate we have 


ten, 19( He? 


) 9(Hem, 4> 
a afte) O(He tl, 4-+3). .18) 


Case B should certainly apply for He*. We then have 


2(Heu Lyx) = 


22», 14(He*) = (1 —X) 3 x, (Her) (7.19) 


where X, which is approximately equal to 0-3, is the fraction of excited state 
recombinations which populate 2s (16). Using data from (7) and (17) we 
obtain 

Lo ( 

2p,1s 

Sr i eee , .20 

a4, 3(He*) (7.20) 
practically independent of T,. The fluorescent efficiency is therefore 


y a. qi(3p°5 35°P) 
R=0° “€ 7: 
O°19 X 1°05 X q(He tl, 4 >3) ( ”" 


and substitution of ¢/(3p °S—3s°P) = 65-3 and q(Hell, 4--3)= 48 gives #=0-43. 
7.5. The recombination contribution to forbidden line intensities.—It is usually 
assumed that collisional processes alone are responsible for forbidden line 


excitation. We are now in a position to check that this is justified. We have 





g(2p2'S->2p2'D) = 2H '820"'P x g (3p 3D-+35 3P). 22) 
X3n*D,389P 
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In Section 4 we obtained an upper limit of 87 x 10-“ cm*sec™ for a 913 9,2 1p: 
For the same temperature, 7',= 10* “K, we have 5, sp ..sp = 154 x 10-4 cm?sec™}, 
With g,(3p *D->35°P) = 3°4 we therefore obtain g,(2p?!5-—2p?!D)<1-9g. Com- 
parison with the observed rate, g(2p?!S—2p?!D)=27, shows recombination to 
be unimportant. 

8. Abundances of O** and O*° 

8.1. The O** abundance.—Recombination of O** ions in the 15s?2s? ground 
configuration will populate O*? levels with configuration 1s*2s*n/. ‘Vhis could 
give the observed Orv line, 37D;.>3?P 3, 43412. A further weak O1v line 
observed by Aller, Bowen and Minkowski (5), 


1572s2p(8P)3d 4F y).->18*2s2p(2P)3p4D.2, 3737, 


would not be expected in the recombination spectrum. It would be worth 
looking for the other members of the multiplet in order to check the identification. 

We use the 37D,.~37P;,. line to obtain the O* abundance. Using 
2p, , 3¢p,, = 12x 107 cm sec™! calculated for Case B and T= 1-78 x 10° K 
(Section 3.5) we obtain V(O**)/ N(H)=0°14 x 107%. ‘The O** abundance would 
be greater by a factor of order 10° if Case A were assumed. Since the O11 triplet 
resonance lines have central optical depths of order 10 (Section 7.1), and since 
the number of O** ions exceeds the number of O** ions by a factor of about 4, 
the optical depths will certainly be very large for the O Iv resonance lines and we 
may expect Case B to be approached but the discussion of fluorescent efficiencies 
(Section 2) suggests that the O*? abundance may exceed that calculated for 
Case B by a factor of 2 or 3. 

We adopt an O** abundance twice as large as that calculated for Case B. 

8.2. The O** abundance.—A single weak O v line is observed, 3! P->3!5S, A 5112. 
Using %3:p.31g =4 * 107 cmsec~! calculated for Case B and 7,=1-7 x 10° “K 
we obtain N(O**)/ N(H)=0°04 x 107%. 

As in the case of O**, we adopt an abundance for O*° which is twice as large 
as that calculated for Case B. 

g. Summary and discussion 

g.1. Summary.—TVable X summarizes our results for the oxygen ion abun- 
dances. ‘The chemical abundance ratio for NGC 7027, N(O)/ N(H) = 1-2 x 107%, 


TABLE X 
10°NCX)/N(H) 1018ay, (X) I.P. (eV) 


0°0035 7°6 
roo} fe) 


13°60 
cee 
54°93 
77°35 
113°86 


I 
O°15 3°5 
o'64 I 
0°25 ‘9 


0:08 *32 f 138°06 


agrees closely with the chemical abundance ratios obtained for atmospheres of 
early type stars (26, 27). In ‘Table X we include estimates of the optical depths 
for photo-ionization of the oxygen ions. ‘These are calculated from 


79(X) = 16 x 10° N(X)a; (X) (9.1) 
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where 1-6 x 10!’ cm is taken to be the radius of the nebula (22) and where the 
ay, are the threshold photo-ionization cross-sections, calculated from formulae 
given in (8) and (20). 

g.2. Discussion of sources of error.—For the total oxygen abundance the main 
uncertainty is in the observed intensities. Only further observational work can 
decide whether we are justified in introducing the correction factor € which 
brings the measured H1, Het and Het! intensities into agreement with theory. 
If this factor is not introduced it is difficult to obtain any value for the O/H 
abundance ratio since we then have no satisfactory theory for the H1I spectrum. 
If we had not introduced the correction factor but had still assumed the 
recombination theory for the strong H1 lines we would have obtained a ratio 
N(O)) N(CH) = 2°6 x 1074 

‘The main uncertainties on the theory side are concerned with the question 
of how closely Case B is approached for Oi, Otv and Ov. Improved 
calculations will be possible when a satisfactory theory has been obtained for the 
eficiencies of the fluorescent mechanism. So far as the total oxygen abundance 
is concerned it is very fortunate that, for reasons explained in Section 4, the 
effective recombination coefficient for O 111 3p *D->3s ’P is not sensitive to whether 
case A or Case B is assumed. ‘The abundance of O**, which is the most abundant 
oxygen ion, should therefore be determined fairly accurately. ‘The abundances of 
O** and O*° are determined a good deal less accurately but, so long as departures 
from Case B are not much greater than we have assumed, this does not introduce 
serious uncertainty in the final chemical abundance ratio. Further improvement 
in the accuracy of the O** and O*° abundances is still very desirable since a 
quantitative study of the ionization equilibrium for oxygen should eventually 
provide important information about the stellar ultra-violet radiation field. 

9.3. The determination of oxygen abundances in other nebulae.—It would be 
of interest to determine O/H abundance ratios for a number of nebulae. ‘The 
ions O** and O*° do not make a large contribution to the total oxygen abundance 
in NGC 7027 and since this object has particularly large abundances of highly 
ionized atoms we may expect that O*' and O*° may be neglected in most 
nebulae. ‘The main need is to determine O** abundances from recombination 
lines in addition to O°, O* and O*? from forbidden lines. The O*% abundance 
is best obtained from the Out multiplet 3p °D->3s°P. ‘The wave-length interval 
3750A to 3800A contains all five lines of the Or multiplet and, in addition, 
three H1 lines and two He lines. In NGC 7027 these ten lines all have 
comparable intensities and this may also be expected in other high excitation 
nebulae. Further careful spectrophotometry of this 50A interval in selected 
nebulae would be of considerable interest. From the relative intensities of the 
Ou lines the fluorescent and recombination contributions may be separated. 
The ratio of Our recombination intensities to H1 intensities then gives the 
O-*3,H abundance ratio and the ratio of Or fluorescent intensities to Hell 
intensities give the efficiency of the fluorescent mechanism. 
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VARIABLE STARS IN NGC 2257 
JF. B. Alexander* 


(Communicated by the Radcliffe Observer) 
(Received 1959 November 30) 

Summary 
From the blink-comparison of plates taken at the Radcliffe Observatory 
over a period of about six years, 28 variables have been discovered in the 
globular cluster NGC 2257 which is on the outskirts of the Large Magellanic 
Cloud. Visual estimates of the brightness of these variables have been made 
on a step basis relative to a sequence. For six variables with a large amplitude 
of light-variation, periods ranging from 0°51 to 0-69 days have been derived, 
and the light-curves obtained are like those of Bailey a and b type RR Lyrae 
variables in the galaxy. It has not been possible from the available material 
to work out the periods of the many variables with a much smaller range of 
light-variation, but at least some of these stars have the gradual rise to 

maximum found in Bailey ¢ type variables. 





1. Introduction.—The globular cluster NGC 2257 (R.A. 6®29™9, Dec. 
~64 27’; 1950) was first regarded as a possible outlying member of the Large 
Magellanic Cloud by Shapley and Mohr (1), who assigned to it an integrated 
photographic magnitude of 12™-4. The degree of central concentration is not 
high (see Plate 1), and on good plates taken with the Radclitfe 74-inch reflector 
it is possible to resolve some stars at the centre of the cluster at a magnitude 


corresponding to the minima of the RR Lyrae variables. No determination of 
the integrated colour of the cluster has yet been made. 

Variable stars were first found in the cluster by Dr A. D. Thackeray in 1953. 
A series of plates extending from September 1953 to May 1959 have been taken 
with a view to the further discovery of variables, and, if possible, to the eventual 
determination of periods and light-curves. There exist 66 plates, mostly of 
30min exposure, taken on Kodak 103a-O (without filter) at the Newtonian 
focus of the Radcliffe 74-inch reflector. Of these, 37 were taken by Dr A. D. 
Thackeray, 25 by Dr A. J. Wesselink, 1 by Dr T. D. Kinman and the remaining 
3 bythe writer. ‘There is also one 30 min exposure in the yellow which is sufficient 
to show that the brightest stars in the cluster are not blue. 

2. The blink-comparisons.—Twelve pairs of plates were examined by the 
writer in the blink comparator which is on loan to the Radclitfe Observatory 
trom the Royal Observatory, Cape. On the basis of this survey, 16 stars were 
considered with certainty to be variables, and a further 11 stars were regarded as 
likely variables. Of this latter group, 7 were afterwards considered to be certain 
variables when the various estimates of brightness described in Section 3 had been 
made. 

Some idea of the degree of completeness of the examination can be obtained 
as follows. The 23 certain variables are divided into two groups according to 
whether the amplitude of the observed light-variations is less than or more than 
eight steps on the scale described in Section 3. The blinking was done in such 
a way that a method described by Van Gent (2) can be applied, and the results 
of this analysis are given in Table I. Group 1 is the large amplitude group. 

* On secondment from the Royal Greenwich Observatory, Herstmonceux. 
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TaBLe | 
Predicted number Number found 
in cluster* in survey 
Group 1 8-0 8 
Group 2 16-6 15 

The probability of the discovery of a variable from the blink-comparison of 
one pair of plates (i.e. the quantity « in van Gent’s notation) is 0-53 for Group 1 
and o-18 for Group 2. 

Although the cluster is not very heavily concentrated towards its centre, it is 
likely that several variables have been missed in the blinking because of crowding, 
but no allowance is made for this effect in any simple statistical method such as 
that used above. ‘This is to some extent confirmed by what follows below. 

The blinking done by the writer was carried out independently of any other 
investigation. Dr A. R. Sandage has also searched the plates for variables, and 
an unpublished list of 19 certain variables found by him contains 5 not discovered 
by the writer. ‘There are therefore 28 stars in all regarded as certain variables in 
at least one of the two surveys. ‘These are indicated by numbers in Plate 1. 
There are four C.P.D. stars in the field of the plates, and these have been used to 
obtain the approximate distances of the variables east and north (AX and AY 
respectively) of an origin near the centre of the cluster. These are given in 
Table I] below. The letters refer to the comparison stars mentioned in the next 
section. 

Tasie II 
Positions of variables in NGC 2257 for 1950 (seconds of arc) 
Variable AX ay Variable AX AY 
+27 +45 15 —14 
—4I —28 16 ~ 7 
+28 —24 17 91 
+99 +28 18 
+51 se 
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Comparison AY 
A 5 +46 
B +70 70 
Cc +47 47 
D +60 63 
E + +42 72 
A,G,. 
an 
variable on one pair of plates. 1 is the total number of pairs which are blinked. If we denote 
by aj the number of variables which are found a total of i times during the course of the m blink 
comparisons, then A, and G, are defined by the relations: 


* Evaluated from Van Gent’s equation N= a is the probability of the discovery of a 
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NGC 2257 in blue light with variable and comparison stars indicated. 
Scale of reproduction: 1":g/mm. 


J. B. Alexander, Variable stars in NGC 2257 
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3. The visual estimates.—Estimates of the brightness of the variable stars 
relative to a sequence were made on all the plates by means of visual inspection 
with an eyepiece. The five sequence stars A, B, C, D and E are shown in Plate 1 
and their positions are given in Table Il. The gap between successive sequence 
stars was taken to be either four or six steps, the interval being kept fixed from plate 
to plate. A star equal in brightness to the brightest sequence star A is assigned 
the number o, and the corresponding values for the other sequence stars are 
given below. 

Sequence Star A B C D E 


Step ° 6 10 16 20 
Estimates have been made at least once on each plate of all the stars regarded as 
certain or likely variables in the blink surveys. In many cases, owing to the 
quality of the seeing, all that can be said is that a variable is fainter than a certain 
sequence star. 

No photoelectric calibration of this sequence exists at present. 

4. Individual results.—Owing to the spacing of good plates, it is very difficult 
to find the periods of the variables from the estimates. One serious problem is 
the difficulty of obtaining good conditions for observing on consecutive nights 
during the summer months at Pretoria. Moreover, plates have to be taken 
when the cluster is at rather a low altitude if a reasonably long series of observations 
is to be obtained in a night. 

Periods have been derived for six variables with large amplitudes of light- 
variation and steep rises giving well-defined maxima. ‘The elements for these 
variables are given below: 

Tas_e III 
Elements of six variables in NGC 2257 
Maximum at ie a 
Variable JD 2436000 + days~! days 
280°29 "$5047 0°54040 
254°43 97767 0°50565 
‘ 280°34 82884 0°54679 
) 281°31 45843 0°68567 
7 254°38 ‘77141 0°56452 
I4 695°26 ‘97287 050688 
The corresponding light-curves for these variables are plotted in Fig. 1. At the 
time when the estimates described in Section 3 were made, in many cases it was 
felt that there was some doubt about the estimate due either to the general quality 
of the images on the plate or to the effects of crowding. Such step-estimates are 
plotted as open circles in the figure. 

It should be borne in mind that these light-curves will differ from the true 
light-curves due to the effect of the finite exposure time. In particular, the 
steeply rising portions of the Bailey a type light-curves will be in reality some- 
what steeper than they appear in Fig. 1. 

Various details for the 28 variables are given below : 

Vi. P=0'54040 days. A Bailey a type light-curve. 

V2. P=0-50565 days. A Bailey a type light-curve. 

V3. The step-estimates vary from 12 to 18. 

V4. <A nearby bright star makes estimation very difficult. 
Amplitude of light-variation appears to be large. 
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P=0-54679 days. ‘There is a fairly bright star nearby, and there is also a 
very close faint star which is only resolved in conditions of excellent seeing. 
The estimates are therefore very difficult, and refer to the total light of the 
variable and its faint neighbour. This should be borne in mind especially 
when the fainter portions of the light-curve are compared with the corre- 
sponding parts of those of the other variables. 

P=0-68567 days. This is the longest period of the six which have been 
derived. ‘The amplitude of the light-variation and the asymmetry of the 
light-curve are both markedly less than for any of the other five variables. 
P=0°56452 days. Nearby star makes estimation rather difficult. A 
Bailey a type light-curve. 
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Fic. 1.—Light-curves of variable stars in NGC 2257. Abscissae: phase with o-2P overlap 
right and left. Ordinates: step-estimates relative to sequence. 


(a) V1; P-'=1-85047. & eo: FF 
(6) V2; P-'=1-97767. (f) V14; P 
(c) V5; P-*=1°82884. (g) Vis; P 
(d) V6; P-'=1-45843. (h) Vis; P 


Lox 5-77141. 
1— 1-97287. 
1— 4:00418. 
1— 3:00418. 
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Step-estimates from 12 to 18. 
Step-estimates from 13 to 18. 

. Step-estimates from 13 to 18. It would appear to have a gradual rise to 
maximum. 

. A close bright star makes estimation difficult. Step-estimates vary from 
10 to 19. ‘The rise to maximum appears to be fairly steep. 

. Step-estimates from 12 to 19. 

Step-estimates from 10 to 17. 

. P=0-50688 days. A Bailey a type light-curve. 

. The observations are satisfied tolerably well both by a period of 0°249739 
days and by one of 0:332870 days (see Fig. 1). It is by no means certain 
that the true period is either of the two possibilities above. 

. Step-estimates from 12 to 18. 

Step-estimates from 13 to 18. 
Step-estimates from 13 to 19g. 
Step-estimates from 14 to 18. 
Step-estimates from 13 to 18. 
Step-estimates from 12 to 17. 

. Step-estimates from 14 to 18. 

. Nearby bright star makes estimation extremely difficult. 

. Estimation difficult. Step-estimates from 14 to 18. 

Step-estimates from (g) to 16. 

. Step-estimates from (6)to18. The light-curve has a steep rise to maximum. 
Two nearby stars make estimation extremely difficult. Amplitude of light- 
variation appears to be large. 

. Very close to a fairly bright star. Certainly a variable of large amplitude. 

5. Discussion.—A recent survey of the red globular clusters in the Large 
Magellanic Cloud has been made by Hodge (3). Although NGC 2257 is just 
outside the region he considers, he remarks that the clusters probably extend 
beyond the limits of his search on the eastward side of the Cloud. 

As no photoelectric calibration exists at present, it is not yet possible to give 
the magnitudes of any stars in NGC 2257. A very rough idea of the brightness 
of the variables can be obtained by comparing the plates of the cluster with those 
of NGC 1466 and NGC 1978 (both in the Large Magellanic Cloud) also taken 
with the Radcliffe 74-inch reflector. Since plates of these two clusters were 
taken in similar conditions to and in some cases on the same nights as NGC 2257, 
a comparison can be made. It seems that the mean magnitude of the RR Lyrae 
variables in NGC 2257 is the same as those in the other two clusters to within a 
third of a magnitude. (The median photographic magnitude of the RR Lyrae 
variables in NGC 1466 and NGC 1978 is 19™-2 (4)). 

Hence from the above there is no contradiction with the view that NGC 2257 
is an outlying member of the Large Magellanic Cloud. 

The position of the 28 variables is such that most of them must be regarded 
as members of the globular cluster NGC 2257. It is possible that one or two 
may be field stars in the Large Magellanic Cloud, but none of the variables is far 
enough away from the cluster centre to be in a position where cluster members 
are very scarce. NGC 2257 must therefore be regarded as being rich in variables. 
Of 72 globular clusters in the galaxy which had been searched for variables prior 
to 1955, in only 13 of them had more than 28 variables been found (5). 
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The six variables for which periods have been derived have light-curves like 
those of Bailey a and 6 type RR Lyrae variables in the galaxy. For the five 
Bailey a type variables with periods between 0-50 and 0-57 days the step-estimates 
of the maxima and minima differ little from star to star, and the amplitudes of 
light-variation are all about 14 or 15 steps (an exception being made for V5 where 
the minimum cannot be determined accurately). ‘The five stars have very steep 
rises to maximum, but, since the apparent slope of the rising portion of the light- 
curve is very sensitive to small changes in the assumed period, no quantitative 
assessments of the asymmetry have been made. 

V6, with a period of 0-69 days, has a smaller amplitude of light variation than 
the other five stars for which periods have been determined. Moreover, the 
much more gradul rise to maximum is typical of the Bailey 5 type variables found 
in the galaxy. 

It is reasonable to assume that the other 22 variables found are all of the 
RR Lyrae type. Although periods have not been deduced for these, all show 
clear indications of a change in brightness during times of an hour or so. The 
median magnitudes of the smaller amplitude members are all about 15 or 16 on 
the step-scale used in the estimates. 

It is interesting to compare NGC 2257 in the Large Magellanic Cloud with 
NGC 121 in the Small Magellanic Cloud (6). The latter is a globular cluster of 
somewhat smaller angular extent but with a much greater central concentration. 
Its known variable star population consists of two long-period variables and 
three RR Lyrae variables with periods between 0-50 and 0-64 days, there being 
apparently none of the small amplitude variables which are very abundant in 
NGC 2257. ‘This difference between two individual clusters is, of course, no 
indication of any fundamental dissimilarity between the globular cluster systems 
of the two Magellanic Clouds. The variable star content of globular clusters in 
the galaxy varies enormously from cluster to cluster, and, in fact, from the available 
evidence neither NGC 2257 nor NGC 121 would appear to possess any extra- 
ordinary peculiarities compared with their counterparts in the galaxy. 
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Summary 


In connection with an Ha atlas of the Southern Milky Way being 
produced at Mount Stromlo, a catalogue has been made of the H 11 regions 
distributed about the galactic equator within galactic latitude limits +15°, 
and extending from galactic longitude 190° to 12°. The results of this 
survey have been compared with those of previous surveys overlapping 
the region investigated. 





In 1957, a survey programme of the Southern Milky Way for the detection 
of Hit regions was initiated at Mount Stromlo Observatory by Professor 
B.J. Bok. This programme was carried out with a Meinel- Pearson §-inch f/1 flat 
field Schmidt which was mounted as a counterweight to the 6-inch Farnham 
refractor used as the guiding telescope. ‘The camera was delivered to Mount 
Stromlo in December 1957, when observations were commenced, the Milky 


Way programme being completed in April 1959. 

Optical quality of photographic images depends in large part upon the quality 
of the objective filters used in the particular survey. ‘The filters used here were 
Chance OR1 with Kodak 103a-E emulsion for the Ha plates, and Chance OY1 
with Kodak 103a-D emulsion for the yellow comparison plates. ‘The glass disks, 
each 8-25 inches in diameter, were 8mm (OR1r) and 5mm (OYr) thick. After 
final adjustment of the camera it was found that the smallest well-exposed images 
on the 103a-E emulsion were 17 microns in diameter. This figure was regarded 
as acceptable in view of the limiting resolution imposed by the emulsion. 

The measured colour transmissions of the filters are given in Table I. 


TABLE I 
Filter Wavelength A Ha 
5600 5800 6000 6200 6400 6563 6600 
Chance OR1,9 mm ___ Transmission ° ° ° ° . 52 61 
Chance OY1, 5 mm (per cent) ° 48 83 88 88 


At wavelengths of 6300A and 6364A, where prominent night sky lines of 
[O1] occur, the transmission of the OR1 filter was less than 2 per cent. 

The survey was carried out in two parts: in 1958, H« and comparison plates 
were obtained on 26 centres of the region along the Lund galactic plane commenc- 
ing at /=1g0 and extending to /=12° through southern declinations, together 
with 26 North and 26 South plates centred at b= +9°. The plate centres are 
7 degrees apart in longitude and the diameter of the coma-free field is 12 degrees. 
The area of the sky surveyed, then, is 194 degrees long and a maximum of 30 
degrees wide. 
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The majority of the plates were obtained by Rodgers, the remainder by 
Campbell. The possible H11 regions were detected by a blink comparison of 
the H« and yellow plates. The plate pairs were exposed to approximately the 
same limiting stellar magnitude, which is near 14 on the visual scale. The 
exposure times averaged 20 minutes for the red plates and 7 minutes for the 
yellow. The blinking technique was preferred to the ‘‘ negative—positive super- 
position ’’ method used by Gum (1), being simpler and one that takes advantage 
of the excellent optical quality of the camera which may otherwise be lost in 
any copying process. The plates were blinked by Whiteoak, Rodgers and 
Campbell, and the positions, dimensions and estimated brightnesses of the 
H11 regions are given in Tables II and III. 

The columns of the ‘Tables are: 

1. Current catalogue number. 
New galactic longitude, /". 
New galactic latitude, 5". 
Old galactic longitude, /'. 
Old galactic latitude, 5". 
. Right ascension, 2959. 
. Declination, 8)95». 
. Dimensions (minutes of arc). 

g. Estimated brightness on a scale of bright (b), medium (m), or faint (f). 

10. The number of the object in other catalogues. 
The tables are followed by a list of remarks conce rning individual nebulae. The 
coordinates have an estimated accuracy of +3 min of arc in right ascension and 
declination. Comparison is made with surveys of Gum (1), Sharpless (2), 
NGC (3), Hase and Shajn (4), and Bok, Bester and Wade (5). NGC or IC 
numbers in parentheses refer to clusters in emission regions. Catalogues in 
column 10 are designated thus: HS: Hase and Shajn; E: Sharpless; BBW: 
Bok, Bester and Wade; and G: Gum. Reduction from equatorial coordinates 
was effected using the computer SILLIAC of the University of Sydney. 

Table I] contains regions with dimensions greater than 4 min of arc, which are 
unambiguously H11 regions. Table III contains objects with diameters less 
than 4 minutes, found to be bright on the Ha plates, some of which are stellar 
or semi-stellar and may include emission-line B stars, planetary nebulae, Wolf- 
Rayet stars or variable stars. The known variables listed in the Variable Star 
Catalogue of P. V. Kukarkin et al. (6) and detected in the survey have been 
omitted from Table III. In the higher latitude fields, only objects definitely 
found to be H11 regions are noted in Table II. 

It is our experience that in this survey we have reached fainter limiting 
emission measures in a given region of the sky than did Gum, primarily because 
of the increase in resolution of our camera. This has resulted in the detection 
of large areas of faint diffuse emission particularly in the region of Vela-Puppis. 
Gum (7) discovered large areas of Hx emission in Vela and Puppis, and Abt, 
Morgan and Strémgren (8) photographed this object at low altitude with an 
8-inch, f/1 Schmidt camera. The nebula as defined by the latter work is inside 
the area of the present survey. We have found that the outlying fragments in 
Antlia and Pyxis are in fact joined to the brighter main body of the nebulae in 
Vela and Puppis, that faint emission fills the shell structure of the brighter part 
of the nebula near 07"50™; — 45°, and that there are extensions to the nebula, 
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to the outer parts in Antlia and Pyxis, near 5 Canis Majoris, and near the galactic 
equator to longitude 249°. The overall dimensions of the nebula then become 
27° x 42°, elongated along the galactic plane and symmetrical about it. Consider- 
able difficulty lies in the adequate description and cataloguing of this nebula due 
to its size and complex brightness distribution. For this reason Table II lists 
only the extensions of the Vela-Puppis nebula found here but not described by 
Gum (9) or in the Yerkes survey of Abt, Morgan and Strémgren. 

To give pictorial representation of these nebulae, Mount Stromlo Observatory 
is preparing a complete Hz atlas of the region surveyed, together with key charts 
giving right ascension, declination and new galactic coordinates based on the 
recent recommendations of the Commission 33b of the International Astro- 
nomical Union. 


Tas_e II 
Table of regions greater than 4° diam 


b1 Zasse 8:93. Dimen- Bright- Comparison 
sions ness 
hm ° ’ (min of arc) 
07 02 —09 30 150X150 NGC 2327; HS 111, 
~>07 14 >—1224 113,114; G2,3 

07 02 —1027 14X14 IC 2177; HS 109; G1 
O7 21°5 —08 30 60x60 
07 16°71 13 09 8x8 


3 
> 


NGC 2359; HS 116; 
G4 
G5 


07 08-0 —1824 12X10 5 
(NGC 2409)?; G6 


07 29°5 -1651 17X15 
06 43 —2420 90X30 
07 29'0 —1854 18x14 
06 52 —2335 60X20 
0733°7. —1842 6x3 
07 30°! —1918 12X10 
07 17'0 —2150 6x4 
07 08 —2500 300 300 
07512 —2615 33X33 
08 13°5 —3542 48x40 
08 17°5 —3600 10X10 
0913°5 —3110 45X45 
08 36°5 —4012 100X100 
08 43°0 —4109 27X27 
08 49°5 —4154 95x80 
0849°7, —4355 30X30 
0857°5 —4333 126 
08 58°5 | —4545 13X3 
0857°5 —4716 40x40 
ego1r'o)|=— — 48 27 8x8 
09148 —4745 8x4 
09224 —5154 9X7 
10 10°0 —5609 16x16 
10060 —5715 I5X15 
1003°5 —5842 25X20 
101470 —5736 15X10 
1022 —5727 90X35 
10 24'5 —5654 12X12 
10 36°1 —5742 12X12 
10 43°5 —§818 15X15 
10 40 —5930 210X210 


G7 
(NGC 2367) 


(NGC 2467); Gog 
G10 
NGC 2579; Gi 


14 
315 

317 

318 

G 20 

Part of NGC 2736 
G 22, 23 24, 

G25 


++++4+4+44 


G 26 


NGC 3199; G28 
(NGC 3247); G29 
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BBW 25500; G 32 

NGC 3293, 3324, 3372; 
IC 2599; BBW 
25500; G 30, 31, 33 

- 2} 10 56 —58 42 210x60 (NGC 3503)? (NGC 


4ooe 
s 


21-1050 ->—61 30 3572)? G 34, a, b, 35, 
o4J)—-1109 —~-—6036 36, 37 

30 105472 —6245 8x8 

o's 1112*5 —6056 170X400 NGC 3603; G 38, a,b 
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TaBLe [I1—continued 


Oi9se 81980 + Dimen- Bright- Comparison 
sions ness 
’ (min of arc) 
260°0 ° : ‘ 7X7 
260°8 : ‘ 180 X 150 
261-4 ° rs 50X50 IC 2872; BBW 26201; 
G 39, 40 
261°8 —- 2 , - 63 3 1§ X15 G4qlI 
262°4 - . 4 E 80 x 80 IC 2944, 8; BBW 
26201; G42 
264°6 , 5 250 X 40 
267°6 
268-7 
269°4 
269°9 
272'°9 
273°9 
275'6 
276°4 
277°0 
278°7 
279'6 
281°! 
288'8 
289°9 
293°9 
294'3 
295°3 
299°6 
300°! 
300°5 
300°6 
300°6 
304°0 
304°! 


1X6 G 43 
5X15 j 
3 G 45 
12 
X13 BBW 27300; G 48a 
x 30 3 BBW 27500; G 48b 
5 BBW 27600; G 48c 
x10 BBW 27700; G 48d 
4 4 
x 50 j BBW 27901 
5 X20 BBW 28100 
X10 f BBW 28801 
‘ F BBW 28900 
BBW 29400a 
BBW 29401 
BBW 29501; G49 
BBW 29902 


i] 
NN I 


Owr Ow 


Oo 


I 
ROW NU 
° 


oOmuMnNUM 


BBW 30000; G 51 


NW NW 


CO wCOWNHOUNOUU OOM CUM 
ono 


Neen Qe han Oounuou OOo 


NGC 6164, 5; G52 

(NGC 6193); BBW 
30402; G53 

G 54 

G5 

(NGC 6231); BBW 
31100; G55 


~ 
nN 
° 


308°6 


3 


wm 


(NGC 6281); BBW 
31500; G57,a,b; 
E 1 

G 58; E2 

NGC 6337; G59; E6? 


, + 


Q 


BBW 31800a 

NGC 6334; BBW 
31800b; E7; HS 
119; G 61, 62, 63, 
64a, b,c 


aw www 
uuu + 
c 


° 
o= 000 
“SION Ou 


E4 
34.18 x 55 (NGC 6357); BBW 
32100; HS 120; E8; 
G 66 
-32 42 x8 (NGC 6383); BBW 
32301; HS 121; 
E11; G67 
~31 36 ex BBW 32300; E10; 
G 68 
—3114 x BBW 32602; E17; 
G 69 
2918 : BBW 32801; E13 
28 49 E 12 
—29 02 : F BBW 32701; E 16 
—28 46 x E 15 
—24 54 x Near NGC 6432 


nw WwW 
amu 
WN 
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TasB_Le [1—continued 


ove ENE 


X1950 


hm 
17 50°7 


17 57 
18 03 


17 59°5 


18 o4°0 
18 13°5 


18 16°0 


18 03°0 
18 09°3 
18 15°0 
18 13'2 


18 03°5 


18 22'0 
18 18°5 


18 18:0 
18 20'0 


18 27°1 
18 23°0 
18 16 


18 22°6 
1818 


18 28:2 
18 28°5 
18 28°6 
18 33°8 
18 34°0 


18 29°2 
18 44°2 
18 42°7 
18 41°38 
19 O12 
18 52°3 


31950 : 
sions 


° 


—2500 85x65 


—2315 
—2400 


90 X 35 
120 X 90 


—2254 16x16 


—2145 
—2025 


120 X 30 
100 X 35 


—1920 60X20 


—1700 
—174I1 
—17 30 
—16 36 


—1412 


— 16 36 
—1600 


—1509 
—1440 


—15 24 
—1451 
=—53.98 


—1309 15x 
—1154 180 


—o9 48 
— 09 09 
— 08 27 
—07 32 
—06 41 


—02 13 a 

03 48 7) 
—02 04 8> 
—0024 12 
+0209 20x 
+06 00 5x 


x 


4 
5 
8 
I 
I 
5 


Taste III 


of possible regions less than 4’ diam. 


bl X1950 


hm 


07 14'0 


07 28°3 
07 45°4 
07 57°2 
08 28°8 


51950 
sions 


’ (min of arc) 


—og 18 4X4 
—18 27 2x6 
—27 13 2x2 
—33 42 1X1 
—38 10 IXI 


Dimen- Bright- 
ness 


’ (min of arc) 


Dimen- Bright- 


Comparison 


f BBW 33201; HS 123; 

E19; G71 

HS 124; E20;G74a,b 

NGC 6523(M8), 6559; 
BBW 33402; HS 
126; E 23, 27, 29, 30} 
G 72,75 

NGC 6514 (M 20); HS 
125; E22; G76 

HS 129; E24; G77a 

NGC 6526; HS 134; 
E 32,33; G77b 

(IC 1283, 4); BBW 
33903; HS 140; E 39; 
G78 


HS 133; E31 

HS 136; E34 

(IC 4701); BBW 
34101; HS 137; 
E35; G79 

BBW 34201; HS 128; 
E25; G8o 


NGC 6618; (IC 4706, 
7); BBW 34202; 
HS 144, 145; E 43; 
G 81a,b 

HS 143 

BBW 34401; HS 146; 
E 44; G82 


HS 148; E 46 

NGC 6611; BBW 
34400; HS 142; 
E 42; G83 

HS 147; E 45 

(NGC 6604); BBW 
34600; E 41; G 84, 


85 
HS 150; E48 


HS 151; E49 

HS 155?; E53 

BBW 35201; HS 154; 
E 52 

HS 152 

HS 159; E57 

HS 158; E56 

HS 156; E55 


Comparison 
ness 


HS 115 
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TaBLe I11—continued 


ja Sr es0 81950 Dimen- Bright- Comparison 


sions ness 
hm ° * (min of arc) 
08 35°0 —-38 21 
08 23°6 —40 12 
08 36°9 — 39 02 
08 44°6 — 38 27 
08 19°6 — 42 33 
08 20°0 — 42 32 
o8 21°! —42 29 
08 22°3 —4255 
08 54°4 —4256 
09 OI'9 48 12 
09 25°5 = 3s $4 
09 29°6 — 56°05 
Il 03°3 —62 12 
12 16°7 —62 41 
12 32°6 —61 39 
12 33° -61 45 
12 45° — 63 34 
12 47° —6118 
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— 3140 
—29 02 
—28 50 
—21 35 
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—2°9 19 06°3 +05 32 
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Notes to Tables II and III 


Horseshoe-shaped bright region 2° in diameter, with 8’ x 5’ concentration at 07h 03.0, 
—12° 12’. Fainter region 150’ diameter, centred on 071 08M, — 09° 30’. 

Circular, connected to 1 by faint emission. 

Circular region 

Contains diffuse outer region. 

Circular region. 

Part of the Vela-Puppis Nebula ? 

Part of Vela-Puppis ? 

Circular region. 
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Notes to Tables II and I1I—continued 


Possibly connected with Vela-Puppis. Possible concentration in 15. 
Loop, possibly connected with Vela-Puppis. 
Circular region. 
Large area with one bright edge. 
Surrounded by diffuse emission; could be associated with 19. 
Part of Vela-Puppis ? 
Region showing structure. 
Two small regions 1’ x 1’. 
Connected with Vela-Puppis ? 
Circular region. 
Almost circular, with more intense area ‘“‘ comma ’’-shaped. 
Group of 3 bright regions—‘‘ S ’’-shaped. 
Filamentary. 
Emission region containing four concentrations at: 08h 57™-5, — 47° 16’ (6’ x 45 o8h 57m-7 
-47° 22’ (3’ X 2’); 08h 58™-0, — 47° 08’ (9’ x 7’); and o8h 58™-4, —47° 20 *(12’ x 9’). 
Bright region possibly associated with diffuse surrounding background emission. 
Circular. 
Stellar-like. 
Bright crescent shape. 
Emission region outlying 7 Carina nebula. 
Possibly associated with 49. 
Bright knot possibly separated from main » Carina nebula by obscuration. 
Bright knot in » Carina nebula separated from main region by absorption. 
Main irregular 7 Carina nebula. Boundaries at 1oh 28™, —58°; 10h 55M, —58° 18’; 
roh 27M, — 60° 48’. 
Outlying 7 Carina, with concentrations at: 10" 57™-5, —61° 00 (40’ x 20’); 1oh 58™-o, 
59° 20’ (70’ x 40’) (filamentary); 11h 08™-0, — 59° 54’ (25’ x 25’); 11h 10™-5, —58° 30’ 
(25’ x 10’). 
Circular region, appears to have central star; not part of » Carina. 
Appears obscuration-bounded and contains bright crescent-shape.' region 50’ X 20’. 
3° diameter loop. 
Outlying part of A Cen nebula. Contains bright regions centred at: 11% 26-5, —62° 21 
(21’ x 18’); 11h 26™-5, —62° 42’ (15’ x 10’). 
Circular region near A Cen nebula. 
Main A Cen nebula, uneven intensity. 
Large band of emission. 
Region has diffuse edges. 
Crescent-shaped. 
Diffuse edges. 
Planetary ? 
Brighter central region (10’ x 6’). 
Circular. 
Possibly obscuration bound. 
Crescent-shaped. 
Possibly obscuration bound. 
Like a planetary in appearance. 
Circular and filamentary. 
Appears to be surrounding a star. 
Appears to be distributed about a faint central star. 
Stellar-like. 
Stellar-like. 
Contains two bright areas: 16" 15™-5, —50° 51’ (20’ x 7’); 16 17™-0, —50° 42’ (12’ K 12’). 
Bright region near 108. 
Irregular intensity variation with bright region at 16! 36™, — 48° 30’ (60’ x 60’). 
Stellar-like. 
Could be connected to 113 by faint emission. 
Near 110, possibly connected. 
Large loop of ionization in region of fainter emission. Boundaries: (164 37™, — 41°50’) to 
(16h 58™, — 41° 30’); and (16h 50™, — 43° 30’) to (16h 45™, — 39° 40’). 
Ring of emission 50’ wide. 
Concentration inside 113. 
Filamentary, with bright concentration at 16 57™-4, — 38° 13° (15’ x 15’). 
Contains dark central rift. 
Loop. 
Faint extensions to 17h 10™, — 35° 30’. Composed mainly of four bright concentrations. 
Region contains the eae t Sco. 
Appears to be a more centrated region of a large area of diffuse emission. 
Concentrations at: 17! 21™-0, —33° 58’ (30’ x 30 oF Lp — 34° 10’ (9 x9’); 
17h 23™-1, —34° 27’ (6' x6’); 17h 23™-9, —34° 03° (5° = 
Crescent-shaped bright region with faint extensions to 17" 3 “Sm, 26°. 
Almost circular. 
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Notes to Tables II and I1I-—continued 
Circular. 
Stellar in appearance. 
Near galactic centre. 
Diffuse edges. 
Perhaps connected to 141. 
Circular. 
Main bright area divided in three, probably by obscuration, with centres: 18h 01™-8, 
24° 12’ (100’ x 45’); 18h 03-0, —23° 30’ (45’ x 45’); 18h 06-5, —23° 50’ (60’ x 50’) 
Almost circular. 
Stellar in appearance. 
Band of emission, one side bounded by heavy obscuration. 
Stellar in appearance. 
Faint region with bright area [18h 14™-7, —19° 40’ (19’ X 19’)] at one end. 
Two concentrations: 185 08-9, —17° 36’ (7x 7’); 18h og!-3, —17° 40’ (8’x 5’). 
Possibly local concentration of general diffuse emission present. Brighter region at 
8h 14™-9, —17° 18’ (40’ X 20’). 
Faint extensions to 18" o1™, —14° 00 (E 25); and 18h 14™, —17° 00. 
Surrounded by diffuse emission. : 
Surrounded by diffuse emission 40’ wide. 
Oval-shaped, possible concentration of general diffuse emission present. 
Appears connected with adjacent regions by diffuse emission. 
Surrounded by diffuse emission. 
Circular region. 
Crescent-shaped. 
Stellar-like. 


The majority of the nebulae listed in Tables II and III show complex structure 
over large areas. In part this complexity is due to overlying absorption. Thus 
in many cases the identification of individual exciting stars from existing spectral 
data is impossible; an example of this is the extremely faint nebulosity against 
which the Southern Coalsack is silhouetted. ‘The dimensions of this nebulosity 
approximately 7 degrees in diameter, are similar to those of the whole concen- 


tration of B stars forming the I Crucis association, and it is probable that all the 
earlier stars of the association contribute to the excitation of the nebula. 
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THE SPECTRUM OF THE CYGNUS (1gN4A) AND CASSIOPEIA 
(23N5A) RADIO SOURCES BELOW 30MC/S 


A.C. B. Lovell and H. W. Wells 
(Received 1960 January 30) 


Summary 

The 250 ft steerable radio telescope and the 220 ft transit radio telescope 
at Jodrell Bank have been used as an interferometer to measure the ratio of the 
intensities of the Cygnus and Cassiopeia radio sources from 16-0 to 26-0 Mc/s. 
The ratio remains constant at 0°53. ‘The results are compared with previous 
measurements. ‘The amount of absorption to be expected in H 11 regions is 
calculated and it is shown that the difference in H 11 content along the lines of 
sight to Cassiopeia and Cygnus cannot be greater than that contained in a cloud 
of extent 4 parsecs and electron density less than 3 electrons per c.c.; which 
is well below the value normally assumed for an average H 11 cloud. 





1. Introduction.—A recent survey by Whitfield (1) has focused attention on 
the existing discrepancies and difficulties in the determination of the spectra 
of the radio sources. At low frequencies in the 10 to 30 Mc/s region experimental 
difficulties are increased because of the interference from long distance ionospheric 
scatter signals, and the possible existence of differential ionospheric absorption 
introduces another uncertainty into the evaluation of the measurements. These 
uncertainties are manifest in the discrepancies which exist in the published spectral 
information below 30 Mc/s for the Cygnus (1gN4A) and Cassiopeia (23N5A) 
radio sources as given by Hey and Hughes (2), Lamden and Lovell (3) and Wells 
(4). Some of these results indicate that absorption in HII regions might be 
responsible for the introduction of a sharp and differential fall in intensity of the 
two sources below 22 Mc/s. The work described in the present paper does not 
aspire to present new determinations of the absolute intensities of the sources in 
this frequency region. ‘The aim was specifically limited to a measurement of the 
ratio of the intensities of the two sources over the frequency range below 26 Mc/s. 

2. Apparatus.—The aerial systems consisted of the 220 ft transit telescope and 
the 25oft steerable radio telescope at Jodrell Bank. In each case the primary 
feeds were broad-band cage type dipoles and reflectors. ‘These were combined 
in a phase switched interferometer, and since the E-W separation of the instru- 
ments is about 305 metres the lobe interval varied between about 2 and 4 degrees 
over the frequency range used. The receiver consisted of two commercial 
Eddystone Type 680X receivers with a common local oscillator. The separate 
intermediate frequency outputs were fed into an interferometer unit which per- 
formed the phase switching, detection, amplification and other necessary functions. 
By switching at the I.F. and using broad-band preamplifiers coupled to the aerials 
through broad-band transformers, it was possible to change the recording 
frequency by simply adjusting the two receivers and balancing to compensate 
for gain cl.anges. Normally the receivers were used at maximum selectivity 
with crystal filters giving a band width of about 1 Ke/s. 
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3. Results.—The apparatus was used in various forms from February to July 
1959. ‘The successive transits of the Cygnus and Cassiopeia sources through the 
aerial beams required at least four hours and the long distance scatter interference 
made such measurements impossible until around midsummer when the sources 
were in transit during the early morning hours. Even so, very few records were 
sufficiently free from interference or other effects to facilitate the determination 
of the ratio. ‘The successful measurements are summarized in Table I. The 
lowest frequency at which a complete interferometer record of both sources was 
obtained was 16 Mc/s on 1959 July 18-19. This record is reproduced in Fig. r. 
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Fic. 1. The interferometer record (a) of Cygnus and (b) of Cassiopeia obtained on 1959 Fuly 
18-19 on a frequency of 16:0 Mc/s. The abscissae are in sidereal time and the ordinates are 
intensity of signal in arbitrary units. Scintillations are prominent on the Cassiopeia record. 
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The ratio of intensities of the Cygnus and Cassiopeia radio sources from 16 to 26 Mc/s 
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Apart from the interference problems, periods of selective absorption invali- 
dated several of the records. For example, during some 16-7 Mc/s measurements 
on 1959 June 28 the fringes associated with Cygnus failed to appear although the 
Cassiopeia source gave a normal record three hours later. On June 29 the 
Cygnus source was obscured for the first half of the transit but then appeared 
normal, whereas the Cassiopeia source was again unaffected. 

4. Discussion.—The new results for the ratio of Cygnus to Cassiopeia given in 
‘Table I are plotted in Fig. 2, together with the already published data on this ratio 
for frequencies below 40 Mc/s. In addition the ratio assessed by Whitfield (1) for 
the range 30 to 350 Mc/sisshown. Although there is no evidence for the marked 
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3. 2. The ratio of intensities of the Cygnus and Cassiopeia sources as a function of frequency. 


Mean ratic assessed by Whitfield (1) for the range 30 to 350 Me/s. 
4 Hey and Hughes (2) Q Adgie and Smith (6) 
W Lamden and Lovell (3) -] Wells (4) 
A Wells (12) x Lovell and Wells (this paper). 
... Theoretical curve for J N? ds=5 X 10% cm-* 


Theoretical curve for f N? ds=10” cm-5, 


differential absorption below 20 Mc/s which appeared to be indicated by the earlier 
results (3), the bias of the collected results is well below the ratio of 0-62 given by 
Whitfield in the 30 to 350 Mc/s range. According to Whitficld, whereas the 
spectral index of Cassiopeia remains constant at — 0-80, that of Cygnus changes 
from —0-66 in this range to — 1-02 above 350 Mc/s, and a further change in this 
low frequency region below 30 Mc/s cannot be excluded by the present results. 
Since Cassiopeia lies at 2 degrees from the galactic plane and Cygnus at 
5 degrees, this result must be compared with the existing evidence for the pre- 
sence of substantial H 1 absorption (Shain (5), Adgie and Smith (6)). ‘The result 
is not necessarily inconsistent, particularly since Shain (7) has shown that even 
near the galactic centre, there can be very little H 11 absorption at 4 degrees from the 
plane. An upper limit for the differential absorption by H 11 regions in the line of 
sight of Cygnus and Cassiopeia can be calculated from these results. Based on 
earlier work of Smerd and Westfold (8) and of Ryle and Scheuer (g), Lamden 
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and Lovell (3) calculated the ratio of intensities before and after absorption in an 
H 11 region at frequency v as 


I= Iyexp | —0'177 x 10%» * 


where N, is the number of electrons per cm* in the region. ‘Thus 
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The effect of such differential absorption on the ratios is plotted in Fig. 2 for 
values of [N2ds=5x 10cm and 10%cm~° respectively together with the 
experimental points. ‘This indicates that the difference in H 11 content along the 
two lines of sight is such that [N27 ds isabout 107°em~°. Foran HI region with 
an extent of 4 parsecs in the line of sight this implies that the electron density must 
be less than 3cm~*. ‘This is well below the figure for an average H11 cloud of 
20cm~* and dimensions 5 parsecs given by Oort (10) which would imply a value 
of 


” 


| N,2ds=6 x 10%4cm~>. 


This result is consistent with some unpublished measurements of Davies (11) 
in which the 250 ft radio telescope was used to measure the continuous emission 
on frequencies of 400 Mc/s and 1390 Mc/s in the neighbourhood of Cygnus and 
Cassiopeia. ‘The comparison of the brightness temperatures on the two 
frequencies implied that any H 11 region in the line of sight of Cassiopeia gave a 
value of [N,2ds<10x10"%cm~° andof <6x 107°cm~® for Cygnus. 


University of Manchestez, 
Jodrell Bank Experimental Station: 
1960 January 29. 
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INFORMATIONAL CONSIDERATIONS IN THE DESIGN OF 
ASTRONOMICAL SPECTROGRAPHS 


E. H. Linfoot 
(Received 1960 February 15) 


Summary 


The problem of designing spectrographs so as to maximize their informa- 
tion rate on random low-contrast spectral distributions is considered. In 
Section 2 it is shown how spectrographs can be discussed within the general 
framework of Fourier optics. In Section 3 the effects of optical aberrations, of 
photographic spread and granularity, and of variation of slit width are briefly 
considered and the implications of the band-limitedness of optical images in the 
problem of allowing for instrumental bandwidth are pointed out. An 
approximate expression for the information rate is obtained in the special 
case of an aberration-free spectrograph. 





1. Introduction.—The basic problem which confronts the astronomical 
spectroscopist is to make the most efficient use of his time. This includes bot. 
observing time, during which spectra are recorded photographically, and time 
spent on extracting information from the photographic plates. 

‘The advent of very large and expensive telescopes and the increasing use of 
rockets for astronomical spectroscopy have made it of great importance to use the 
observing time as efficiently as possible, even if this should involve the introduction 
of more elaborate reduction techniques. For example, a 10 per cent increase in 
the efficiency with which. information is recorded on the photographic emulsion 
carried by a large rocket spectrograph would justify the purchase of many hours of 
machine time on a microsecond computer to extract this information and recode 
it into assimilable form. 

In the provisional and incomplete discussion which follows, it will be indicated 
how the problem of making the best use of observing time can be approached in a 
more systematic way than has been customary among astronomers in the past. 
We restrict the discussion to the case in which the cost of observing time is so 
large compared with the cost of reduction that the latter can be neglected. This 
limiting case is a specially favourable one for analytical discussion because it can 
be treated in a fairly simple way by an application of Shannon’s information theory. 
Before making the application, we first consider how spectrographs fit into the 
general framework of Fourier optics. ‘This will be examined in the next section. 


2. Notation and preliminaries 

Let S be the plane of an incoherently lit optical object, (x, y) Cartesian co- 
ordinates in S. An optical system images a small region F (the working field) 
near the originin Sontoaplane S’. (x’, y’) are rectangular coordinates in S’, but 
we keep the scale of these coordinates at our disposal. In the later discussion 


(Section 3) S will be the plane of the spectrograph slit, S’ that of the photographic 
receiving surface. 
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2.01 Point-response functions.—A point source situated at (x, y) in S, emitting 
quasi-monochromatic light of wavelength A, isimaged on to S’ by the optical system. 
We can describe the energy distribution in the image by a point-itmage function 
w(x’, y’; x,y; A), normalized so that for all x, y, A 


| [" w dx’ dy'=1. (2.1) 


We regard x’, y’ as the argument-variables of the function w and x, y, A as para- 
meters. The Fourier transform 


— © 
13 x,y; A)=Flw]= {| e-2riqur'+ vy gy dx’ dy’ (2.2) 
- 6) 


will be called the local frequency response of the optical system. It is, of course, 
closely related to the well known ct-function 7(u, v; A) of an isoplanatic system, but 
here the system is not restricted to be isoplanatic. ‘The formal relationship 
between 7 and 7 will become clear later (Sections 2.2, 2.3). We regard the spatial 
frequencies u, v relative to the (x’, y’)-coordinate system as the argument-variables 
of the function 7, and we regard x, y, A as parameters. 

By (2.1) and (2.2) 

7(0, 0; x, y; A)=1 (2.3) 

for all x, y, A. 

Each of the functions w, 7 suffices to describe the imaging by the system of 
incoherently illuminated objects lying in the working field F. By Fourier’s 
inversion theorem, applied to (2.2), 


ex 
w | e2tituc’+ cu) + du dv. (2.4) 
— x 


. 


Thus any property of the optical imaging which can be expressed in terms of w 
can be expressed with equal completeness in terms of 7. 

Evidently 7, like w, depends on the choice of coordinates in S and S’ as well as 
on the imaging process. In particular, we easily see from (2.2) that a shift of 
(x’, y’)-origin to the point x’ = €, y’ = y multiplies 7 by e?"("F*+"™, 

2.1 Gauss optics of a centred system in terms of 7.—In this case we take the 
(x, y)-origin and the (x, y’)-origin on the optic axis of the system and the 
coordinate axes Ox’, Oy’ parallel to Ox, Oy respectively. On the Gauss approxi- 
mation, the image of a luminous object point situated at (x, y) in S is a sharply 
detined point at (x’, v’) in S’, where 


x’ = Mx, y' = My. (2.5) 
If the scale of the rectangular coordinates x’, y’ is the same as that of the 
coordinates (x, y), 1J= M(A)is the magnification at wave-length. Inany case, we 
have 
w(x’, y'; x, v3 A)=8(x' — Mx, y’— My), (2.6) 
where 4 denotes the Dirac 6-function, and 


#(u, Ul; x, Vt r) — e@—2ni(ur+ ry)M(A) (2.7 
In discussing monochromats, it is often worth while to choose the scale of the 
(x’, y’)-coordinates so that M@=1. Then, in the Gauss approximation, the image 
of a bright object point at (x, y) in S is situated at the point x’ =x, y’=y in S’. 
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This gives 
w(x’, y’ 3x, 3A) =8(x’ — x,y’ —y) (2.8) 
7(u, 0; %, y; A) =e Mert OD), (2.9) 
We can describe these (x’, y’)-coordinates as paraxially matched to the (x, y)- 
coordinates. 

2.2 Diffraction optics of a centred isoplanat; relation between 7 and 7.—In a 
centred isoplanat the image of a bright point (x, y) is a light-distribution centred 
on the point (2.5) and radially symmetrical about that point. In other words, the 
point-image function w(x’, y’; x, y; A) has the special form w(x"? + y”*; A), where 

x" =x' — Mx, y"=y'— My (2.10) 
are so-called local coordinates relative to the point (Mx, My) in S’. The factor 
M = M(A) represents the magnification at wave-length A when (x’, vy’) are ordinary 
Cartesian coordinates. In a system of this kind the ct-function at wavelength A, 
represented by the function 7(u, v; A), satisfies the equation 


t(u,v; A) “= 


(° a0 
| e~ Pur’ + ey") gn(x"2 + y"2; A) dx” dy” 
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en 


— e2ri(ux + we 


| e~tekur'+ oy en(x’, y's x,y; A) dx’ dy’, 
~ -@ 


by (2.10), 
= etiurt+ ow (4,0; x,y; A). (2.11) 


From the remark at the end of Section 2.01 we see that r(u, v; A) is the value which 
7(u,v; x,y; A) would assume if the (x’, y’)-origin were shifted to the point which is 
the ‘ideal image’ of (x, y), namely the point x= Mx, y’= My. From (2.7) and 
(2.11) we see that in the approximation represented by Gauss optics 

t(u,v; A)=1 
forall u,v, A. 


2.3 Extension to more general systems.—The last result can be immediately 
extended to the more general case of a system, not necessarily centred, working in 
an isoplanatism patch 4 throughout which the optical distortion can be regarded as 
negligibly small. 

The familiar concept of optical distortion presupposes that some procedure 
has been laid down for specifying which point (x’, y’) in S’ corresponds optically 
toa given point (x, y)in S. In Section 2.2 the symmetry of the spread function w 
provided a basis for such a procedure; we took the centre of the diffraction image 
of (x, y) to be the point of S’ corresponding ‘ideally ’ to (x, y) or, to put it simply, 
matched against (x, y). This procedure is easily seen to be equivalent, in the 
case of symmetrical point-images, with the classical one based on ray tracing. 

In the general case it seems better to depart from the classical procedure and 
to take the mean centre (£, 7) of the ray-theoretic image as the point of S’ matched 
against (x,y). ‘Then (2.11) is replaced by the equation 

7 (u,v; A) =e*hus+ oD Fy, 0; x,y; A) 
= e2milué +e) FT a], (2.12) 
where F'[w] is the Fourier transform (2.2), and where the quantities € = (x, y; A), 
7 =7(x,y; A) depend on x, yand A. 
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A theoretical advantage of the present choice of matching is that with this 
matching the limiting rate of change, near the (u, v)-origin, of the ct-function 
T,(u,v; A) depends only on the aperture of the system and on the wavelength A. 
An inconvenience (which it shares with the classical matching procedure by 
principal rays) is that in general the matching is slightly different in different wave 
lengths A. 

Let 7,(u, v; A) denote the ct-function of a system with the same aperture as the 
given system, but free from aberrations or defocusing. ‘Then it can be shown 


that 
Li. alii 4 (2.13) 
it las | 


Here [¢7/cu] ,, stands for 
lim — 7(u,0; x,¥; A), 
un» CU . 


and similarly for the others. - If the stop is circular, 


OT OT 2A 
=~ =| — -——., 2.14) 
| =| 0 || 0 7 Sin % ose 


where sin z is the numerical aperture of the system*. 
A distortion-free system in the usual sense is one in which, at the so-called 
principal wavelength A = A,, 7 satisfies the condition 


=| (; ; ix) (3 “| Me* (x+iy)—p—1q, 2.15) 
2771 ou ax’ 0 


where the real numbers p, q, the positive factor 17, and the real azimuth & remain 
constant as (x, y) varies in F. 


For, by (2.13) this property of = ensures that 
f= M.(xcos% —ysiny) + 

\ ; : )+P (2.16) 
7 = M.(xsings + ycosy)+q 


and a change of axes in S’ converts (2.16) into the form 


&= Mx, n= My. (2.17) 
In general, colour aberrations of the system may be present, and may show them- 
selves in a small, slow variation of M, us, p,qwithA. In practical cases it is often a 
permissible approximation to assume that (2.15) remains valid for suitably chosen 
M,w, p,q; that is, that the system is distortion-free in each A separately. When 
this is the case, it follows from the definition of isoplanatism that in any isoplanatism 
patch A the values of MW and ¢ can be regarded as independent of A to a sufficient 
approximation. ‘This facilitates the use of a single (x’, y’)-coordinate system for 
the different values of A. 

The usual arguments can now be applied to prove the basic theorem that, in 
quasi-monochromatic light of wavelength A, the intensity distribution in the 
image of an incoherently lit, extended object occupying an isoplanatism patch A 
can be represented, with acceptably small error, as the convolution of the ‘ideal’ 
image of this object with the spread function w,(x",y"; A). The proof can be 


* The last result is due to H. H. Hopkins (1956). 
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extended to systems with the usual small amounts of distortion, but we shall not 
go into further detail here. 


2.4 Spectrographs.—In a spectrograph we choose the (x’, y’)-coordinates 
so that the »’-axis is parallel to the dispersion. Then 7 still satisfies a condition of 
the form (2.15), but now the variation of p with A is no longer small and slow; 
it is smooth and monotonic, and its A-derivative is a measure of the dispersion. 
Conversely, an optical system of which the local frequency response 7 both falls off 
relatively slowly in size as \/42+ 2 increases from zero and satisfies a condition 
(2.15) with the properties just described can be used as a spectrograph; it is 
sufficient to take as object a suitably oriented, incoherently illuminated slit. 

In astronomical spectrographs the illumination of the slit is partially coherent, 
but in some cases we can make the simplifying assumption that the slit may, without 
serious error, be replaced by an incoherent object of the same size and shape. 
This is generally true in prime focus spectroscopy and in wide-slit spectroscopy 
using high dispersion and fast Schmidt cameras. In the case of prime focus 
spectroscopy with an F 4 paraboloid, for example, the radius of the first Airy dark 
ring is about 5A, or about 10°4inch. ‘Thus the effects of coherence on the intensity 
distribution in the image of a slit o-oo1 inch wide may be expected to be small. 
The same could not safely be said of coudé spectrographs working, with narrow 
slits, at focal ratios in the neighbourhood of F/30 and we therefore exclude them 
from the present discussion. 


3. Information in photographic spectra.—I\n an astronomical spectrograph, we 
try to use a slit wide enough to receive nearly all the ‘* wanted ’’ light, a high-resolu- 
tion grating to reduce the overlap of slit images in different wavelengths, a fast 
camera, and an emulsion of suitable speed. 

Even if optical aberrations are left out of the discussion, the choice of camera 
speed requires some thought. Of course, an increase in camera speed reduces the 
time needed to obtain a suitably dense photographic image, but what really matters 
in present circumstances is whether the increase in information rate implied by this 
reduction in exposure time outweighs the reduction in information rate which 
results from the greater relative effect of photographic spread and noise (granu- 
larity) in the smaller picture. Actually, the use of a faster camera of the same 
aperture implies larger optical aberrations, and this means a greater reduction in 
the information content of the photographic spectra due to the interaction 
between aberrations and photographic noise. (See Fellgett and Linfoot 1955, 
p- 371; Black and Linfoot 1958, Fig. 12. In Section 5 of the latter paper, F/2 
should read F/4 throughout; on p. 583 line 3 the values 20 and 5 megabits/cm? 
should read 3-5 and 1 respectively; and on p. 539, bottom line, + 104 should 
read +32u.) We need a means of estimating quantitatively the relative import- 
ance of these two effects on the information rate, in order to be able to make the 
best choice of camera speed for a given slit-width, a given grating and a given 
emulsion. ‘The best choice obviously changes as progress in optical design makes 
new types of system available; for example, the invention of the Schmidt camera 
made the ‘ best’ speed much faster than before because it improved the relation- 
ship between speed and aberrations in available optical systems. ‘The result was 
an increase in the information rate obtainable with a given slit-width, a given 
grating, and a given emulsion. 
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The choice of the photographic emulsion involves similar considerations. A 
faster emulsion means in practice one with larger silver halide grains, similar or 
larger spread, and more granularity. Anincreaseinspread or granularity lowers, in 
a calculable manner, the information content of the photographicimage. ‘Thus both 
exposure time and information gained are less witha faster emulsion. It is their 
ratio which we wish to optimize. If the spread of the ‘aerial’ image of the slit in 
each separate wavelength A is small in comparison with the photographic spread, 
this ratio (the information rate) is only weakly dependent on the speed in a range 
of emulsions of the same equivalent quantum efficiency; that is, the same signal- 
to-noise ratio. Butif the emulsion spread is small compared with the spread of the 
aerial image, then the gain in speed obtained by using a coarser emulsion outweighs 
the loss in information due to the larger spread and granularity; the information- 
rate increases with the emulsion speed if the signal-to-noise ratio in the developed 
emulsion remains constant. 

Quantitative discussion of photographic image quality depends essentially 
on being able to treat the effects of spread and of granularity quantitatively within a 
single analytical framework. Spread can be handled by the use of Fourier trans- 
forms; the discussion of granularity requires the use of probability theory. The 
particular branch of probability theory needed for this purpose is entropy theory or, 
what amounts to the same thing, information theory. 

The discussion of the best choice of slit width is complicated by the question 
of image decoding. From measurements of the photographic density along the 
spectrum we should like to be able to infer the energy density along the ‘ideal’ 
spectrum which the system would give in the absence of diffraction, optical 
aberrations, non-linearity of photographic response, and photographic noise. 
Because an optical system has a finite spatial-frequency pass band, this cannot be 
done. What we can infer is the probability distributions of those coefficients in 
the Fourier transform of the ‘ideal’ spectrum which belong to spatial frequencies 
lying in the pass-band. ‘This is the logically correct method of ‘allowing for 
instrumental band-width’. It requires the use of a fast computer to carry it 
through, however. To testa theory (about line profiles, for example) the logically 
correct procedure is to calculate what values a predicted line profile should give 
to the above-mentioned Fourier coefficients when observed with the system, and 
to compare the calculated values with those obtained by carrying out a Fourier 
transformation on the ‘observed’ intensity distribution obtained by measuring 
the photographic densities along the spectrum and allowing for photographic 
response, which of course varies with A as well as being non-linear. 

Widening the slit leaves the mean values of the above-mentioned inferred 
probability-distributions unaltered, but increases their probability-spreads. ‘Thus 
it decreases the information content. At the same time, it reduces the exposure 
time. Again, we need a means of deciding when their ratio, the information rate, is 
maximized. 

In the remainder of this section we try to indicate briefly the lines on which it 
seems that a useful attempt to deal with these problems could be made. We base 
the attempt on the assumption that the discrimination rate of a spectrograph, here 
defined as the effective number of distinguishable spectra which it can record per 
hour when used to observe a normalized random low-contrast set of spectral 
distributions, is a reasonable indicator of its potential usefulness. The normaliza- 
tion just referred to consists in prescribing for the working wavelength range (A,, Ag) 
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(1) the statistical mean total surface brightness in the slit aperture, (2) the statistical 
r.m.s. deviation of the spectral ‘shape function’ p(A) from its mean value in 
(A, Az). The effective number of distinguishable spectra per hour has for its 
logarithm (to base 2) the statistical mean information rate in bits per hour (compare 
Fellgett and Linfoot 1955). 

We consider a spectrograph with negligible optical aberrations, working over a 
short wavelength range (A,, A.) which, however, is large compared with the resolu- 
tion limit*. The rectangular grating acts as an aperture stop. The slit is taken 
as an incoherently lit object; as already noted, this necessitates excluding some 
cases of practical importance. A simplified 1-dimensional treatment is possible 
in the aberration-free case. 

The geometrical image x’ of an object point x emitting A-light can be taken as 
given by the equation 

A—Ag = a(x’ — Mx) =ax", (3.1) 


where Ay = (A, + A,)/2 and x’ = MxinA,-light. Evidently « is the linear dispersion, 
M the geometrical magnification for fixed A. We can assume photographic 
sensitivity constant over the short A-range, and photographic spread function 
w,(x’; x; A)=w,(x") for all x in the slit, A in the short A-range. In dealing with 
low-contrast images, the photographic process can be treated as linear. 

Then the photographic image is a convolution (o*w*w,*h)(x), where 
ao = a(x") is equal to 1 in the geometrical image of the slit, to o outside it; A=h(x") 
describes the spectral density function of the object geometrically mapped on to 
the photographic surface S’ in accordance with (3.1); w=w(x"; A) is the point- 
image function of the camera with the edge of the grating as aperture stop. 
Evidently 


h(x") =p(A), (3.2) 
where x” and ) are connected by (3.1). If J(x’) denotes the intensity function of 
the image freed from emulsion noise (that is, the statistical mean image-intensity | 
function over a large number of repetitions of the photographic observation of a 
single object), we have by the Fourier product theorem 


FU) = Flo) Fw] Flo) FTA). (3-3) 
Because the wavelength range (Aj, A,) is short, the ct-function F[w] is approxi- 
mately given by the equation 
Fre ba 
=1-Pylu| (Iul< Z-) 


7 ( oe (3-4) 


=0 


where F on the right of (3.4) is the focal ratio of the converging pencils and u 
the spatial frequency in S’. Fora slit —a<x<a of width 2a, 


Flo] = = | ce trius” de" _ sin2ma'u (3.5) 


2ra'u 


where a’ = Ma. 


* A longer wavelength range can be treated by dividing it up into short ranges and considering 
these separately. 
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By arguments similar to those of Fellgett and Linfoot (1955) it now follows 
that the information rate per unit A in the short section (A,, A,) of the spectrum is, to 
a fairly good approximation, proportional to 


UPA 2/sin 27Mau\? 4 
if tor {3 +(1—FA|ul) (ia) | 7(u) | xt du, (3.6) 
where / denotes the emulsion speed (measured in some arbitrarily fixed unit), M 
the geometrical magnification of the slit image on the emulsion, 7, = F[w,] the 
acceptance factor of the emulsion near the wavelength Ay, N the emulsion noise 
power level, P the statistical mean structure power in the low-contrast spectral 
shape-function p(A). The constant of proportionality is independent of the slit 
width 2a, of the linear dispersion «, of the slit image magnification M, of the 
camera focal ratio F, of the emulsion, and of the mean spectral contrast in the 
object set, provided it is small. (But P/N need not be small.) 

The conclusion that when the object set is spectrally random the information 
rate per unit A is independent of the dispersion is not really paradoxical. It corre- 
sponds to the fact that, with a spectrally random object set, the information 
density per unit x’ at a given mean density level is independent of the linear 
dispersion, because a linear magnification of a random set does not change its 
statistical structure. Doubling the linear dispersion therefore doubles the 
information content per unit A and doubles the exposure time needed to obtain this 
information. 

The expression (3.6) does not, of course, purport to give a numerical measure 
of the practical usefulness of the system under ordinary observatory conditions. 
But it provides an indication of the changes which should be made in the slit-width, 
the camera speed, and the emulsion type in order to convert an orthodox spectro- 
graph with very small aberrations into one designed to work in the special con- 
ditions where as much information as possible has to be recorded ina limited time, 
and where a fast electronic computer is available to extract this information. 
The slit-width, camera speed, and emulsion type should be varied within their 
disposable ranges so as to maximize the value of the expression (3.6) for a value of 
P near to that expected for the set of objects to be examined. 

An extension of the analysis to spectrographs in which appreciable aberra- 
tions are present can be based on the formalism in Section 2. In this more general 
case the expression (/— FA, |u|)? in the integrand of (3.6) is replaced by | 7(u) |?, 
where 7(u) is the local optical contrast transfer function of the system at the wave- 
length Ap, evaluated for spatial frequencies along the length of the spectrum. 
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